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INTRODUCTION 


The  study  of  vakes  behind  objects  is  a  challenging  problem  in  fluid 
mechanics.  Even  for  the  least  congilicated  cases  involving  incompressible 
flow  past  simple  body  shapes,  there  does  not  exist  a  unified  theory  that 
provides  a  satisfactory  solution  to  the  problem.  This  is  because  wake 
flow  may  be  considered  as  a  Juxtaposition  of  a  number  of  flow  phenomena  - 
l.e.,  free  shear  layers,  turbulence,  and  flow  separation  —  that  have  yet 
to  be  \inderstood  properly  on  an  individual  basis.  In  view  of  the  intrac¬ 
table  nature  of  the  wake  pi’Oblera,  it  is  not  surprising  that  'until  recently 
research  in  this  area  was  somewhat  fragmentary. 

_  The  advent  of  ballistic  missiles  and  other  reentry  vehicles,  however, 
has  provided  the  impetus  for  a  concerted  effort  to  study  the  wake  problem, 
particularly  in  the  hypervelocity  flow  regime.  The  motivation  lies  in 
the  need  to  monitor  reentry  vehicles  through  the  observation  of  their 
trail,  that  is,  to  relate  the  dimensional  characteristics  of  the  vehicle 
to  the  obsei^ble  quantities  in  the  wake,  such  as  length  of  trail  or 
radar  cross  section.  Research  activities  along  these  lines  have  therefore 
increased  markedly  during  the  past  ten  years.  This  is  evidenced  by  the 
fact  that  in  a  survey  article  on  recent  work,  Lykoudis  (1966)  cites  well 
over  a  hundred  references.  Yet  despite  these  efforts,  there  still  re¬ 
mains  a  lack  of  thorough  understanding  of  the  basic  flow  mechanisms  within 
the  wake.  This  state  of  affairs  attests  to  the  partic’olarly  congjlex  na- 
tiire  of  the  problem  in  the  hypeivelocity  regime. 

To  appreciate  the  essential  features  of  the  hypervelocity  wake, 
consider  the  flow  field  for  a  sphere  as  sketched  in  Fig.  1.  The  regions 
that  dominate  the  picture  are  •'  ’le  recirculation  zone  and  the  viscous  wake 
downstream  of  it.  The  origins  of  the  viscous  wake  can  be  traced  to  the 
boundary  le^rer  over  the  front  half  (approximately)  of  the  sphere.  Near 
the  merldlcui  9  =  «/2  (the  angle  9  defined  in  Pig.  l),  the  bo\mdary 
layer  separates  to  form  an  annxjlar,  free  shear  layer  which  subsequently 
coalesces  near  the  wake  neck  and  becomes  the  viscous  wake.  The  flow 
within  the  viscous  wake  is  Initially  laminar  but  becomes  turbulent 
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somevhere  downstream.  The  position  at  which  transition  occurs  is 
dependent  on  the  Reynolds  number.  Outside  the  viscous  wake,  the  flow 
field  ccMii’onn?  ebscuLially  to  the  inviccid  field  aro’jied  p  >'emlsphere- 
cylinder  model. 

In  the  foregoing  model  of  the  wake  flow,  seveDral  distinctive  problem 
areas  can  be  discerned;  (l)  the  separation  of  the  boundary  layer  and  its 
subsequent  role  as  a  free  shear  layer;  (2)  the  complex  flow  structure  of 
the  "dead-water"  region  or  recirculation  zone;  (5)  the  viscous-inviscid 
interaction  of  the  free  shear  layer  with  the  outer  flow  as  veil  as  viscous 
interaction  with  the  recirculation  zone;  and  (!<■)  the  problem  of  laminar- 
to-turbulent  transition  and  the  description  of  the  transport  mechanisms 
In  the  viscous  wake  when  t\irbulence  sets  in.  In  all  of  these  areas  the 
effects  of  conpresslblllty  and  of  heat  transfer  to  the  body  must  be  taken 
into  account  in  the  hypervelocity  case.  The  theoretician  is  thus  con¬ 
fronted  with  many  complex  problems  in  his  attempt  at  finding  wake  solutions. 
At  the  experimental  level,  the  dlfflcxilty  with  which  data  can  be  obtained 
under  hypervelocity  flow  conditions  makes  it  hard  to  achieve  a  compre¬ 
hensive  verification  of  theoretical  solutions. 

In  an  attempt  five  years  ago  to  £^ln  insight  into  the  wake  problem, 
an  experimental  program  was  initiated  at  Stanford  University  to  study  the 
wake  flow  behind  a  clrcvilar  cylinder  in  a  spark-heated  wind  tunnel.  The 
original  program  called  for  mapping  the  velocity  field  in  the  wake  with 
"tracer-spark"  equipment  as  well  as  for  Pitot-pressure  surveys  with  con¬ 
ventional  methods.  The  tracer-spark  technique  for  velocity  measurement 
was  first  used  for  low-speed  flows  by  Bomelburg  (1958);  an  extension  of 
the  method  to  the  hypervelocity  flow  regime  was  made  by  Kyser  (1964).  In 
addition  to  the  velocity  and  Pitot-pressure  measurements,  exploratory 
studies  were  also  planned  to  determine  the  feasibility  of  using  the  phe¬ 
nomenon  of  resonance  scattering  to  detect  lamlnar-to-turbulent  transition 
in  the  viscovus  wake.  To  achieve  resonance  scattering  in  the  nitrogen 
stream  of  the  tiuinel,  a  nitrogen  laser  was  to  be  used  as  the  light  source. 


The  presence  of  the  recompression  shock  wave  modifies  somewhat  the 
invlscld  flow  field. 
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The  exploratory  studies  shovjed,  however,  that  the  laser  originally 

intended  for  the  resonance-scattering  experiments  did  not  have  sufficient 

pcjwci  uuLput  for  uac  with  Cvjwiiier^ially  available  photorecording  fi?.ms  or 

plates.  The  transition  studies  were  therefore  abandoned.  As  to  the 

mapping  of  the  velocity  field,  it  was  found  that  the  gradients  in  the  wake 

* 

were  too  severe  for  the  tracer-spark  technique  to  apply.  As  a  result  of 
these  instrumentation  limitations,  the  wake-studies  program  was  reduced 
essentially  to  the  surveying  of  Pitot  pressure. 

The  measurement  of  wake  Pitot  pressure  in  itself  is  a  fairly  routine 
"idertaking.  Detailed  surveys  have  been  made  by  McCarthy  (1962),  for 
example,  for  the  case  of  a  circular  cylinder  at  a  Mach  number  of  5*8*  In 
McCarthy's  work,  heat  transfer  was  negligible  since  the  tests  were  con¬ 
ducted  in  a  steady-flow  tunnel.  Smith,  Kramer,  and  Brown  (1966),  on  the 
other  hand,  measured  the  Pitot  pressure  behind  a  cylinder  under  conditions 
where  heat-transfer  effects  were  in^iortant.  In  this  Instance,  the  tests 
were  performed  in  a  short -duration  tunnel  at  hypervelocity  conditions 
(free-stream  Mach  nianber  approximately  20 ),  and  the  data  were  obtained 
with  a  conventional  rake  fixed  with  respect  to  the  tunnel. 

The  experiments  to  be  reported  herein  were  performed  at  flow  con- 
dltlcas  similar  to  those  of  Smith  et  al.  Differences  exist,  however,  in 
the  type  of  model  used  and  in  the  technique  employed  to  survey  the  Pitot 
pressure.  In  place  of  the  conventio'.al  stationary  rake  typical  of  short- 
duration  testing  facilities,  use  was  made  of  rapidly  moving  spring- 
driven  probes  capable  of  recording  continuous  pressiure  data  along  a 
traversed  path  within  the  wake.  This  was  done  to  increase  the  data- 
gathering  capability  of  the  available  pressure-recording  system  and  to 
improve  the  quality  of  data  that  could  have  been  achieved  with  stationary 
rakes.  To  cover  the  entire  wake  region,  two  traversing  probes  were  \ised, 
one  for  axial  traverses  and  the  other  for  traverses  across  the  wake. 
Because  of  the  short  tunnel  running  times,  these  probes  were  made  to  be 


Research  efforts  on  this  technique  were  then  directed  toward  studying 
the  details  of  the  sparking  process  itself.  This  was  done  to  estab¬ 
lish  the  measuring  accuracy  that  can  be  expected  when  the  technique  is 
used  for  measuring  free-stream  velocity  (see  Kyser,  1967)* 
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capable  of  high  speeds  of  traverse  (35  fps  and  60  fps,  respectively).  In 
addition  to  the  rapidly  driven  probes,  a  second  innovation  was  the  imple¬ 
mentation  of  a  technique  that  permitted  testing  in  the  wake  of  a  support- 
free  sphere.  A  description  of  these  systems  and  the  associated  apparatus, 
together  with  a  more  detailed  discussion  of  the  motivation  for  the  use  of 
traversing  probes,  is  given  in  Chapter  2. 

The  use  of  rapidly  traversing  probes,  though  facilitating  the 
gathering  of  data.  Introduced  an  Instrumentation  problem  of  unexpectedly 
serious  proportions.  After  a  considerable  portion  of  the  testing  had 
been  con^leted,  it  became  apparent  that  the  pressure-recording  system  used 
in  conjunction  with  the  traversing  probes  did  not  have  a  fast  enough  re¬ 
sponse  to  cope  with  the  rapidly  changing  pressure  signal  encountered 
during  a  traverse  across  the  wake.  As  a  result,  the  recorded  pressure 
output  exhibited  considerable  time  lag  and  distortion  of  shape.  In  an 
attempt  to  resolve  this  problem,  a  combined  experimental  and  analytical 
study  was  made  to  ascertain  the  response  characteristics  of  the  system, 
comprised  of  a  pressure  transducer,  carrier  amplifier,  and  galvanometer. 
With  the  r  jsponse  characteristics  known,  it  is  possible  to  use  the  solution 
of  an  analogous  problem  in  the  theory  of  sampled-data  control  systems  to 
devise  a  scheme  for  correcting  the  galvar.ometer  output  for  response  lag. 

A  discussion  of  the  response  problem  and  an  outline  of  the  subsequent 
correction  scheme  is  given  in  Chapter  5. 

A  second  difficulty  in  the  experiments,  also  recognized  only  after 
most  of  the  testing  had  been  performed,  was  the  seriousness  of  the  tunnel 
soiarce-flow  effect  on  the  test  data.  It  became  apparent  that  because  of 
the  large  downstream  distances  Involved  in  wake  testing,  the  expanding 
stream  in  the  conical  nozzle  could  significantly  affect  both  the  magnitude 
of  the  Pitot  pressure  and  the  rate  of  growth  of  the  wake.  The  suspected 
soiarce-flow  effect  on  the  magnitude  of  the  Pitot  pressure  has  been  veri¬ 
fied  (in  the  Invlscld  portion  of  the  wake)  by  the  use  of  results  from  the 
method  of  characteristics,  and  this  is  explained  in  Chapter  4.  iEly  means 
of  a  comiblnatlon  of  available  source-flow  solutions.  Including  some  method- 
of -characteristics  results,  an  approximate  scheme  is  then  devised  to 
correct  both  the  magnitude  of  the  measured  Pitot  pressure  and  the  lateral 
wake  dlmensicm. 
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The  salient  points  in  the  final  overall  data-reduction  procedure 
are  discussed  in  Chapter  These  Include  the  essential  steps  in  the 
?ppllcetl''n  of  the  norTectl on  Bchemes  —  for  both  response  lag  and  source- 
flow  effects  —  to  th(  experimental  data.  Representative  Pitot-pressxare 
results,  after  appropriate  correction,  are  also  presented  and  discussed. 
The  corrections  are  quite  large  in  seme  cases  and  lead  to  significant 
changes  in  the  shape  of  the  Pitot-pressure  profiles.  Because  of  the  lack 
of  data  that  are  free  of  errors  due  to  response  lag  and  source-flow 
effects,  it  is  not  possible  to  make  a  direct  assessment  of  the  accuracy 
of  the  corrections.  The  scheme  for  correcting  response  lag  is  checked, 
however,  by  applying  it  to  a  constructed  test  case.  Also,  the  magnitude 
of  the  source-flow  correction  compares  favorably  with  the  results  from 
the  method  of  characteristics.  On  the  basis  of  these  indirect,  partial 
checks,  it  is  believed  that  the  corrections  represent  an  improvement  on 
the  raw  data. 


2.  APPARATUS  AND  TEST  CONDITIONS 


2.1  Wind  Tunnel  and  Model 

The  experiments  vere  performed  in  the  Stanl'ord  University  spark- 
heated  vind  tunnel;  shown  schematically  in  Fig.  2.  This  tunnel,  vhich 
uses  nitrogen  as  the  test  gas,  has  a  useful  run  time  of  about  30  milli¬ 
seconds.  The  nominal  test-section  Mach  nximber  is  17.  It  is  possible  to 
vary  free-stream  conditions  such  as  density  and  temperature  by  appropriate 
choice  of  throat  size,  initial  gas  density  in  the  arc  chamber,  and  amount 
of  energy  imparted  to  the  gas  through  the  arc  discharge.  One  can,  for 
exan^le,  achieve  a  tenfold  variation  in  the  Reynolds  number.  A  discussion 
of  the  range  and  vaiiatiun  of  Reynolds  number  and  stagnation  temperature 
encountered  in  the  experiments  is  given  in  Section  2.7  under  "Test 
Conditions."  A  detailed  description  of  the  tunnel  and  its  operation 
appears  in  the  vork  of  Karamchetl,  Vali,  and  Vincent!  (1961). 

The  model  is  a  2-inch-diameter  sphere,  a  commercially  available  steel 
ball  bearing.  Its  location  vith  respect  to  the  tiannel  differs  according 
to  vhich  portion  of  the  vake  is  being  studied.  The  alternative  of  fixing 
the  sphere  location  and  changing  the  position  of  the  probe  is  not  fea¬ 
sible,  since  the  probe  units  ere  heavy  and  bulky,  and  once  they  are  in¬ 
stalled  in  the  tunr.:.!,  th^y  remain  in  their  specified  position  of  instal¬ 
lation.  In  the  actual  arrangement,  the  sphere  is  placed  in  the  portion 
of  the  nozzle  upstream  of  the  nominal  test  section,  vhich  is  in  fact  a 
portion  of  the  nozzle.  In  the  course  of  the  tests,  eight  different 
stations  vere  used.  These  are  illustrated  in  Fig.  5,  vhich  also  contains 
a  schedule  of  their  distance  from  the  vertically  traversing  probe. 

Since  the  flow  in  the  nozzle  is  continuously  expanding,  free-stream 
conditions  at  the  different  stations  are  not  the  same,  even  for  identical 
conditions  in  the  test  section.  A  more  ccxnplete  discussion  of  the  lon¬ 
gitudinal  variation  in  stream  properties  is  given  in  Section  4.2. 

2.2  Model  Suspension  and  Release 

The  task  cf  obtaining  vake  data  behind  three-dimensional  models  in 
a  vind  tunnel  is  inlierently  difficult.  This  results  from  the  fact  that 
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one  cannot  support  the  model  from  the  rear,  as  is  usually  done  when 
measuring  some  quantity,  such  as  surface  pressure,  on  the  model  itself. 
Several  attempts  have  been  niade  at  various  ways  of  supporting  three- 
dimensional  models  by  means  of  thin  wires  or  rods  (see,  for  example, 
Demetriades  and  Bauer,  I966).  None  of  these  methods  has  been  proven  to 
be  unquestionably  reliable,  however.  A  different  approach  is  used  in  the 
work  of  Vas,  Murman,  and  Bogdonoff  (1965).  In  this  instance  the  model, 
e  sphere,  is  magnetically  suspended  in  hha  st.rftam,  so  that  the  resulting 
wake  flow  is  genuinely  free  of  support  interference. 

The  present  study  is  similar  to  the  work  of  Vas  et  al.  in  that  a 
three-dimensional  model  is  used.  There  is,  however,  a  major  difference 
in  the  method  of  suspension.  In  conti*ast  to  the  steady-state  wind  tiuinel 
used  by  Vas  et  al.,  the  tunnel  en^jlcyed  for  the  present  work  has  the 
relatively  short  running  time  of  30  milliseconds.  In  this  case,  if  the 
model  is  suddenly  released  (i.e.,  left  unsupported)  in  the  tunnel  at  the 
start  of  a  run,  it  will  not  fall  appreciably  from  the  effects  of  gravity 
over  the  duration  of  the  run.  Even  if  the  vertical  drop  of  the  model  is 
not  con^iletely  negligible,  it  can  be  calculated  from  the  usual  time- 
distance  relation  for  free  fall.  The  problem  of  testing  behind  a  support- 
free  model  is  thus  merely  one  of  devising  a  means  of  releasing  the  model 
Just  before  flow  is  established  in  the  tunnel. 

A  system  capable  of  this  type  of  release  for  the  special  case  of  a 
sphere  is  shown  in  Fig.  4.  Before  the  run,  the  sphere  is  suspended  by  a 
wire  spotwelded  to  the  sphere.  At  the  start  of  the  run,  the  wire  is  given 
a  sudden  pull,  whereupon  the  weld  breaks  and  the  sphere  is  released  into 
the  stream.  If  the  weld  is  weak  enough,  and  if  the  impulse  supplied  by 
the  pull  is  sufficiently  large  to  break  the  weld  instantaneously,  only  a 
small  fraction  of  the  energy  (of  the  pulling  motion)  will  be  transmitted 
to  the  sphere.  The  little  that  is  transmitted  does  not  raise  the  sphere 
appreciably,  by  virtue  of  the  sphere’s  large  inertial  mass  (about  1.2  Ibm. ). 
An  experiment  was  performed  to  verify  the  fact  that  the  rise  is  negligible. 
The  sphere  was  suspended  between  a  light  source  and  a  photocell  so  that 
any  movement  of  the  sphere  affected  the  amount  of  light  entering  the 
photocell.  It  was  thus  possible  to  relate  the  photocell  output  to  the 
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verticle  position  of  the  sphere,  rfith  this  system,  it  vas  found  that  the 
rise  of  the  sphere  vas  indeed  negligible. 

Considerable  effort  vas  required  to  find  a  spotvelding  process  that 
provides  a  veak  enough  veld.  The  veld  must  also,  of  course,  be  strong 
enough  to  support  the  veight  of  the  sphere.  An  inqsortant  factor  in 
controlling  the  spotveld  vas  found  to  be  the  material  of  the  vire.  Another 
question  to  be  considered  in  the  choice  of  the  vire  is  the  tensile  strength 
of  the  vire  relative  to  the  strength  of  the  veld.  Although  it  is  desii’able 
to  use  a  thin  vire  to  achieve  a  veak  veld,  an  excessively  thin  vire  tends 
to  stretch  and  break  vhen  the  pxill  is  app'’ied.  After  some  experimentation, 
it  vas  found  that  the  best  results  vere  obtained  vith  a  relatively  thick 
tinned  copper  vire  (gage  no.  22).  In  subsequent  tests,  the  spotvelding 
process  vas  carefully  controlled. 

The  mechanism  for  providing  the  sudden  pull  on  the  suspension  vire 
is  also  shown  in  Fig.  4.  A  heavy  mass  is  set  in  motion  by  the  energy 
stored  in  the  con5)ressed  spring,  and  the  momentum  of  this  moving  mass  is 
used  to  supply  the  pull.  The  mass  is  initially  restrained  from  motion  by 
a  locked  cam,  vhich  is  released  tlirough  the  action  of  a  rotary  solenoid. 

The  triggering  assembly  is  similar  to  the  one  in  the  verti:ally  traver¬ 
sing  probe  and  is  described  in  detail  in  Section  2.4.  A  predetermined 
amount  of  slack  is  provided  in  the  music  vire  that  connects  the  mass  to 
the  suspension  vire.  The  mass  is  thus  permitted  to  acquire  a  pre¬ 
determined  speed  before  it  is  made  to  break  the  veld  and  vithdrav  the 
suspension  vire.  This  speed  is  necessary  so  that  the  vire  vill  clear  the 
flow  field  rapidly.  In  addition,  the  moving  mass  then  exerts  a  large 
enough  impulsive  force  on  the  veld  to  break  it  instantaneously. 

As  the  mass  travels  along  the  guide  tube,  it  trips  two  microswitches 
that  are  installed  inside  a  slot  cut  along  one  side  of  the  tube  (see  Fig. 
4).  The  signal  from  the  first  switch  indicates  the  instant  at  vhich  the 
spotweld  is  broken,  thereby  supplying  the  information  needed  to  calcvilate 
the  free  fall  of  the  sphere.  The  second  switch  is  located  three  inches 
from  the  first  and  thus  indicates  the  instant  at  vhich  the  tip  of  the 
suspension  vire  is  three  inches  from  the  point  of  the  veld  (or  four  inches 
from  the  tunnel  axis,  since  the  radius  of  the  sphere  is  one  inch).  At 


this  instant,  the  tip  of  the  vire  has  cleared  the  inviscic  core  flow  of 
the  tunnel  —  the  inviscid  core  flow  having  a  maximian  radius  of  about  four 
inches  (see  Section  4.2)  —  and  the  wire  is  assumed  to  have  negligible 
effect  on  the  flow  field  around  the  sphere.  Ey  means  of  the  second  switch, 
it  is  thus  possible  to  ensure  that  the  pressure  probing  is  done  only 
after  a  clean  flow  field  has  been  established. 

2.3  Traversing  Pressure  Probes 

With  a  wind  tunnel  of  the  steady-flow  type,  it  is  common  practice 
to  make  traverses  with  various  tjpes  of  probes  within  a  given  region  of 
interest.  For  a  short-duration  tunnel,  however,  conventional  motor-driven 
probes  are  too  slow.  To  survey  the  Pitot  pressure,  for  example,  one  may 
then  resort  to  the  use  of  a  stationary  rake  with  an  array  of  impact  tubes. 
An  exanple  of  this  procedure  is  the  work  of  Smith,  Kramer,  and  Brown 
(1966),  where  a  stationary  rake  was  vised  to  survey  the  Pitot  pressure  in 
the  wake  of  a  circular  cylinder,  aJiso  in  a  spark-heated  tunnel.  The  dis¬ 
advantage  of  such  a  fixed  rake  is  that  a  large  amount  of  instrumentation 
must  be  available.  A  still  more  serious  difficulty  lies  in  the  fact  that 
the  Inpact  tubes  cannot  be  too  closely  spaced,  or  shock  interaction  bet¬ 
ween  adjacent  tubes  will  result.  Because  of  this  limitation,  it  is 
difficult  to  use  data  from  a  given  run  to  resolve  the  detailed  variation 
in  regions  where  gradients  are  high,  as  for  example  in  the  region  of 
transition  from  the  inviscid  to  the  viscous  portions  of  the  wake,  or  in 
the  viscous  wake  itself.  It  is,  of  course,  possible  to  shift  the  rake 
slightly  from  one  run  to  the  next  and  then  combine  the  results  from 
sev;  ’•al  rvins .  Unfortunately,  in  the  case  of  the  spark -heated  tunnel  the 
inherent  lack  of  precise  run-to-run  repeatability  precludes  this  method. 

The  situation  is  worse  when  an  axial  vauriatlon  is  to  be  studied.  In  this 
case,  since  one  cannot  place  inpact  tubes  one  behind  another,  only  one 
data  point  can  be  obtained  with  each  run,  and  poor  repeatability  then 
presents  an  even  more  formidable  problem.  In  view  of  these  difficulties, 
it  was  decided  to  extend  the  traversing-probe  technique  to  cope  with  the 
short  run  time  of  the  spark-heated  tunnel. 

To  estimate  the  necessary  probe  speed,  we  note  that  for  the  probe  to 
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con^jlete  a  traverse  of,  say,  9  Ini.  ss  vithin  the  tunnel  run  time  of  30 
mi.'  iseconds,  the  minimum  probe  speed  must  be,  ideally,  25  fps.  In 
prei-  e  it  is  not  possible  to  coordinate  or  time  everything  precisely, 
and  a  more  realistic  minimum  is  about  35  fps.  Furthermore,  it  is  desirable 
to  make  even  faster  traverses  because  the  flow  in  a  spark-heated  tunnel 
de'  ays  with  time.  To  corap.lete  the  traverse  in  less  than  10  milliseco.ids 
would  require  a  speed  of  about  100  fps.  Such  speed  is  difficult,  if  not 
impossible,  to  attain  with  a  motor-driven  probe.  One  must  resort  instead 
to  pneumatic  systems  or  spring-driven  devices. 

For  the  present  measurements  of  Pitot  pressure,  two  spring-driven 
probes  were  designed  and  built,  (l)  a  vertically  traversing  probe  that 
gives  transverse  profiles  across  the  wake,  and  (2)  an  axially  traversing 
probe  capable  of  surveyint  the  longitudinal  •'variation.  These  are  described 
in  the  follc.jing  sections. 

2.4  Vertically  Traversing  Probe 

General  Description;  The  bulk  of  the  results  to  be  presented  were  obtained 
with  the  vertical li'  traversing  probe,  shown  mounted  on  the  tunnel  test 
section  in  the  schematic  drawing  of  Fig.  5.  This  probe  uses  one  spring 
for  providing  the  driving  force  and  one  for  braking.  The  traversing 
portion  of  the  probe,  sh.  wn  in  Fig.  6,  consists  of  a  length  of  aluminum 
tubing  and  an  aluminum  shell,  rigidly  pin-movnted  at  one  er  5hat  houses 
two  wafer-type  pressure  transducers.  The  tube,  in  addition  to  being  the 
driven  shaft,  also  provides  a  protected  passageway  for  the  transducers' 
electriceil  and  vacuum  connections.  A  steel  collar  pin-mounted  to  the 
tube  at  about  the  half-way  mark  serves  three  functions;  (i,  it  is  the 
"catch"  through  which  the  action  of  the  driver  and  braking  springs  is 
transmitted  to  the  driven  probe;  (2)  it  serves  as  a  guide  for  the  entire 
traversing  assembly  to  slide  vertically  inside  the  outer  brass  tube  that 
encases  the  probe  and  the  springs;  and  (3)  its  lower  flat  surface  engages 
the  locked  cam  when  the  probe  is  being  cocked.  A  set-screw  inserted  in 
the  collar  slides  in  a  vertical  slot  in  the  brass  tube,  thereby  curbing 
rotation  of  the  probe  about  its  own  axis.  The  brass  tube  is  bolted  onto 
a  heavy  base  plate,  on  which  is  also  mounted  the  support  stand  for  the 
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trigger  assembly.  When  installed  in  the  tunnel,  the  bulk  of  the 
traversing-probe  unit  remains  outside  the  tunnel,  the  only  parts  entering 
inside  being  the  transducer-carrier  ho\ising  and  a  portion  of  the  aluminum 
guide  tube.  Since  it  is  difficult  to  make  the  unit  completely  air-tight, 
a  sheet-metal  hood  vith  an  O-ring  attached  is  placed  over  the  unit  during 
tunnel  runs.  The  enclosure  thus  becomes  a  part  of  the  tunnel  vacuum 
chamber. 

Driving  Mechanism;  The  maximum  design  probe  speed  was  set  at  100  fps. 

In  addition,  the  useful  acceleration-free  traverse  of  the  probe  was 
specified  to  be  8  inches. 

To  design  a  driver  spring  tiiat  will  give  these  results,  certain  physical 
limitations  had  to  be  taken  into  accoimt.  For  exninple,  the  spring  was 
to  be  made  from  music  wire,  and  the  largest  music  wire  that  is  available 
commercially  has  a  diameter  of  O.I80  inches.  Furthermore,  at  the  loca¬ 
tion  where  the  traversing-probe  xanit  is  instedlLed,  the  tunnel  inner 
diameter  is  I6  inches.  If  the  probe  is  to  have  an  acceleration-free 
traverse  of  8  inches  centered  about  the  tunnel  axis,  the  probe  and  the 
driver  spring  must  con^jlete  their  acceleration  in  the  first  4  inches  of 
motion.  This  iiig)lles  that  the  initial  congjression  of  the  spring  should 
not  appreciably  exceed  4  inches. 

The  theory  and  formulas  used  to  design  the  spring  are  the  same  as 
those  already  reported  for  the  case  of  the  axially  traversing  probe  (see 
AUigood,  Kyser,  and  Tsao,  1965)  and  will  not  be  elaborated  here.  It 
will  suffice  to  mention  that  the  method  is  based  on  the  surge-wave-theory 
approach  of  Maier,  an  account  of  which  is  given  in  Chironis  (1961).  The 
finished  spring  has  a  coll  outer  diameter  of  1.25  inches  and  is  9*5  inches 
long  when  fully  extended.  It  consists  of  21  turns  of  0.l80-inch-diameter 
mvsic  wire,  first  coiled  into  form  and  tnen  heat  treated  to  increase  the 
yield  point  of  the  steel.  Without  heat  treatment,  the  wire  will  not  be 
stiff  enough  to  give  probe  speeds  close  to  100  fps. 

After  the  spring  s  designed  and  built,  the  congjression  required 
to  give  a  particular  probe  speed  was  found  experimentally  during  bench 
testing.  For  all  subsequent  tunnel  runs,  the  probe  was  operated  at  a 
speed  of  60  fps,  obtained  with  a  spring  compression  of  5*5  inches.  This 


11 


speed  is  short  of  the  design  value  of  100  fps.  This  was  primarily  the 
result  of  the  fact  that  the  mass  of  the  driven  assembly  came  to  exceed 
the  design  figure  when  it  was  found  necessary  to  reinforce  the  cushioning 
and  support  details  ror  the  transducer  wiring. 

Loading  the  Probe;  The  probe  is  loaded  by  first  'engaging  the  steel  collar 
to  the  locking  cam  and  then  compressing  the  spring  by  means  of  the 
threaded  brass  plunger  (refer  to  Fig.  5)*  To  prevent  the  spring  from 
interfering  with  the  probe  after  its  release,  the  spring  is  welded  onto 
the  plunger.  Thus,  when  the  spring  is  being  contpressed,  it  tends  to 
rotaoe  with  the  plimger.  Large  bearing  loads  result  between  the  spring 
and  the  collar,  and  it  was  found  that  of  several  commercially  available 
ball  bearings  tested  none  could  handle  these  loads-  The  difficulty  was 
resolved  through  the  use  of  a  specially  made,  flat  teflon  ring  bearing, 
sketched  in  Fig.  6. 

Trigger  Assembly;  The  mechanism  for  releasing  the  probe  is  illustrated 
in  Fig.  7*  The  cam-and-sear  combination  is  controlled  by  a  Ledex  25“ 
rotary  solenoid  with  an  axial  stroke  of  0.6  inches. 

jt  is  essential  that  the  probe  traverse  be  synchronized  with  the 
tujinel  run.  This  is  accoii5)lished  by  integrating  the  switch  closure  (of 
the  d.  c.  circuit)  into  the  circuit  for  starting  the  tvinnel.  Through  the 
use  of  an  appropriate  nianber  of  relays,  one  can  control  the  timing  to 
within  5  milliseconds.  It  was  found,  however,  that  if  the  probe  were 
left  in  the  loaded  position  for  a  prolonged  time  before  triggering,  the 
various  links  and  joints  would  stick  slightly.  The  delay  time  between 
switch  closure  and  probe  release  can  thus  Increase  by  as  much  as  15  to 
20  milliseconds,  and  even  though  the  switch  is  closed  at  the  specified 
instant  relative  to  tunnel  firing,  the  probe  is  released  too  late  to 
complete  its  traverse  within  the  tunnel  run.  Care  was  therefore  taken 
to  lead  the  probe  just  before  making  the  run.  In  addition,  the  solenoid- 
activating  voltage  was  increased  from  the  rated  50  volts  to  90  volts  to 


These  remarks  concerning  the  necessity  of  synchronization  and  the  means 
of  achieving  it  apply  also  in  the  case  of  the  release  of  the  sphere. 


provide  a  higher  force  from  the  rotary  arm.  When  thes.’  precautions  vere 
taken,  it  vas  possible  to  program  the  probe  traverse  into  the  timnel 
running  sequence  to  within  the  desired  accuracy  of  5  milliseconds. 

Position  Transducer;  The  position  of  the  traversing  probe  relative  to 
the  timne?  is  indicaced  by  a  combination  of  a  light  source  and  photocell. 

The  physical  arrangement  is  shown  schematically  in  Fig.  8(a),  and  Fig. 

8(b)  is  a  diagram  of  the  associated  circuitry.  An  aluminum  ring,  attached 
to  the  brass  tube  near  the  base,  holds  the  light  and  the  photocell  in 
diametrically  opposed  positions.  Small  holes  (l/l6-inch  in  diameter)  are 
drilled  at  half -inch  and  one-inch  intervals  through  the  aluminum  traversing 
tube.  Snug-fitting  Incite  rods  are  inserted  into  these  holes  and  glued  in 
place. 

The  position  transducer  functions  as  follows:  Whenever  one  of  the 
lucite  rods  lines  up  with  the  idealized  light  path,  the  photocell  senses 
the  light  directly.  When  the  wall  of  the  aluminum  tubing  is  in  the  way, 
however,  only  reflected  light  reaches  the  photocell.  The  signal  from  the 
photocell  therefore  exhibits  a  series  of  spikes  corresponding  to  the 
position  of  the  lucite  rods,  and  can  be  used  as  a  position  indicator. 

Figure  9(a)  shows  an  oscillograph  record  of  the  photocell-circuit  output. 

To  establish  the  actual  position  correspondence  of  the  individual  spikes, 
one  such  spike  imist  be  calibrated  by  direct  measurement.  It  is  then 
possible  to  obtain  a  time-position  plot  from  the  o  cillograph  trace.  An 
exanqjle  is  shown  in  Fig.  9(^)>  where  it  can  be  seen  that  the  probe  reaches 
its  peak  speed  after  the  first  3  to  4  inches  of  travel.  This  speed  of 
approximately  60  fps  is  maintained  for  the  8  inches  of  traverse,  within 
which  the  data  recorded  by  the  pressure  transducers  are  considered  reliable. 

Pressure  Transducers  and  Associated  Wiring;  The  probe  carries  two  pressure 
transducers  of  the  variable -reluctance  wafer  design  origir  .ted  at  AEDC 
(see  Smotherman,  i960).  These  are  housed  in  the  aluminum  shell  illus¬ 
trated  in  Fig.  6,  with  the  wafer  plane  verticaJ.  so  as  to  minimize  any 
effects  that  the  probe's  acceleration  might  have  on  the  displacement  of 
the  transducer  diaphragm.  To  obtain  the  fastest  response  possible  with 
the  transducers,  their  pressure-side  inlets  are  used  directly  as  Pitot- 
pressure  orifices. 


Electrical  eind  refersnce-vacuum  connections  to  the  transducers  must 
be  routed  through  the  aluminum  tube.  The  size  of  the  electric  cables 
hence  dictates  not  only  their  own  weight  but  eilso  the  s'^’e  and  weight  of 
the  tube.  It  is  therefore  iTi5)ortant  to  use  the  smallest  cables  available. 

The  present  choice  was  a  Microdot  coaxial  cable  with  an  ou:  side  imeter 
of  0.068  inches. 

Some  difficulties  were  encountered  in  finding  adequate  means  of  anchoring 
the  cables  to  the  shaft.  It  was  thought  at  first  that  potting  the  entire 
transducer-carrier  shell  would  enable  the  soldered  joints  to  withstand 
the  large  in^iact  loads  that  are  involved  during  acceleration  and  braking. 
Potting  eilone  was  not  sufficient,  however,  because  the  cabling  that  hangs 
outside  the  shaft  has  a  tendency  to  whip  about  during  the  in5)ulsive  start 
and  stop  of  the  probe  motion.  The  situation  was  finally  remedied  by 
clan5)ing  the  cables  to  the  shaft  at  the  top,  and  in  addition  the  potting 
was  extended  to  include  the  lower  portion  of  the  shaft.  The  potting 
coirpound  u'^ed  is  a  Genereil  Electric  RTV  102  adhesive.  The  conpound  proved 
to  be  an  ideal  potting  medium  because  it  sets  into  a  rubbery  substance 
that  provides  a  firm  cushioning  svpport;  at  the  same  time  it  does  not 
harden  into  a  solid  that  easily  cracks  upon  inpact. 

The  pressure  transducers  were  patterned  after  the  initial  design  of 
Sraotherman  (i960).  They  are  not  immune  to  external  forces,  as  is  a  later 
acceleration-conpensated  version  of  the  same  basic  transducer.  After 
several  runs  with  the  probe,  it  became  apparent  that  the  probe's  random 
lateral  vibrations  during  descent  had  a  noticeable  effect  on  the  diaphragm 
deflection,  as  manifested  by  a  low-frequency  noise  in  the  transducer  out¬ 
put.  The  noise  level  could  be  high  as  5  to  10^  of  the  actual  pressure 
signal.  To  correct  the  output  for  this  vibrational  noise,  the  upper 
transducer  was  converted  into  an  accelerometer  as  follows:  A  short  length 
of  tygon  tubing,  with  the  front  end  blocked  by  a  thin  meteil  disc,  was 
slipped  over  the  transducer  intake.  A  tiny  puncture  is  made  in  the  wall 
of  the  tygon  tube  to  permit  the  passage  of  air  during  static  ceilibration. 

With  this  arrangement  the  upper  transducer  experiences  no  Pito;  pressure 
during  a  tunnel  run  but  does  register  the  noise  due  to  the  lateral  vi¬ 
bration.  Assuming  that  the  two  transducers  experience  identical  vibration 
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and  that  their  sensitivity  to  acceleration  forces  is  given  hy  static 
calibration,  one  can  then  subtract  the  noise  from  the  ou'.put  of  tl  e  lower 
transducer,  leaving  a  clesin  pressure  signal. 

2.5  A.:ially  Traversing  Probe 

A  detailed  account  of  the  design,  construction,  and  operation  of 
the  axially  traversing  probe  has  already  been  given  by  AUigood  et  al. 
(1965 )•  The  tnecbandcnl  ’iTf'iln-r  to  those  of  the  vertically 

traversing  probe  and  need  not  be  elaborated  upon  here.  An  assembly  draw¬ 
ing  of  the  probe  is  given  in  Fig.  10.  A  photograph  of  the  probe  mounted 
in  the  tunnel  is  shown  in  Fig.  11  (the  nozzle  of  the  tunnel,  including 
the  test  section,  is  rolled  back  to  the  left). 

A  redesigned  Pitot-pressure  rake  has  since  replaced  the  one  shown 
in  the  photograph.  The  details  of  the  new  ralce  are  shown  in  Fig.  12. 

The  spacing  of  the  Pitot  tubes  has  been  modified  to  give  an  improved 
survey  of  the  inner  viscous  wake.  The  size  of  the  tubes  has  also  been 
modified  on  the  basis  of  the  following  considerations  of  low-Reynolds- 
number  effects,  data  resolution,  and  structural  strength.  The  parameter 
controlling  the  effect  of  low  Reynolds  numbers  is  the  ratio  of  the  inner 
to  the  outer  diameter  of  the  tube.  The  effects  are  minimized  by  maxi¬ 
mizing  this  ratio,  that  is,  by  minimizing  the  wall  thiclcness  for  a  given 
outside  diameter.  To  obtain  good  data  resolution,  it  is  desirable  to  use 
as  small  a  tube  as  possible  consistent  with  adequate  strength  from  the 
structural  standpoint.  The  Pitot  tubes  as  finally  selected  have  an  out¬ 
side  diameter  of  0.054  inches  and  a  wall  thickness  of  0.002  inches. 

2.6  Pressure  Instr\anentation 

To  calciilate  the  flow  conditions  in  the  test  region,  measurements 
are  needed  of  the  tiiae  history  of  the  arc -chamber  pressure  and  the  test- 
section  Pitot  pressure  (see  Section  2.7). 

The  arc-chamber  pressure  is  measured  by  a  bonded- strain-gage  trans¬ 
ducer  fabricated  at  Stanford.  The  sensing  element  is  basically  a  four- 
arm  resistive  Wheatstone  bridge.  The  transducer  has  a  range  of  0-20,000 
psig.  It  is  calibrated  periodically  against  a  Heise  gage  that  has  a 
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range  of  0-30,000  psigt  Betveen  such  static  ceilibrations,  the  linear 
behavior  of  the  1  ansducer  is  assumed  to  hold.  The  overall  sensitivity 
of  the  amplifier-galvanometer  system  is  checked  from  run  to  run,  however. 
This  is  done  by  means  of  a  standard  resistor  (in  the  amplifier  circuit) 
that  serves  as  an  extemeil  calibrator. 

The  Pitot  pressure  in  the  test  section,  as  well  as  in  the  wake  flow, 
is  measured  with  a  two-eirm  variable -reluctance  transducer  with  a  range 
of  0-5  psi  differential.  Two  tj^jes  of  gage  are  used;  one  is  of  the 
wafer  design  of  AEDC  but  fabricated  at  Stanford,  and  the  other  is  a 
heavier  Hidyne  unit.  All  Pitot-pressure  transducers  are  calibrated 
statically  before  each  run  against  manometers  containing  either  mercury 
or  oil  of  unit  specific  gravity. 

The  transducers  are  used  in  conjunction  with  20-kilocycle  carrier 
amplifiers  (CEC  Type  1-127).  The  output  is  recorded  on  a  CEC  Type  5-11^ 
recording  oscillograph.  The  recording  galvanometers  are  flat  to  1000  cps 
in  the  case  of  the  arc -chamber  pressure  and  to  120  cps  or  I85  cps 
(galvanometer  option)  in  the  case  of  the  Pitot  pressure. 

The  oscillograph  is  also  used  to  record  the  photocell-circuit  output 
after  it  has  been  an^jlified  by  a  Sanborn  150-1500  d.  .  amplifier.  Ouher 
data  outputs  (viz.,  information  pertaining  to  sphere  release,  etc.)  are 
channeled  directly  to  the  oscillograph. 

2.7  Teat  Conditions 

The  calculation  of  flow  conditions  in  the  test  region  is  based  on  the 
work  of  Grabau,  Hunphrey,  and  Little  (1961)  and  is  i)erformed  on  an  elec¬ 
tronic  computer.  The  range  of  temperatures  in  which  the  calculations  are 
\alid  has  been  extended  with  the  incorporation  of  a  nitrogen-properties 
subroutine  due  to  Smith  (1962). 

The  essential  input  information  required  for  the  calculations  consists 
of  the  arc-chamber  pressure  and  the  stagnation  pressure  at  the  station 
where  the  flow  conditions  are  desired.  The  arc-chamber  pressure  is 
readily  obtained;  the  stagnat''on  pressure  at  the  nose  of  an  unsupported 
sphere  cannot  be  obtained  experimenteilly,  however.  Furthermore,  the 
test  region  is  too  small  to  accommodate  a  separate  probe  to  monitor 


equivalent  infonnation,  i.e.,  the  Pitot  pressure  in  the  free  stream  at 
the  same  longitudinal  station.  It  is  therefore  necessary  to  use  infor¬ 
mation  taJcen  from  a  separate  calibration  run  in  vhich  the  time  history 
of  the  Pitot  pressure  is  obtained  with  a  stationary  probe  positioned  in 
place  of  the  sphere.  It  has  been  found  that  for  a  particular'  set  of 
operating  conditions,  the  ratio  of  arc-chamber  pressure  to  test-section 
Pitot  pressure  is  more  repeatable  than  either  pressure  individually. 
Consequently  instead  of  taking  the  Pitot  pressui-c  directly  from  the 
calibration  run,  ve  adopt  the  ratio  of  the  tvo  pressures  for  use  in  the 
sphere  runs. 

A  consideration  is  now  given  of  the  variation  of  test  conditions 

encountered  in  the  experiments.  The  flow  quantity  that  will  be  used  as 

the  primary  parametric  variable  is  the  Reynolds  number.  Re  ,  based 

in^Q 

on  the  diameter  d  of  the  sphere  and  on  free-stream  conditions  at  the 
nose.  Three  sets  of  operating  conditions  have  been  used  to  obtain  three 
basic  Reynolds  numbers.  Within  each  of  the  three  sets  of  conditions, 
some  variation  in  Reyno-Tls  number  arises  from  the  following  causes;  (l) 
the  tunnel  flow  decays  in  time,  and  the  Reynolds  number  decreases  during 
a  given  probe  traverse;  (2)  free-stream  conditions  at  the  nose  of  the 
sphere  w.'e  different  for  different  positions  of  the  sphere  since  the 
tunnel  stream  expands  along  the  nozzle;  and  (3)  overall  tunnel  flow 
conditions  vary  somewnat  from  run  to  run  because  of  lack  of  exact  repeat¬ 
ability. 

Typical  flow  conditions  are  summarized  in  Tables  2.1  and  2.2  at  the 
end  of  this  section.  The  results  of  Table  2.1  pertain  to  a  f'lxed  test 
station,  the  variable  being  the  operating  conditions.  In  contrast,  the 
results  of  Table  2.2  reflect  the  variation  in  conditions  caused  by  a 
shift  of  the  test  station.  In  both  tables,  the  time  for  which  the  data 
are  given  is  20  milliseconds  sifter  the  start  of  the  t\mnel  run,  which 
corresponds  roughly  to  the  middle  of  the  run. 


Table  2.1 

Free-Stream  Conditions  at  Station  A 


Total  Temperature 

“K 

Mach  Number 

Reynolds  Number 
(based  on  d  =  2") 

5,730 

16.0 

2.l8  X  10^ 

2,570 

17.7 

7.40  X  10^ 

1,'^50 

16.2 

1.77  X  10^ 
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Table  2.2 


Variation  of  Free-Stream  Conditions  vith  Location  of  Test  Station 


Station 

Total  Teit^i. 

“K 

Mach  Number 

Reynolds  Number 

(based  on  d  =  2" 

A 

1,270 

16.31 

2.35  X  10^ 

B 

1,270 

16.16 

2.40  X  10^ 

C 

1,270 

16.00 

2.46  X  10^ 

D 

1,270 

15.84 

2.52  X  105 

E 

1,270 

15.68 

2.58  X  10^ 

F 

1.270 

15.51 

2.64  X  10^ 

G 

1,270 

15.55 

2.70  X  10^ 

H 

1,270 

15.03 

2.81  X  10^ 
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5.  TRANSDUCER-G.AI.VM(M?rSR  RESPONSE  AND  CORRECTION 

5.1  Preliminary  Remarks 

When  pressure  instrumentation  was  first  chosen  for  the  Stanford 
spark-heated  tunnel,  it  was  decided  that  an  oscillograph  would  be  used 
for  recording  pressure  information.  The  oscilloscope,  despite  its  superior 
response  and  versatility  in  other  respects,  was  passed  over  in  deference 
to  the  oscillograph’s  capability  of  recording  as  many  as  ?6  channels  of 
data  on  a  single  time  base.  To  use  oscilloscopes  to  provide  the  equiva¬ 
lent  capacity  would  not  have  been  economically  possible.  Also,  the  task 
of  operating  uhe  many  osciLloscopes  simultaneously  during  calibration  or 
tunnel  run  would  have  been  difficult.  Furthermore,  data  reduction  is 
easier  and  more  accurate  on  the  oscillograph.  These  factors  all  contri¬ 
buted  towards  the  selection  of  the  oscillograph  as  the  basic  instrument 
for  recording  pressure  data. 

3.2  Response  Considerations 

For  the  recording  of  routine  arc-chamber  pressure  and  of  Pitot 

pressure  in  the  test  section  (measured  with  a  stationary  probe),  the 

choice  of  the  oscillograph  was  satisfactory  because  the  different  types 

of  galvanometer  that  were  available  turned  out  to  be  adequate  for  their 

respective  assignments.  For  exanqjle,  in  the  case  of  arc-chamber  pressure, 

where  good  response  is  essential,  one  can  use  a  relatively  fast  galvanometer 

(flat  to  1000  cps).  On  the  other  hand,  for  the  Pitot-pressure  transducer, 

where  a  substantial  level  of  high-frequency  noise  in  the  amplified  output 

requires  a  slower  galvanometer  (flat  to  125  or  l80  cps),  response  is  not 

* 

of  prime  concern  when  only  the  test-section  Pitot  pressure  is  involved. 

When  the  traversing  probes  were  Introduced,  the  same  combination  of 
galvanometer  and  Pitot-pressure  transducer  was  taken  over.  This  led  to 
a  response  problem  that  was  recognized  only  after  a  substantial  amount  of 

* 

This  difference  in  response  requirements  stems  from  the  fact  that  the 
time  history  of  the  arc-chamber  pressure  is  the  more  critical  of  the 
two  pressure  inputs  in  the  calculation  of  test-section  flow  conditions. 


data  had  already  been  taken.  The  problem  Is  particularly  acute  vlth  the 
vertically  traversing  probe,  because  the  probe  speed  Is  high  and  large 
gradients  are  present  In  the  transverse  Pitot-pressure  profiles.  For 
exan^jle,  for  the  probe  speed  of  60  fps  used  and  for  an  ovei'all  system 
response  lag  of  -  say  -  one  millisecond,  the  entire  profile  would  appear 
to  be  shifted  off  center  by  about  eight -tenths  of  an  inch.  In  addition, 
since  large  gradients  are  present  In  the  profile.  Its  shape  would  be 
distorted,  and  the  pressure  trace  actually  recorded  would  lose  some  of 
the  symmetry  normally  to  be  expected  In  such  a  profile.  In  other  appli¬ 
cations,  for  Instance  the  axially  traversing  probe,  response  requirements 
are  not  so  exacting.  The  discussion  that  follows  will  pertain  only  to 
the  vertically  traversing  probe. 

When  It  became  apparent  that  the  response  of  the  pressure-recording 
system  was  not  fast  enough  for  work  with  the  vertically  traversing  probe. 

It  was  also  recognized  that  It  would  be  difficult  to  Improve  on  the  respons? 
by  replacing  con^jonents  of  the  existing  system  with  other  readily  available 
equipment.  One  possibility  was  to  replace  the  galvanometer  with  an  oscil¬ 
loscope.  This  was  not  attenipted,  however,  because  previous  experience 
showed  that  the  anqjllfier  output  exhibited  considerable  high-frequency 
noise,  and  to  obtain  a  clean  pressure  signal  would  have  required  that  a 
low-pass  filter  be  inserted  ahead  of  the  oscilloscope.  From  the  response 
standpoint,  the  filter  would  then  effectively  reduce  the  oscilloscope  to 
the  level  of  a  galvanometer,  thereby  defeating  the  purpose  of  using  the 
oscilloscope  in  the  first  place. 

Besides  the  slow  response  of  the  galvanometer,  there  Is  also  the 
question  of  the  response  behavior  of  the  transducer  Itself.  The  problem 
of  transducer  response  has  been  studied  by  Smotherman  and  Maddox  (1963), 
who  determined  the  rise  time  for  various  configurations  of  the  wafer  trans¬ 
ducer.  This  work  is  not  directly  applicable  to  the  Stanford  version  of 
the  basic  AEDC  design,  however,  since  some  differences  exist  In  the 
dimensions  and  construction  details,  and  these  control  the  response  char- 
acterlstl's  of  the  transducer.  Furthermore,  Smotherman  and  Maddox  did 
not  study  the  response  of  the  combined  transducer-an^jllfler-galvanometer 
system,  which  is  of  interest  her«.  It  was  decided  therefore  that  a 
separate  experimental  study  should  be  conducted  to  obtain  the  desired 
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infoimation.  When  the  response  characteristics  of  the  overall  system 
T^  ere  known,  one  could  hope  to  devise  an  analytical  correction  scheme 
that  could  Chen  be  applied  to  the  recorded  pressure  output. 

The  experimental  determination  of  the  system  response  characteristics 
is  based  on  the  theory  of  linear  systems.  It  is  therefore  appropriate 
to  introduce  some  of  the  fundame«:tal  concepts  before  proceeding  t  )  the 
results  of  the  experimental  study. 


^•3  System  Transfer  Function 

The  transducer-galvanometer  combination  may  be  considered  as  a 
linear  system  provided  an  approximation  is  made  regarding  the  transducer 
itself  (see  Appendix  A).  In  w.'jat  follows,  the  assumption  of  a  linear 
system  will  be  adopted  throughout. 

The  response  characteristics  of  the  transducer-galvanometer  system 
is  specified  if  the  system's  transfer  fvuiction  is  known.  In  the  theory 
of  linear  systems,  a  transfer  function  G(s)  is  defined  (see,  for 
example,  Clark,  1962)  as  the  ratio  of  the  Laplace  transform  of  the  out¬ 
put  to  the  Laplace  transform  of  the  input,  i.e., 


G(s) 


P  4.(s) 

out 

p,  (s; 
in 


*{Pout(‘h 


(3.1) 


where  s  is  the  Laplace  variabl:,  and  and  Pp^^(t)  are  re¬ 

spectively  the  input  and  output  to  the  system,  as  shown  in  Fig.  15(a). 

■;n  the  present  case,  Pj^jj('t)  and  pQy^(t)  refer  specifics lly  to  the 
pressure  input  and  output.  When  the  transfer  function  is  known,  the  system 
o\  .put  to  a  specific  input  is  the:'  gi’.'en  by 


=^'^{G(s)-P^^(s)) 


(5.2) 


A  convenient  way  of  obtaining  the  transfer  function  is  to  consMer 
the  response  of  the  system  to  a  unit  step- input  (whose  Laplace  trensform 
is  l/s),  the  output  in  this  case  being  known  as  the  "transient  response". 
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If  the  Laplace  transform  of  the  trant^ient  response  is  denoted  by  X(s), 
then  from  Eq.  (5.1)  the  transfer  function  is  sin^jly  sfXC?)].  This 
relationship  in  fact  provides  a  means  for  experimental  determination  of 
the  transfer  function.  That  is,  one  can  in5)ose  r  step  function  on  the 
system  and  record  the  corresponding  output;  the  output  can  then  be  com¬ 
pared  vith  theoretically  calculated  transient-response  curves  (given  in 
standard  textbooks  on  automatic  control  systems )  to  yield  the  appropriate 
function  for  X(s),  vhich  ir.  turn  gives  the  transfer  function.  This 
procedure  is  not  difficult  if  the  order  of  the  system  is  known.  For  a 
system  of  con5)onents  in  series,  an  alternative  approach  is  to  find  the 
transfer  function  for  each  congjonent;  the  overall  transfer  function  is 
then  siii5)ly  the  product  of  the  individual  transfer  functions  (see  Clark, 

1962,  p.  50). 

To  find  the  xransfer  function  for  the  transducer-galvanometer  system, 
the  experimental  approach  was  used.  Briefly,  a  step  increase  in  pressure 
of  known  magnitude  was  generated  at  one  end  of  a  simple  shook  tube  by 
means  of  shock -wave  reflection.  The  tran  ducer  positioned  at  this  end 
thus  experienced  ■  ftep  input  of  pressure.  The  transducer  output  and 
galvc nometer  output  were  then  used  to  find  the  transfer  function  of  the 
system  components.  The  experimental  details  and  the  intermediate  steps 
used  in  arriving  at  the  composite  transfer  function  G(s)  for  the 
ti'ansdu^er-galvanometer  system  are  outlined  in  Appendix  A.  The  final 
result  for  G(s)  was  fo\md  to  be 


koj  ,  -s& 

G(s)  .  °  ^ ^  , 

s(s  +  203  ts  +  0)  ) 

'  n^  n' 


(5.5) 


where  o)^  =  undan^jed  natural  frequency  of  galvanometer, 

^  =  damping  ratio  of  galvanometer, 

k  =  correctibn  factor  ■'•o  transducer,  and 
5  =  approximate  rise  time  of  transducer. 

The  values  of  the  system  parametric  constants  o)^,  k,  and  6  are 
also  given  in  Appendix  B. 
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The  determination  of  the  transfer  function  defines  the  response 
characteristics  of  the  transducer-galvanometer  syetesi  It  does  not, 
however,  provide  directly  a  means  for  calculating  the  system  response  to 
a  pressure  input  of  arbitrary  shape,  such  as  that  encountered  in  the  wake. 
This  follows  from  the  fact  that  in  general  an  arbitrarily  varying  time 
function  cannot  be  expressed  analytically,  and  its  Laplace  transform  -  if 
it  exists  at  all  —  cannot  be  calcx'J.ated  in  closed  form.  The  problem  of 
finding  the  output  to  an  input  of  arbitrary  shape  has  an  approximate 
solution,  however,  provided  a  sampling  process  is  introduced  at  the  input 
and  at  the  output.  We  shall  first  show  that  the  data- reduct ion  procedure 
in  the  present  case  does  introduce  a  fictitious  data-saiiipling  process. 

The  resulting  problem  is  then  amenable  to  treatment  by  the  theory  of 
sarapled-data  control  systems.  One  of  the  basic  mathematical  devices  of 
this  theory,  the  z-transform,  is  introdw.oed  in  Section  5.5,  which  also 
gives  an  account  of  the  application  of  the  theory  to  the  present  case  of 
the  pres  sure -recording  system.  After  we  have  establishes',  the  procedure 
for  calculating  the  sampled  output  from  the  input  and  the  transfer  func¬ 
tion,  we  then  consider  the  inverse  problem  of  deducing  the  input  when  the 
output  is  given.  The  solution  is  based  on  a  method  of  successive  approxi¬ 
mations  and  is  presented  in  Section  5*6. 

5.4  Analog  in  Sarapled-Data  Theory 

In  the  data-reduction  scheme,  output-pressure  information  is  obtr  ned 
by  reading  deflections  off  the  oscillograph  record  at  regular  time  inter¬ 
vals,  typically  of  one -half  millisecond  for  runs  with  the  vertically 
traversing  probe.  This  procedxire  may  be  likened  to  sampling  the  pressure 
output  by  inserting  after  the  galvanometer  a  switch  that  closes  momentarily 
every  T  seconds,  T  being  the  sampling  period.  An  equivalent  sampling 
process  would  be  to  block  the  galvanometer  light  beam  from  the  photo- 
recording  paper  at  all  tines  except  for  a  flash  every  T  seconds.  This 
analogy  with  a  sarapled-data  system  is  illustrated  in  Pig.  15.  (For  the 
fundamentals  of  saupled-data  theory,  see  Kuo,  19^5,  or  Pagazzini  and 
Freuiklin,  1958. )  The  physical  transducer-galvanometer  system,  whose 
block  diagram  was  shown  in  Fig.  15(a),  is  a  system  for  continuous  signals. 
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The  switch  that  slmiilates  the  saitqpllng  process^  the  so-called  "ideal 
sampler",  is  depicted  In  Fig.  13(b)  together  with  a  graphical  re  >  3- 
sentatlon  of  Its  sampling  function  In  the  form  of  a  unit-impulse  train. 

The  block  diagram  for  the  analogov  sampled-data  system  Is  shown  In 
Fig.  15(c).  (The  function  of  the  "hold  circuit"  placed  ahead  of  the 
pressure  transducer  will  be  explained  In  the  next  section. ) 

We  note  that  In  addition  to  uhe  Insertion  of  a  sampler  switch  at  the 
output  In  Fig.  13(c)  —  which  Is  the  appropriate  representation  of  the 
process  of  reading  off  discrete  valuas  -  a  similar  switch  Is  Inserted  at 
the  Input  as  well.  This  Input  sampler  Is  Introduced  to  complete  the 
analogy,  since  In  a  san^jled-data  system  the  presence  of  this  sampler  Is 
mandatory.  This  Is  because  the  fvindamental  problem  In  sampled-data  theory 
Is  to  find  the  appropriate  manner  of  sampling  the  Input  so  that  one  needs 
to  feed  to  the  system  only  the  sampled  sequence  Instead  of  the  continuous 
Input  while  still  being  able  to  extract  approximately  the  same  Information 
at  the  output.  Thus  In  a  sampled-data  system  the  Input  sampler  Is  the 
essential  switch;  the  synchronous  switch  at  the  output  Is  Inserted 
simply  to  facilitate  the  analysis.  In  the  analog  for  the  transducer- 
galvanometer  system,  on  the  other  hand,  the  output  switch  Is  the  more 
meanlngfxil  one  because  It  Is  at  the  output  that  discrete  values  are  read 
off;  the  Input  sampler  Is  Included  so  that  we  may  draw  an  exact  analogy 
with  the  sampled-data  problem  and  use  the  available  results  from  the 
theory.  With  the  Introduction  of  the  fictitious  sampler  Into  the  trans¬ 
ducer-galvanometer  system,  however,  the  question  arises  as  to  whether  the 
sequence  of  Input  sampled  by  the  switch  Is  an  accurate  representation  of 
the  continuous  Input  at  the  sampling  Instants.  Thus,  after  the  analogy 
has  been  drawn,  the  subsequent  problem  Is  twofold;  (l)  to  find  the  appro¬ 
priate  sampling  procedure  for  the  transducer-galvanometer  system,  and 
(2)  to  show  that  for  this  system  there  is  Indeed  a  correspondence  between 
the  sampled  Input  and  the  continuous  Input. 

To  find  the  appropriate  sampling  procedure  In  the  present  Instance, 
use  was  made  of  the  test  case  Involving  the  transient  response  of  the  system. 

*  . 

This  test  case  has  been  considered  In  the  last  section,  Eq.  (5.5) 

and  the  paragraph  preceding  It. 


Since  the  transient  resptmse  vas  determined  experimentally,  it  vas 
possible  to  choose  the  san^jling  procedure  that  gave  the  best  agreement 
vith  experimental  results,  and  this  vill  be  discussed  in  the  next  section. 
To  show  that  there  is  correspondence  between  the  sampled  input  and  the 
continuous  input  (the  second  part  of  the  problem),  the  following  plau¬ 
sibility  argument  is  used.  In  sarapled-data  theory,  the  sarapled-sequence 
representation  of  the  input  is  by  definition  exact,  irrespective  of  the 
sar^jling  procedure.  The  basic  tenet  of  the  theory  states  that  with  the 
proper  choice  of  a  saraplirg  process,  the  sampled-output  information  then 
accorately  represents  the  continuous  output.  Now,  in  the  analogous  case 
for  the  transducer-galvanometer  system,  the  roles  of  input  and  output  are 
Interchanged.  At  the  output  the  sampled  sequence  is  known  to  correspond 
to  the  continuous  signal,  since  the  sampling  describes  in  actuality  the 
process  of  reeding  discrete  values  off  a  continuous  trace  on  the  oscil¬ 
lograph  record.  Now,  by  virtue  of  part  (l)  in  the  foregoing,  it  can  be 
assumed  that  the  appropriate  8aii5)llng  procedure  is  ujed.  We  thus  infer, 
by  an  application  in  reverse  of  the  basic  tenet  of  sanqpled-data  theory, 
that  the  saii5)led  input  does  represent  adequately  the  continuous  input 
at  the  particular  sanqpling  instants. 

5.5  The  z -Transform  Method 

In  the  theory  of  sampled-data  control  systems,  the  basic  mathe¬ 
matical  tool  is  the  z-transform.  The  z-transform  is  actually  a  short¬ 
hand  notation  for  the  reciprocal  of  the  Laplace  transform  of  a  unit 
iii5)vilse  occurring  at  t  =  T,  that  is, 

^{6(t-T)}  =  e‘^®  =  z"^  .  (5.4) 

Similarly,  the  Laplace  transform  of  the  unit-impulse  train  of  Fig.  15(b), 
which  is  basic  to  the  data-o.  mpllng  process,  is  expressible  in  terras  of 
a  polynomial  in  z  ^  as  follows: 


X  it  8(t-nT)}  =  f  e-"^^  =  £  •  (5-5) 

n=0  n=0  n=0 
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Since  8  sampled  sequence  (with  sampling  period  T)  of  any  function 

00 

f(t)  can  be  represented  by  f(nT)[&(t-nT)J,  it  follows  from  Eq. 

n=0 

(5.5)  that  the  Laplace  transi jrm  of  the  sampled  input  can  be  expressed  as 


^  Pin(nT)6(t-nT))  =  £  Pin(nT)z""  =  .  (3.6) 


n=0 


The  corresponding  expression  for  the  sampled  output  is 


^out(^^  =  Z  Pout("^^^“”  • 


n=0 


(3.7) 


It  is  shown  in  standard  texts  on  sanpled-data  systems  that  the 
functional  relationship  between  Pj^j^(z)  and  is  completely 

analogous  to  the  corresponding  case  of  Eq.  (3»l)  for  continuous  signals. 
Thus,  for  sanpled-data  systems  the  relationship  is 


in' 


(3.8) 


where  G(z),  the  z-transforra  of  the  transfer  function,  is  eM-.luated 
according  to  the  formula  (see  Ragazzini  and  Franklin,  1958,  p.  57) 


G(z) 


residue 

poles  of 
cKs) 


G(s) 


T  Ts  - 
1-e  z 


(3.9) 


The  function  G(z)  is  normally  expressed  in  terms  of  the  quotient  of  two 
polynomials  in  z“^.  Further  simplification  is  possible  by  long  division 
to  yield  a  single  polynomial  in  z”^. 

Ideally,  Eq.  (3.8)  used  in  conjunction  with  Eqs.  (3.6)  and  (3.?)  for 
and  P(^^(z)  gives  a  method  for  calculating  PQy^(z)  vhen 
G(z)  and  Pj^jj(z)  are  given  (by  equating  coefficients  of  like  powers  in 
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z  ^).  However,  if  the  sampled  input  were  fed  directly  into  the  linear 
system,  some  of  the  original  Information  between  sanpllng  instants  would 
be  lost,  resulting  in  an  Inaccurate  output.  It  is  therefore  necessary 
to  use  some  form  of  lnterp)olation  scheme  to  supply  additional  information 
between  the  sanipllng  Instants.  This  is  done  by  the  fictitious  hold 
circuit,  or  "data  hold".  Inserted  between  the  input  sampler  and  the  trans¬ 
ducer  (refer  to  Fig.  13(c)).  The  function  of  the  hold  circuit  is  to 
reconstruct  the  sampled -input  seq^uence  and  transform  it  into  a  piecewise 
continuous  or  piecewise  smooth  data  curve  before  feeding  the  information 
into  the  linear  system.  Interpolators  of  various  orders  can  be  used, 
depending  on  the  accuracy  desired  and  on  the  shape  of  the  partlc’ular  data 
curve  under  consideration.  A  detailed  discussion  of  the  subject  of  hold 
circuits  is  given  by  Jury  (1964). 

In  the  present  Instance,  the  most  suitable  hold  circuit  is  a  first- 
order  interpolator  having  the  transfer  function  (see  Jury,  1964,  p.  222) 

Q(z,s)  =  .  (3.10) 

Ts 

Because  this  transfer  fvmctlon  contains  the  Laplace  variable  s,  the 
function  must  be  included  in  the  right-hand  side  of  Eq,.  (3*9)  before  the 
z-transform  of  the  overall  transfer  function  can  be  evaluated.  For 
clarity,  the  overall  transfer  function  that  Includes  the  first-order 
interpolator  is  denoted  by  S(z,s),  that  is, 

S(z,s)  =  Q(z,s)  G(s)  ,  (3»ll} 

where  G(s)  for  the  transducer-galvanometer  system  is  given  by  Eq.  (3.3). 
The  algebraic  details  of  evaluating  t>(z),  the  z-transform  of  (j(z,8), 
from  Eq.  (3*9)  are  outlined  in  Appendix  B.  The  final  result  is 


^  -lu.  -2 
m  j^z.+m  +m^z  +mgZ 

-1.  ■  ^2 


8(z)  =4 


(3.12a) 


•  •  • 


(3.12b) 


or. 


= 


-1 


^2" 


-2 


where  T  =  cd^T,  6  =  and  X,  m  nig,  dg  are  all  non- 

dliTjnsional  fimctlons  of  the  system  parameters.  The  series  representation 
of  S(z)  in  Eq.  (3.12b),  found  by  straightforward  division  from  Eq. 
(3.12a),  is  in  a  form  convenient  for  later  use.  The  coefficients  r 
r^,  r^,  r^, . . .  have  been  calculated  with  an  electronic  computer  and  are 
listed  in  Appendix  B.  An  inspection  of  the  magnitudes  suggests  that  the 
series  is  convergent. 

In  the  analogy  of  the  transducer-galvanometer  system  with  a  sanpled- 
data  system,  we  have  assiimed  that  a  suitable  sampling  rate  has  been  used. 
Furthermore,  no  immediate  Justification  was  given  for  choosing  an  inter¬ 
polator  of  a  particular  order  as  ^he  appropriate  data  hold.  These  two 
aspects  of  the  problem  are  interrelated,  and  several  combinations  of 
sampling  period  and  interpolator  were  tested  to  establish  the  most  sat¬ 
isfactory  combination.  The  criterion  used  to  gage  the  relative  success 
of  a  particular  combination  was  to  determine  how  closely  the  calculated 
sampled  output  matches  the  experimentally  recorded  output  corresponding 
to  the  unit-step  pressure  input  generated  by  shock-wave  reflection. 
Interpolators  of  zeroth,  first,  and  second  order  were  tried  and  the 
sampling  period  was  varied  from  0.1  to  0.5  milliseconds.  The  best  results 
were  obtained  with  the  first-order  interpolar  of  Eq.  (3.10)  and  a  sampling 
period  of  0.25  milliseconds.  Now,  j.  .-reduction  procedure  for 

wake  i\ins  (discussed  in  Chapter  5)^  the  time  interval  betweoi:  data  read¬ 
outs  was  0.5  milliseconds.  To  use  the  foregoing  results  with  a  sanpllng 
period  of  0.25  milliseconds,  additional  data  had  to  be  obtained  at  points 
midway  between  any  two  adjacent  data-readout  instants.  This  additional 
information  was  obtained  by  interpolating  the  0.5-milllsecond  data.  The 
interpolation  was  carried  out  with  Lagrange's  interpolation  coefficients 
(see  Henrici,  1964,  pp.  l83-l85);  either  a  three-point  or  a  four-point 
formula  was  used,  depending  on  the  portion  of  the  output  trace  considered. 
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5.6  The  Indirect  Correction  Scheme 

With  the  inclusion  of  the  data  hold  in  the  overall  transfer  function, 
Eq.  (5.8)  must  now  he  modified  to  read 

■  (3.13) 

The  problem  at  hand  is  to  find  when  ^(z)  are 

given.  Conceptually,  one  is  therefore  inclined  to  seek  an  inverse  of  the 
transfer  function  so  as  to  calculate  directly  from 

according  to 

Pij,(z)  =?}‘^(z)  •  pQut^z)  •  (3.14) 

An  examination  of  S(z)  reveals,  however,  that  the  inverted  function, 
if  evaluated  exactly,  will  be  unbounded  at  t  =  0.  The  use  of  Eq.  (5.l4) 
would  then  lead  to  anomalous  results. 

To  understand  the  behavior  of  ^”^(z)  about  t  =  0,  consider  the 
physical  interpretation  of  the  series  representation  of  d(z)  in  Eq. 
(5.12b).  The  response  lag  in  the  transducer-galvanometer  system  is  due 
mainly  to  the  galvanometer,  which  is  a  second-order  system.  The  output 
from  the  galvanometer  corresponding  to  a  step-input  thus  approximates 
that  calculated  for  a  second-order  system;  namely,  at  t  =  0  the  out¬ 
put  curve  has  zero  deflecticu  as  well  as  zero  slope  (see  Fig.  l4).  For 
zero  deflection,  the  coefficient  of  z^  in  “wst  vanish.  It 

follows  from  Eq.  (5.15)  that  ideally  the  coefficient  r^  in  Eq.  (5.12b) 
would  be  zero  as  well.  In  this  ideal  case,  it  is  easy  to  see  that  the 
inverse  transfer  function  S  ^(z)  would  indeed  be  unbounded  at  t  =  0. 

In  actuality  the  coefficient  r^  is  small  but  not  zero,  owing  to  the 
slight  distortion  of  information  inherent  in  the  data-sampling  and 


The  function  2r(z)  in  the  present  study  involves  the  first-order 
interpolator  Q(z,s)  of  Eq.  (5-10).  In  general,  however,  ?Kz) 
may  involve  a  data  hold  of  any  order. 


reconstruction  process.  The  behavior  of  the  inverse  function  is  never¬ 
theless  unstable  about  t  =  0,  and  this  instability  in  turn  infects  the 
entire  series. 

An  inversion  of  ^^(2)  war  in  fact  atten^Jted  by  first  taking  the 
reciprocal  of  the  right-hand  side  of  Eq.  (3.12a)  and  then  dividing  to 
obtain  a  series  in  z  The  series  thus  obtained  is  divergent,  however, 
as  contrasted  with  the  convergent  series  for  ^^(z)  of  Eq.  (3.12b).  The 
divergence  of  the  series  for  S”^(2)  is  an  indication  of  the  instability 
of  the  inverse  function. 

With  the  failure  of  the  direct  calculation  of  Pj^jj(2)  according  to 

Eq.  (3.13),  an  iterative  scheme  was  devised  by  which  Pj^jj(2)  can  be 

calculated  indirectly  from  Eq.  (3.13)  through  a  series  of  successive 

approximations.  Briefly,  an  initial  guess  is  made  for  P.  (2),  say 
( 0 }  f  0  ^ 

P)  ^(z),  and  the  output  P'  /(z)  corresponding  to  this  guessed  input 

out  f  0  ^ 

is  found  from  Eq.  (3.13).  By  comparing  ^^^(2)  with  the  actual  output 
P„„4.(z)^  it  is  possible  to  derive  a  correction  term  that  provides  an 
in5)rovement  on  ^^i®  forms  the  basis  for  an  iterative  procedure 

that  can  be  repeated  \mtil  the  desired  accuracy  is  obtained. 

The  derivation  of  the  correction  scheme  is  lengthy;  the  details  are 
given  in  Appendix  C.  In  brief,  by  considering  the  magnitude  of  the 
coefficients  r^  of  S(z)  together  with  the  functional  relationship  of 
Eq.  (3.13),  it  can  be  shown  that  an  improvement  on  the  ith  approximation 
to  Pj^jj(z)  at  t  =  nT,  for  n  >  3,  is  given  by 


p(i+l) 

'^in,n 


'  out, n+3 

6 


p(i) 

out^n+3 


f 


(3.15) 


The  inclusion  of  a  first-order  interpolator  introduces  some  incon¬ 
sistency  in  S(z)  about  t  =  0.  Since  the  interpolator  operates  on 
adjacent  data  points,  at  t  -  C  it  interpolates  the  point  t  =  -T  as 
well.  This  results  in  the  physically  unrealistic  term  of  r_^2‘*'^  in 
6(2)  in  addition  to  ascribing  a  nonzero  value  to  r^. 
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where  '•'h®  output  calculated  from  with  the  use  of 

Eq.  (5. 13)*  For  n  =  0,1,2  the  corresponding  expressions  are  slightly 
different  and  are  given  in  Appendix  C.  Convergence  of  the  correction 
scheme  is  not  particularly  rapid,  since  it  is  based  primarily  on  a  linear 
analysis.  The  iteration  procedure  is  simple,  however,  and  as  many  as  5 
to  15  iterations  can  be  performed  rapidly  on  an  electronic  computer. 

To  test  the  accuracy  of  the  foregoing  correction  scheme,  it  was 
applied  to  a  case  in  which  both  the  input  and  the  output  sequences  are 
constructed  and  are  therefore  known  sequences.  This  verification  procedure 
consists  of  the  following  steps:  (l)  construct  a  sequence  of  input  data 
points  that  simulates  the  type  of  Pit^t-pressure  data  that  one  might  expect 
to  obtain  in  the  experiments;  (2)  calculate  the  corresponding  output 
according  to  Eq.  (5.13);  and  (5)  apply  the  correction  scheme  to  this 
output  in  an  attempt  to  recover  the  constructed  input. 

The  constructed  sequences  of  input  and  output  data  from  steps  (l) 
and  (2)  are  shown  in  Fig.  15.  Even  though  the  data  as  used  in  une  actual 
calculation  are  discrete  points  at  the  sampling  instants,  they  are  never¬ 
theless  shown  in  this  figure  as  piecewise  continuous  curves.  This  is 
done  to  irprove  visual  representation.  The  input  cur/e  is  made  up  of 
three  segments,  corresponding  to  an  idealized  representation  of  the  Pitot- 
pressure  variation  in  the  three  distinct  flow  regions  traversed  by  the 
probe  in  going  from  the  tunnel  wall  to  the  wake  axis.  The  first  region 

is  the  tunnel  boundary  layer,  in  which  the  Pitot  pressure  increases 
* 

linearly.  Upon  crossing  the  bow  shock  wave  into  the  inviscid  portion  of 
the  wake,  the  Pitot  pressure  undergoes  an  initial  jvmp,  after  which  the 
pressure  decreases,  rapidly  at  first  but  flattening  out  gradually  as  the 
viscous  wake  is  approached;  a  parabolic  segment  is  used  to  represent 
this  variation.  In  the  third  region,  the  viscous  wake,  the  Pitot  pressure 

In  the  case  where  the  bow  shock  wave  intersects  the  tunnel  boundary 
layer  downstream  of  the  probing  station,  there  is  a  region  between  the 
boundary  layer  and  the  shock  wave  in  which  the  Pitot  pressure  is  con¬ 
stant.  This  region  need  not  be  represented  sepsrntely,  however,  since 
a  line  denoting  constant  pressure  is  essentially  the  same  as  the  line 
representing  the  linear  pressure  increase  in  the  boundary  layer. 


is  taken  to  be  given  by  a  second  parabolic  segment  whose  minimum  corres¬ 
ponds  to  the  axis  of  the  wake.  The  calculated  output  corresponding  to 
this  input  exhibits  the  essential  featiores  of  the  measured  Pitot-pressure 
profiles  (see  Chapter  5  for  coii5)arison). 

In  the  application  of  the  correction  scheme,  the  initial,  guessed 
input  (shown  as  a  dashed  curve)  was  taken  to  be  essentially  the  output 
curve  shifted  backwards  in  time  by  an  appropriate  amount.  A  modification 
was  provided  at  the  jump  discontinuity  by  using  a  backwards  linear  extra¬ 
polation.  The  corresponding  input  data  as  finally  reconstructed  are 
shown  in  the  figure  as  discrete  points.  (To  avoid  crowding  the  figure, 
we  have  plotted  the  input  at  intervals  of  4t  only. )  From  a  comparison 
of  the  reconstructed  input  with  the  original  input,  we  observe  that  there 
is  good  agreement  between  the  two  sets  of  data. 

A  discussion  of  the  application  of  the  correction  scheme  to  the 
measured  Pitot -pres sure  profiles  is  deferred  to  Chapter  b,  when  the  overall 
data-reduction  procedure  is  described. 


4.  SOl;.iCE-FIiOW  EFFECTS  Alffi  CORRECTION 

If.l  Consequence  of  Free-Stream  Nonuniformities 

Tne  Stanford  cpark-heated  tonne]  has  a  conical  nozzle  tha1  includes 
a  portion  nominally  used  and  designated  as  the  test  section  (see  Fig.  2). 
For  the  present  experiments,  however,  the  actual  test  region  includes 
this  test  section  together  with  a  considerable  portion  of  the  nozzle  up¬ 
stream  (see  Section  2.1  and  Fig.  3).  In  what  follows  it  is  therefore 
appropriate  to  use  the  term  "noz?.le"  to  denote  the  entire  test  region. 

The  lindlsturbed  stream  in  the  nozzle  is  nonuniform  in  two  respects. 

The  flow  through  the  nozzle,  starting  from  the  tunnel  throat,  approxi¬ 
mates  tliat  from  a  three -dimens loral  point  source.  Thus  none  of  the  stream¬ 
lines  are  parallel  to  one  another,  and  the  flow  has  transverse  non- 
uniformity.  The  second  aspect  of  the  nonuniforniity  is  actually  a  conse¬ 
quence  of  the  flrr.t,  though  it  is  useful,  to  consider  it  separately.  As 

the  gas  flows  through  the  diverging  nozzle  it  expands,  resulting  in 

* 

longitudinal  gradients  in  the  various  flow  quantities. 

When  a  model  is  placed  in  such  a  nonuniform  stream,  the  flow  field 
around  it  differs  from  that  for  t.ie  same  body  immersed  in  a  uniform 
stream.  Mxthematlcally  speaking,  the  difference  (for  supersonic  flow) 
stems  from  the  different  upstream  boundary  conditions  for  the  two  cases. 
Since  for  most  flow  fields  tiiere  is  a  shock  wave  emanating  from  the  nose 
of  the  model,  this  difference  in  the  upstream  boundary  condition  can  be 
acc 'unted  for  at  the  upstream  side  of  the  shock  wave.  The  flow  field 
around  the  model  accomo-iates  this  different  set  of  upstream  bo-mdary 
conditions  by  making  suitable  adjustments  in  the  shape  of  the  shock  wave, 
with  -.he  result  that  the  shape  of  the  shock  ive  in  source  flow  differs 
from  that  in  uniform  flow  for  the  same  model.  It  is  useful  to  bear  in 


An  exception  to  this  statement  is  found  in  the  velocity.  In  a  hyper¬ 
velocity  tunnel,  the  sp«ed  of  flow  approaches  the  highest  attainable 

value  of  V,  =  (2h.  early  in  the  expansion  process  and  does  not 
1  t^ 

increase  appreciably  thereafter.  Thus  i'  is  permissible  to  assume,  for 
the  sake  of  ease  of  coiJ5)utation  and  analysis,  that  the  imdirturbsd 
stream  velocity  is  constant  in  the  test  region  under  consideration. 
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mind  that  once  this  change  in  shock-vave  shape  has  been  taken  into  consi¬ 
deration,  it  will  bo  incorrect  to  consider  in  addition  that  the  down¬ 
stream  flow  field  is  also  expanding  locally  within  a  diverging  nozzle. 

The  reason  is  that  the  boundary  pertinent  to  the  downsti'eam  flow  field  is 
the  enveloping  bow  shock  wave  and  not  the  nozzle  wall  (or  strictly 
speaking  the  edge  of  the  boundary  layer  that  exists  along  the  nozzle  wall), 
and  the  effects  of  the  source  flow  associated  with  the  diverging  nozzle 
have  already  been  accounted  for  in  the  shape  adjustment  of  the  shock 
wave.  The  sit^iation  changes,  however,  after  the  shock  wave  intersects 
the  nozzle  boundary  layer.  In  this  situation,  th '  flow  in  the  region 
bounded  by  the  reflected  shock  wave  and  the  edge  of  the  boundary  layer  is 
dii’ectly  affected  by  th'  diverging  geometry  of  the  edge  of  the  boundary 
layer. 

In  the  present  experiments,  the  wake  Pitot  pressure  is  among  the 
flow  q.uantities  affected  by  the  tunnel  source  flow.  To  use  the  pressure 
data  for  comparison  w^ch  any  theoretical  results  worked  out  for  the  case 
of  uniform  flow,  the  data  first  must  be  corrected  for  so’jrce-flow  effects. 
Such  correction  would  also  be  inportant  if  the  data  were  to  be  applied 
to  examples  of  actual  flight,  for  which  a  uniform  atmosphere  pre^^ails. 

Ve  shall  consider  the  effects  or.  .y  insofar  as  they  change  the  shape  of 
the  shock  wave  and  not  after  interaction  between  shock  wave  and  boundary 
layer  has  occurred.  In  the  following  sections,  we  discuss  the  einalytical 
basis  for  approximate  schemes  used  to  correct  the  measured  Pitot  pressure 
for  source-flow  effects. 

4.2  Tunnel  Source  Flow 

Before  proceeding  to  the  main  portion  of  the  analysis,  it  is  appro¬ 
priate  to  consd.der  the  nature  and  magnitude  of  the  source  flow  in  the 
Stanford  tunnel. 

The  flow  in  the  tunnel  icay  be  regarded  as  consisting  of  an  inviscia 
core  flow  and  a  thick  hypersonic  boundary  layer  on  the  tunnel  wall.  If 
the  thickness  of  the  boundary  layer  along  the  nozzle  is  a  constant  fraction 
of  the  local  nozzle  radius,  then  the  edge  of  the  boundary  layer  will  form 
a  cone  withi-i  the  nczzle  cone,  the  apexes  of  both  cones  being  at  the 
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throat  of  the  nozzle.  In  this  case,  the  inviscid  flow  will  be  that  due 
to  a  point  source  situated  at  the  throat.  Theoretical  results  supported 
by  experiment  show,  however,  that  the  edge  of  the  boundary  layer  is  not 
a  conical  surface.  The  boimdeoy  layer  thus  has  a  contouring  effect  on 
the  inviscid  flow,  auid  the  magnitude  of  the  local  source-flow  effect  is 
different  from  that  associated  with  the  ia  ..  case  of  a  point  source 
positioned  at  the  throat.  We  shall  present  a  comparison  of  theory  with 
experiment,  erfter  which  the  magnitude  of  the  locelL  source-flow  effects 
will  be  assessed. 

The  problem  of  boundary-layer  development  in  a  hypersonic  nozzle  has 
been  studied  analytically  by  Rasmussen  and  Karamcheti  (1966).  Results 
from  this  work  (Fig.  6)  can  be  used  to  deduce  the  variation  of  the  area 
ratir  Aj,/A  with  the  longitudinal  distance  L  along  the  nozzle.  Here 
the  effective  area  of  the  inviscid  flow,  A  is  the  tunnel  cross- 
sectional  area,  and  L  is  meac’ired  from  the  throat.  When  the  variation 
of  the  area  ratio  A^/A  is  known,  it  is  possible  -co  calculate  the 
variation  of  other  flow  quantities.  For  example,  to  find  biie  ratio  of 
the  Pitot  pressure  at  any  two  points  along  the  length  of  the  nozzle,  one 
starts  by  writing  the  Pitot  pressure  p^  us 

p^  =  f(p,T,V)pV^  .  ;4.1) 

In  this  equation  the  function  f(p,T,V)  is  essentially  a  constant  for  a 
vide  range  of  conditions  when  the  Mach  number  is  high.  Furthermore,  the 
velocity  may  be  assiuaed  to  be  constant  along  the  length  of  the  nozzle 
iinder  consideration  (see  footnote  on  p.  34).  Making  use  of  mass  conser¬ 
ve  tion,  we  then  have,  for  two  pojjits  a  and  b. 


This  is  an  idealization  based  on  the  fact  that  the  throat  diameter  is 
very  small  relative  to  the  other  dimensions  of  the  conical  noz7.1e. 
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a  a 

(4.2) 

vhere  A^/A^  is  calciilated  directly  from  the  nozzle  geometry. 

It  is  now  possible  to  compare  the  theoretical  results  with  experi¬ 
mental  date  and  this  is  done  in  Fig.  l6.  Two  theoretical  curves  are 
shown.  The  solid  curve  is  a  plot  of  Eq.  (4.2)  in  which  the  ratio  A^/A, 
deduced  from  the  theory  of  Rasmussen  and  Karamcheti,  varies  from  0.4l6 
to  0.475«  The  second  theoretical  curve  (dashed  line)  is  based  or  the 
assumption  that  the  Pitot  pressure  varies  inversely  as  the  tunnel  irea 
A  instead  of  the  effective  area  A„.  This  is  equivalent  to  saying  that 

£t 

in  Eq.  (4.2)  the  ratio  A^/A  is  tsiken  to  be  constant,  that  is,  that  the 
edge  of  the  boixndary  layer  forms  a  conical  sui'face  with  apex  at  the  nozzle 
throat.  The  difference  between  the  two  sec.  of  results  indicates  that  the 
inviscid  flow  in  the  tunnel  as  predicted  by  the  theory  of  Rasmussen  euid 
Karamcheti  is  different  from  that  in  the  ideal  case  of  a  source  sitiated 
at  the  nozzle  throat.  Also  shown  in  Fig.  l6  are  experimental  points 
deduced  from  data  obtained  in  a  jing.le  run  with  the  axially  traversing 
probe.  The  good  agreement  between  the  experimental  results  and  the 
theory  of  Rasmussen  and  Ksramcheti  verifies  the  earlier  statement  tlia'*  the 
boundary  layer  has  a  contouring  effect  on  the  inviscid  flow.  Admittedly 
the  range  in  which  the  theory  is  substant iated  by  the  data  is  sanewhat 
limited,  owing  to  the  relatively  short  traverse  (9  inches)  of  the  probe. 

In  view  of  the  close  agreement,  however,  we  shall  adopt  the  results  of 
Rasmussen  and  Karamcheti  for  present  purposes. 

The  Pitot-pressure  correction  scheme  worked  out  in  later  sections 
(Sections  4.5  and  4.6)  is  based  on  cheory  tiiat  pertains  to  a  siinple  source 
flow.  To  render  the  theory  applicable  to  the  present  situation  where  the 
stream  is  not  a  simple  source  flow,  the  following  approximation  is  made 
in  regard  to  the  stream  nonuniformity.  It  is  assumed  that  the  local 
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nonuiilform  stream  can  be  replaced  by  a  virtual  source  that  gives  at  the 
station  of  interest  the  same  longitudinal  flow  gradients  as  those  calcu¬ 
lated  from  the  theory  of  Rasmussen  and  Karamchetl  (see  Fig.  lT(a)  for 
illustration).  A  characteristic  length  associated  with  this  virtual 
source  is  the  distance  from  the  origin  of  the  source  to  the  nose 

of  the  model,  that  is,  to  the  station  of  Interest.  This  length  is  used 
to  define  a  "source-flow  parameter"  where  r^  is  a  cheiracter- 

istlc  dimension  of  the  model,  here  taken  to  be  the  radius  of  the  sphere. 

The  source-flow  parameter  gives  a  measvire  of  the  seriousness  of  the  source- 
flow  effects  on  the  flow  field  around  the  model. 

The  problem  is  now  to  find  a  relation  that  expresses  in  terms 
of  the  longitudinal  gradients.  With  this  relation,  one  can  then  use  the 

results  of  Rasmussen  and  Karamchetl  to  find  the  variation  of  L  with  the 

o 

axial  distance  L  from  the  nozzle  throat.  To  derive  the  required  relation, 
consider  the  sketch  shown  in  Fig.  lT(b).  The  two  points  a  and  b  with¬ 
in  the  vlrtusl  source  flow  are  situated  at  distances  L  and  (L  +/CSL  ) 

o  ^  o  o 

from  the  origin  of  the  vlrtioal  source.  From  Eq..  (4.2),  we  have 


which  leads  to 


In  the  limit  as 
relation 


(p')i.  A  /L  +Z5L  \-'d 
'^o'b  a  _  (  o  o\ 


L_+i£L_  \-2 


(4.5) 


1  + 


TF 


(4.4) 


b  approaches  a,  Eq.  (4.4)  reduces  to  the  required 


Since  Eq.  (4.2)  deeCLt  with  the  tunnel  stream,  the  appropriate  area  to 
use  was  the  effective  area.  In  the  case  of  the  virtual  source  flow, 
however,  it  is  not  necessary  to  use  the  "effective"  subscript. 
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The  quantity  (dp'/dL)  is  evaluated  gi-aphically  from  the  Pitot -pres  sure 
curve  deduced  from  the  results  of  Rasmussen  and  Karamcheti  and  shown  in 
Fig.  l6.  The  variation  of  thus  calculated  is  given  in  Fig.  l8  as  a 
function  of  L.  Shown  for  comparison  is  the  straight-line  relationship 
=  L  tliat  would  prevail  if  the  edge  of  the  boundary  layer  were  a  conical 
surface  with  apex  at  the  throat.  It  can  be  seen  that  the  distance 
varies  from  20  inches  at  L  =  80  inches  to  32  inches  at  L  =  110  inches. 
Thus  the  origin  of  the  virtual  source  is  at  a  dl  tance  of  20  to  32  i..ches 
from  the  throat.  Considered  in  terms  of  the  overall  values  of 
involved  (100  to  140  inches),  this  gradual  shift  of  12  inches  is  a  rela¬ 
tively  minor  change.  Thus  the  approach  used  here  of  representing  the 
local  nonunifom  stream  by  an  idealized  point  source  appears  to  be  rea¬ 
sonable. 

The  description  of  the  source  flow  is  now  complete.  In  the  next 
section  a  discussion  is  given  of  the  available  solutions  to  problems  of 
flow  around  objects  placed  in  a  source  flow.  We  then  turn  to  the  specific 
problem  of  source-flow  effects  on  waives  and  begin  with  a  general  descrip¬ 
tion  of  the  wake  behind  a  sphere.  Particular  empliasis  is  placed  on  showing 
that  some  of  the  available  source-flow  solutions  cein  be  applied  to  the 
wake  problem.  The  theory  of  the  correction  scheme  for  wake  flow  is  then 
treated  in  Sections  4.5  and  4.6. 

4.3  Available  Source-Flow  Solutions 

Source-flow  effects  were  recognized  to  be  significant  only  after  a 
considerable  number  of  hypervelocity  tvinnels  had  been  put  into  operation. 

It  WC.8  then  found  that  substantial  errors  could  result  from  the  conduct 
of  tests  in  a  diverging  stream.  This  realization  stimulated  a  consid¬ 
erable  amount  of  analyticeil  work  aimed  at  finding  ways  for  correcting  the 
experimental  data. 

The  first  published  work  on  source-flow  effects  is  due  to  Baradell 
and  Bertram  (i960).  They  use  the  method  of  characteristics  to  calculate 
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the  flow  field  about  a  blunted  flat  plate  in  the  flow  from  a  point  source. 
The  solution  for  the  nose  region  is  eui  approximation,  since  the  body- 
surface  here  is  replaced  by  a  so-called  "sonic  wedge".  Several  other 
studies  make  use  of  simple  Newtonian  theory  to  estimate  the  effects  on  the 
surface  pressure  on  slender  bodies.  This  is  done  by  Burke  and  Bird  (1962) 
and  Whitfield  and  Norfleet  (I962)  for  wedges  and  cones,  and  by  Hall  (1965), 
who  includes  the  Busemann  centrifiogal  term  for  flow  over  curved  surfaces. 

The  work  of  Hall,  in  particular,  extends  the  application  of  the  simple 
Newtonian-flow  model  to  generalized  free-stream  nonuniformities.  Hall 
also  treats  the  problem  of  a  blxinted  slender  body  (typically  a  hemisphere- 
cylinder  combination)  by  performing  a  perturbation  analysis  on  a  flow 
model  based  on  the  detached-shock-wave  theory  of  Cheng  (I96O).  The  source- 
flow  perturbations  of  the  shock-wave  location  and  surface  pressure  are 
thereby  estimated. 

A  different  approach  is  used  by  Bemdt  (I962),  who  adopts  the  constant- 
density  approximation  to  estimate  the  effect  on  the  shock-wave  stand-off 
distance  and  the  stagnation -point  velocity  gradient  for  a  blunt  nose. 
Following  a  suggestion  by  Bemdt,  Meyer  (1965)  used  the  blast-wave  analogy 
to  treat  the  problem  of  a  hemisphere-cylinder  model  in  the  flow  from  a 
point  source. 

More  recently,  Sa-vage  (I966)  has  used  the  Newtonian  successive- 
approximation  scheme  of  Cole  (195?)  to  calculate  the  effect  on  the  flow 
field  about  a  slender  cone.  This  method  of  solution  represents  a  more 
systematic  approach  to  the  problem  than  does  the  straightforward  pertur¬ 
bation  method  based  on  the  slii5)le  Newtonian-flow  model.  The  results  of 
Savage,  worked  out  to  second  order  in  the  Newtonian  expansion,  are  not 
significantly  different,  however,  from  those  of  Whitfield  and  Norfleet, 
and  of  Hall. 

As  regards  numerical  solutions,  in  addition  to  the  early  work  of 
Baradell  and  Bertram,  there  are  two  recent  studies  by  Eaves  and  Lewis 
(1965)  and  Inouye  (1966a).  Saves  and  Lewis  consider  a  hemisphere -cylinder 
model  in  which  source-flow  effects  are  included  in  the  method-of- 
characteris-Mcs  calculations  but  neglected  in  the  solution  for  the  sub¬ 
sonic-transonic  flow  in  the  nose  region.  That  is,  the  blunt-body  solution 
is  for  a  uniform  parallel  flow,  whereas  in  the  supersonic  portion  of  the 
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flow  field  the  upstream  boundary  condition  at  the  shock  wave  corresponds 
to  the  case  of  nonuniformities  resulting  from  a  point  source.  Also 
included  in  this  analy..ls  ore  viscous  effects  over  the  model  and  their 
interaction  with  the  inviscid  floii?  field.  The  work  of  Inoviye,  on  the 
other  hand,  treats  exclusively  the  inviscid  problem.  The  solution  is 
exact  in  that  source-flow  effects  are  accoimted  for  in  both  the  subsonic- 
transonic  region  and  the  supersonic  region.  The  solution  for  these 
regions  is  based  respectively  on  the  blunt-body  solution  of  Van  Dyke 
(1958)  and  on  the  method  of  characteristics.  The  model  configurations 
considered  Include  a  spherical  nose,  ellipsoidal  nose,  spherically  blunted 
cone,  and,  of  particular  interest  to  the  present  investigation,  a  hemi¬ 
spherical  nose  with  a  cylindrical  afterbody. 

To  see  which  of  the  aforementioned  solutions  are  applicable  to  the 
current  problem  of  source-flow  correction  we  first  review  briefly  the 
general  character  of  the  wake  behind  a  sphere.  This  is  done  in  the  next 
section,  which  also  gives  a  preliminary  discussion  of  how  the  pertinent 
solutions  '.Jill  be  used. 

4.4  The  Wake  Behind  a  Sphore 

For  the  present  study,  we  are  not  cuncemed  directly  with  the  near¬ 
wake  region  upstream  of  the  wake  neck.  The  nea:'  wake  will  be  considered 
only  in  terras  of  its  gross  effect  on  the  wake  downstream;  the  complicated 
structure  of  the  recirculation  zone  will  be  ignored  for  the  most  part. 

The  remainder  of  the  wake  (downstream  of  the  neck)  is  assumed  to  consist 
of  two  distinct  regions,  the  outer  inviscid  wake  and  the  inner  viscous 
wake  (see  Fig.  l). 

The  flow  in  the  outer  inviscid  region  consists  primarily  of  the 
familiar  entropy  layer  and  shock  layer.  The  overall  picture  of  the 
inviscid  wake  behind  the  sphere  is  comparable  to  the  aft  portion  of  the 
flow  field  around  a  hemisphere-cylinder  model.  This  similarity  is  based 
on  the  following  considerations: 

(1)  The  inviscid  flow  field  around  a  blunt-nosed  body  is  dominated  by 
the  bow  shock  wave,  which  in  turn  is  dictated  by  the  shape  of  the  nose. 

The  sphere  and  the  hemisphere-cylinder  have  identiceil  nose  shapes. 
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(2)  Even  if  the  two  flow  fields  are  not  the  same  in  every  detail,  the 
overall  picture  of  the  inviscid  flow  is  affected  only  slightly  by  the 

differences.  To  illustrate  this, 
we  construct  an  Inviscid  flOT>? 
field  for  the  sphere  in  which 
the  supersonic  stream  comes  to¬ 
gether  at  the  rear  (see  accom- 
panying  sketch  ).  By  construction 
the  space  that  is  taken  up  by  the 
cylindrical  afterbody  in  the  case 
of  the  hemisphere-cylinder  is 
now  occupied  by  fluidj  the  two 
flow  fields  are  therefore  some¬ 
what  different.  This  difference  does  not  significantly  affect  the  rest 
of  the  inviscid  flow  field,  however.  This  is  because  the  fluid  near  the 
longitudinal  axis,  being  ijart  of  the  entropy  layer,  is  of  low  density. 

In  addition,  the  volume  difference  that  is  under  consideration  is  rela¬ 
tively  small.  Thus  the  details  of  the  afterbody  geometry  do  not  play  a 
crucial  role  in  determining  the  overall  inviscid  flow  field.  In  the 
actual  flow  behind  the  sphere  tne  existence  of  the  viscous  wake  modifies 
the  above  picture  in  any  case.  The  viscoiis  wake,  may  be  thought  of  in  a 
sense  as  replacing  the  cylindrical  afterbody.  The  amount  of  fluid  that 
flows  through  the  viscous  wake  is  small,  however  and  the  inviscid  wakes 
may  still  be  considered  similar  to  each  other. 

(3)  Insofar  as  the  inviscid  wake  is  concerned,  the  recorapression  shock 
wave  emanating  from  the  neck  region  has  been  shown  by  Feldmain  (I961)  to 
have  negligible  effect  and  will  be  disregarded  in  the  analysis  to  follow. 
Here  again  the  equivalence  of  the  two  wakes  is  substantiated. 

In  view  of  the  similarity  of  the  inviscid  wakes  for  the  sphere  and 


The  conical  region  behind  the  sphere,  drawn  in  dashed  line,  may  be 
interpreted  as  being  either  a  dead-water  region  or  a  solid  afterbody. 
This  extension  in  boundary  is  necessary  because  otherwise  the  stream 
would  have  to  be  ttimed  through  90“  at  the  rearmost  point  on  the 
sphere,  which  is  not  possible  for  a  supersonic  stream. 
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the  hemisphere -cylinder,  all  the  source-flow  solutions  worked  out  for 
the  latter  case  can  be  applied  to  the  inviscid  wake  of  the  sphere.  In 
particular  the  solutions  of  Meyer  (I965),  Hall  (1965),  and  Inouye  ( 1966a, b) 
will  be  used  simultaneously,  as  well  as  first-order  blast-wave  theory. 

Wliere  possible  the  approximate  solutions  of  rte.ll  and  the  results  from 
blast-wave  theory  are  preferred  because  of  their  inherent  simplicity. 

The  work  of  Inouye  is  invaluable,  however,  in  that  it  provides  a  means 
for  checking  the  accuracy  of  the  other  methods. 

Turning  to  the  inner,  viscous  portion  of  the  wake,  we  note  that  it 
has  its  origin  in  the  boundary  layer  on  the  sphere.  The  boundary  layer 
separates  near  the  meridian  cp  =  3t/2  (<p  defined  in  Fig.  l),  despite  the 

existence  there  of  a  favorable  pressure  gradient  in  the  external  inviscid 
flow.  Separation  in  this  instance  is  caused  by  high-pressure  signals 
that  originate  in  the  doimstream  coii5)ression  region  near  the  wake  neck  and 
propagate  upstream  throiogh  the  subsonic  portion  of  the  wake  and  boundary 
layer.  The  boundary  layer  after  separation  becomes  a  free  shear  layer 
in  which  mixing  takes  place  with  the  external  inviscid  flow.  Viscous- 
inviscid  intei-action  then  characterizes  the  compression  process  that 
accompanies  the  tiirning  of  the  shear  layer  back  to  parallel  with  the  wake 
axis.  The  resultant  inner  wake  is  relatively  hot,  and  viscous  effects 
are  then  important.  A  short  distance  after  the  annular  shear  layer  has 
coalesced  into  the  viscous  wake  proper,  the  flow  becomes  supersonic 
tliroughout.  Since  the  flow  in  the  viscous  wake  is  of  the  bovindary- layer 
type,  the  static  j>ressure  at  a  given  longitudinal  location  is  viniform 
and  equal  to  the  vai^ue  in  the  external  inviscid  flow  at  that  location. 

One  may  thus  think  of  an  invisc id-pressure  "envelope"  svurrounding  the 
viscous  wake.  In  general  the  development  of  the  viscous  wake  depends  both 
on  this  inqjosed  pressure  field  and  on  the  diffusion  processes  present. 

For  the  conditions  of  the  current  tests,  however,  the  pressure  field  has 
the  more  important  effect,  as  will  be  shown  in  Section  4.6.  To  the  extent 
that  this  is  true,  the  toviscid  solutions  for  the  hemisphere-cylinder 
model  are  again  applicable  in  the  analysis  of  the  viscous  wake.  In 
addition,  since  the  viscous  wake  originates  in  the  boundary  layer,  a 
consideration  of  so\irce-flow  effects  on  the  viscous  wake  must  necessarily 
include  the  related  effects  on  the  boundary  layer  as  well.  To  study 
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these  latter  effects,  we  need  the  Inviscid  pressure  distribution  over 
the  nose  region,  which  the  numerical  work  of  Inouye  (1966a)  alone  pro¬ 
vides,  since  the  studies  of  Meyer  and  Ilall  do  not  treat  the  nose  region 
in  detail  but  only  account  for  the  cianulative  effect  of  the  nose  shape 
through  the  overall  nose  drag. 

We  now  proceed  to  a  more  detailed  consideration  of  the  source-flow 
cori-ection  schemes  in  the  next  two  section. 

4.5  Source-Flow  Correction  for  Inviscid  Wake 

It  was  pointed  out  in  the  last  section  that  the  inviscid  wake  behind 
a  sphere  is  essentiaily  equivalent  to  the  aft  portion  of  the  inviscid 
flow  field  of  a  hemisphere-cylinder  model.  The  latter  model  in  a  source 
flow  has  been  studied  by  several  investigators  (see  Section  4.3).  In 
these  studies  the  flow  variable  of  primary  interest  has  been  the  static 
pressure  on  the  surface  of  the  model.  Tbie  results  show  that  a  source  flow 
has  considerable  influence  on  the  magnitude  of  the  surface  pressure,  as 
weil  as  the  static  pressure  in  general..  To  illustrate  this,  the  normal¬ 
ized  pressure  difference,  (pg  -Pg  )/p2  ,  as  calculated  from  unpublished 

results  of  Ino'oye  ,  is  presented  in  Fij.  19  as  i  function  of  the  normal¬ 
ized  radial  distance  f  =  r/R.  Here  pg  denotes  the  static  pressure  in 
the  wake  (that  is,  downstream,  of  the  bow  shock  wave),  the  subscripts 
"u"  and  "s”  refer  to  uniform  flow  and  source  flow  respectively,  and  R 
is  the  distance  from  the  axis  of  syrametiy  to  the  bow  shock  wave.  Tne 
results  are  for  an  x/d  of  5  (x  and  d  defined  in  Fig.  l),  which  is 
the  farthest  station  considered  by  Inouye.  The  values  of  pg  for  source 

flow  are  taken  from  the  case  with  r. /L  =  0.02,  the  smallest  value  of 

b'  o 

the  source -flow  jiarameter  studied  by  Inouye.  Fo„’  the  particxilar  case 

considered,  the  q  lantity  (pg  -pg  )/pg  ranges  from  0.35  to  0.40.  Thus 

u  s  a 

the  difference  in  the  static  pressure  is  an  order  of  magnitude  larger  than 
the  source-flow  parameter  and  should  not  be  ignored. 


The  author  wishes  to  thank  Mr.  M.  Inouye  of  .<biies  Research  Center  for 
providing  these  results. 
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For  the  present  investigation  the  quantity  of  more  direct  interest 

is  the  vaJce  Pitot  pressure  p^.  This  can  be  deduced  from  the  static 

pressure  and  the  Mach  number,  the  latter  being  sQ-so  available  from 

Inouye's  unpublished  results.  The  fractional  difference  (p'  -p'  )/p* 

u  °s  °u 

thus  calculated  is  plotted  along  with  the  static -pressure  results  in 
Fig.  19-  Although  the  difference  here  is  not  as  large  as  for  the  static 
pressure,  it  is  still  from  ten  to  fifteen  times  the  value  of 
the  present  experiments  the  source-flow  parameter  assumes  values 

up  to  0.01.  On  the  basis  of  dependence  on  alone  (if  we  assume 

a  linear  dependence)  errors  of  10  to  15  percent  may  thus  be  anticipated 
in  the  Pitot  pressure  in  the  most  serious  case.  The  quantity  is 

not  the  only  controlling  parameter,  however;  source-flow  effects  also 
Increase  with  the  distance  x/d.  In  the  present  tests  the  station  with 
the  largest  value  of  is  at  an  x/d  of  about  15,  compared  with 

x/d  =  5  for  the  case  considered  in  Fig.  19*  Because  of  the  added  source- 
flow  effect  due  to  the  larger  value  of  x/d,  the  discrepancy  in  Pitot 
pressure  can  thus  be  expected  to  be  substantially  larger  than  10  to  15 
percent. 

Besides  affecting  the  magnitude  of  the  Pitot  pressure,  the  source 
flow  also  changes  the  lateral  extension  of  the  wake.  It  is  thus  desirable 
to  find  a  means  of  correcting  for  both  these  effects.  An  immediately 
obvious  approach  is  the  direct  use  of  method-of -characteristics  results. 
The  calculations  that  this  method  entails  are  involved  and  lengthy,  how¬ 
ever,  and  we  seek  to  develop  an  approximate  scheme  that  makes  use  of 
relatively  simple  analytical  expressions.  In  this  derivation,  Inouye’s 
more  accurate  method-cf -characteristics  results  will  be  used  either  to 
check  the  approximate  theory  eii5)lQyed  or  to  justify  the  assur^jtions  that 
are  made. 

We  shall  first  outline  the  essential  steps  in  the  correction  method, 
after  ;.’hlch  the  var-ious  points  that  require  explanation  will  be  discussed 
in  detail. 

(l)  Up  to  the  point  where  r  is  0.7  to  0.8,  the  variation  of  the 
total  pressure  with  r  is  assumed  to  be  the  same  for  uniform  flow  and 
source  flcr.;. 
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(2)  At  a  given  longitudinal  station,  the  variation  of  the  static > 

pressuT-'!  ratio  PgCrVPg  vith  r  is  assumed  to  be  luichanged  from  uni- 

e 

form  flow  to  source  flow.  Here  Pg  *=  PoC^g)  Is  the  pressure  at  the 

e 

surface  of  the  cylindrical  afterbodj'",  which  in  the  case  of  the  wake  behind 
a  sphere  is  equivalent  to  the  hypothetical  ’'surface "  separating  the  viscous 
emd  inviscid  wakes,  or  the  "edge"  of  the  viscous  wEike.  Since  the  ratio 
P2(r)/pg  caiuiot  be  deduced  from  experimental  data,  we  use  the  results 

from  first-order  blast-wave  theory,  which  compare  favorably  with  the 
method-of-cl'aracteristics  results. 

(3)  The  longitudinal  variation  of  che  surface  pressure  for  uniform  flow, 

(pg  )^,  is  given  satisfactorily  by  first-order  blast-wave  theory .  The 
e 

corresponding  pressure  in  source  flow,  (p^  )  ,  if  obtained  by  combining 

e  ® 

blast-wave  theory  with  a  simple  correction  fonmila  due?  to  Hall  (1963). 

An  alternative  approach  that  applies  to  the  vertically  traversing  probe 

is  to  deduce  (po  )  at  the  particular  station  from  experimenual  results, 
e 

and  then  use  Hall's  correction  formula  to  obtain  (pg  )  .  In  either  case, 

with  the  aid  of  (2)  the  transverse  distributions  Pg  (r)  and  Pg  (r) 

u  ‘^s 

are  ascertained  Independently  of  the  method-of-characteristics  results. 

(4)  With  Pg  (?)  known  from  (3),  we  then  use  the  measvired  Pitot-press\jre 

s 

distribution  p'  (r)  to  deduce  the  distribution  of  the  total  pressure 
s 

(p+  (r)]  •  This  is  done  throvigh  the  use  of  routine  normal-shock  relations, 
tg  s  . 

5y  virtue  of  condition  (l),  the  determination  of  [p+  (r)]„  also  gives 
•the  total -pressure  distribution  in  imiform  flow,  [p^ 

(5)  With  the  distribution  [p.  (r)]  from  (4)  and  Pg  (r)  from  (3), 

the  reqviired  Pitot-pressure  distribution  in  uniform  flow,  p’  (r),  can 

u 

be  aeduced  accordingly. 

(6)  The  jateral  distance  in  the  Inviscid  wake  in  source  flow,  r^,  is 
corrected  by  a  simple  scaling  down  of  the  distance  with  the  ratio 
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vhere  R  is  the  distance  from  the  vake  axis  to  the  shock  vave. 


In  the  foregoing  outli.ie,  steps  (U)  and  (5)  are  straightfonard  euid 
need  no  further  expl-jiation  at  this  point.  .Steps  (l),  (2),  (5)^  and 
(6)  on  the  other  hand  req.uire  elaboration.  Since  the  theory  in  support 
of  these  steps  dravs  from  the  work  of  Meyrr  (1963),  Hall  (1965),  and 
Inouye  (1966a, b),  we  now  describe  these  solutions  in  brief  and  note  some 
of  the  results. 

Meyer  solved  the  inviscid  problem  of  a  hemisphere-cylinder  model  in 

a  source  flow  by  application  of  the  blast-wave  euialogy.  Formulating  the 

problem  in  the  manner  of  the  second-order  blast -wave  theory  of  Sakurai 

(1953,  195^)  and  accounting  for  source-flow  effects  through  the  upstream 

boundary  condition  at  the  bow  shock  wave,  Meyer  showed  that  source-flow 

effects  come  into  play  only  in  the  second-order  terms,  where  they  appear 

through  a  blast-wave  sovirce-flow  parameter  k  defined  by 

s 


In  this  definition  is  a  function  of  7  alone  and  is  C.877  for 
7  =  l.k  (see  Sakirrai,  1955),  Cj^  is  the  drag  coefficient  of  the  nose, 

is  the  Mach  n^mLer  upstream  of  the  shock  wave  at  the  nose,  r^  is 
n 

the  nose  radius  of  the  hemisphere-cylinder,  and  is  the  distance  from 
the  nose  to  the  origin  of  the  source.  With  the  source-flow  parameter 
k  as  the  parametric  variable,  solutions  are  obtained  by  namerical  inte- 
gration  for  the  pressxire,  density,  radial  component  of  velocity,  and 
shock -wave  shape. 

The  second-order  results  from  Meyer's  work  that  are  of  interest  to 
the  present  study  are  those  for  static  pressure  and  shock-wave  shape. 
Although  none  of  these  results  are  actually  employed  in  the  correction 
scheme,  they  are  useful  in  the  conpariscn  of  the  relative  accuracy  of 
various  theories.  It  is  through  this  comparison  that  we  are  able  to  show 


Later  (in  Section  5»5)  when  the  correction  scheme  is  applied  to  ewtual 
test  data,  some  of  the  details  of  these  In^o  steps  will  be  discussed. 


that  the  second-order  results  do  not  represent  a  significant  improvement 
over  fiist-order  blast-vave  theory.  Furthermore,  the  first-order  results, 
vhich  are  adopted  for  use  in  the  correction  scheme,  are  contained  in  the 
second-order  results.  For  these  reasons  ve  sumtnarize  the  pertinent  re¬ 
sults  from  Meyer's  second-order  solution. 

Meyer's  result  for  the  static  pressure  can  be  vritten  as 


P(x,r;kg) 


(0)/*x 

— {1  +  (|r(r;k  )e(x;k  )}  , 

6(x;k^)  ® 


vhere  x  =  x/d,  €  is  the  pen:«meter  in  the  blast-vave  expansion,  and 

and  are  respectively  the  first-  and  second-order  coefficients 

in  the  expansion.  The  quantity  e  is  given  by 


vith 


8  _  X 


n  n 


and 


vM  ■ 

n 


(4.8) 


(4.9) 


(4.10) 


Source-flow  effects  enter  into  Bqs.  (4.8)  and  (4.9)  through  the  parameter 
X,  =  -1.992  -  2.695  k„  ,  for  7  =  1.4  .  (4.11) 

X  B 

/Q>  ^ 

'  ^  in  Eq.  (4.7)  depends  on  r  and  7  but  not  on  k„ 

0 


The  function  g 


and  Is  given  by 


y+1  _  1  -1 

r  "7 


-  1(72+27-1) 


(4.12) 


The  quantity  ’Krjk  )  is  not  expressible  in  closed  form  and  must  be 

obtained  by  numerical  integration.  Meyer's  results  for  \|f  are  shown  in 

Fig.  20  fui  values  of  k  ranging  from  0  to  1.6. 

s 

It  is  now  possible  to  derive  a  simplified  expression  for  the  surface 
pressure.  Since  there  is  practically  no  radial  gradient  in  the  pressure 
near  the  body  surface  —  this  is  confirmed  by  blast-wave  theory  —  it  is 
customary  in  balst-vave  applications  to  use  the  pressure  at  the  axis 
(r  =  0)  as  the  sxirface  pressure.  The  variaoion  of  with  k^  for 

this  special  case  is  viven  by  the  linear  relation 


=  -0*662  -  5.596  k^ 


(4.13) 


When  Eqs.  (4.8)  through  (4.13/  are  used  in  Eq.  (4.?)/  the  resultant 

pressure  can  be  expanded  in  the  parameter  2X^(j o/70p)^^S^  /x.  To  second 

u 

order  the  fellow log  is  obtained  for  the  surface  pressure: 


—  =  0.0671  —2  +  0.355 

Pi  0  - 


-  0.669  k  , 

s' 


for  7  =  l.k  .  (4.1k) 


The  value  of  for  the  hemispherical  nose  is  taken  to  be  1.0. 

In  first-order  blast-wave  theory  (see  Lukasiewicz,  I96I),  the  surface- 
pressure  result  is  sinqjly  the  first  term  of  the  right-hand  sidf*  of  Eq. 
(4.14): 


— S  =  0.0671  Vc,, . 

Pi  ^  X 
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(4.15) 


The  first-order  result  for  the  twunsverse  vailatlon  of  the  static 
prisssure,  vill  also  be  needed.  This  Is  most  conveniently  given 

by  the  ratio  P2(r)/P2  *  which  we  have  from  Eq.  with  the 

second-order  term  omltled; 


\  g^°^(0j7) 


ik.lC) 


The  right-hand  side  of  this  equation  Is  easily  calculated  from  Eq.  (4.12). 

As  to  the  t'tudy  of  Hall  (1963),  source-flow  solutions  for  the  hemi¬ 
sphere-cylinder  model  are  obtained  explicitly  for  surface  pressure  and 
shock-wave  location.  In  this  work  the  quantity  Is  used  as  the 

source-flow  parameter.  We  will  have  occasion  to  use  the  ratios 

(Pn  )  /(P/>  }  Sind  H  /R  ,  for  which  Hall's  results  are 
e  e 


(4.17) 


(4.18) 


We  are  now  In  a  position  to  present  the  supporting  evidence  and 
provide  the  specific  formulas  for  steps  (l),  (2),  (3),  and  (6)  In  the 
correction  procedure  outlined  earlier. 

(1)  In  the  Invlscld  flow  behind  the  bow  shock  wave,  the  total  pressure 

p.  along  any  given  streamline  Is  a  constant.  Its  value  changes  from 

^2 

streamline  to  streamline,  however,  because  the  fluid  along  the  various 

streamlines  has  traversed  portions  of  the  bow  shock  wave  that  have 

different  angles  of  inclination.  Nevertheless,  the  value  of  p.  Is  a 

^2 

function  only  of  the  angle  of  Inclination  of  the  portion  of  the  shock 
nave  traversed  by  the  fluid  In  question.  This  Is  because  the  toted, 
pressure  ahead  of  the  shock  wave,  p.  ,  Is  the  same  everywhere  In  the 

h. 
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free  stream  in  both  the  \inlform  flow  and  the  source  flow.  Now,  at  a 
given  longitudinal  station  where  the  shock-wave  inclination  is  not 
■•Arge,  the  streamlines  that  are  not  too  close  to  the  shock  wave  are  al¬ 
most  horizontal.  The  fluid  here  has  crossed  the  shock  wave  at  points 
considerably  farther  upstream.  In  particular,  if  consideration  is  liioited 
to  the  entropy  layer  next  to  the  body,  the  fluid  there  has  traversed  the 
segment  of  the  bow  shock  wave  that  is  nearly  normal  to  the  flow  axis. 

The  slope  at  corresponding  points  of  this  segment  does  net  dJi'fer  signif¬ 
icantly  between  uniform  flow  and  soiree  flow.  This  can  be  confirmed  by 
comparing  values  of  (dR/dx)  and  (dR/dx)  calculated  from  Eq.  (^+.9)* 
This  suggests  that  the  distribution  p.  (r)  is  the  same  for  both  uniform 

flow  and  source  flow  if  the  portion  of  the  flow  field  near  the  shock  wave 
is  excluded  from  consideration. 

This  conjecture  is  in  fact  supported  by  the  method-of-characteristics 

calculations  of  Inouye  (1966b),  whose  unpublished  results  have  been  used 

to  obtain  the  ratio  (p.  )  /(p.  )  given  in  Fig.  21.  Although  the  ratio 

tends  to  behave  erratically  near  the  shock  wave,  it  is  nevertheless  close 

to  unity  for  a  considerable  portion  of  the  wake.  If  the  difference 

[(p+  ),."(P+  )  I/Cp*  )  is  "to  l®ss  than  5^6  for  the  case  under  con- 
”2  ^  ^2  ^  ^2  ^ 

sideration,  then  the  total-pressure  distribution  may  be  assumed  to  be  the 
same  up  to  r  of  about  0.7*  Within  this  range  the  assumption  is  quite 
good,  considering  the  fact  that  the  corresponding  dilference  in  the  wake 
Pitot  pressure  can  be  as  large  eis  50^. 

(2)  The  normalized  static-pressure  distribution  calcvilated 

e 

from  Inouye 's  method-of -characteristics  results,  is  compared  in  Fig.  22 
for  uniform  flow  and  source  flow.  The  agreement  in  this  case  is  good 
all  the  weiy  to  the  shock  wave,  although  one  cannot  take  full  advantage 
of  this  fact,  in  view  of  the  limited  range  of  validity  of  the  assumption 
concerning  the  identity  of  totsJ.  pressure  in  step  (l).  Shown  in  the  scune 
figure  are  points  calculated  from  Meyer's  second-order  blast-wave  result 
of  Eq.  (J+.7)>  for  both  un?'"  *  'low  and  source  flow,  and  from  the  first- 
order  blast-wave  result  of  Eq.  (4.l6).  There  is  good  egreement  among 
the  five  sets  of  results,  except  for  the  second-order  results  for  source 
flow  in  the  range  of  r  between  O.85  and  1.0,  which  is  in  emy  case 
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outside  the  reutige  of  interest.  For  the  poi^tiou  of  the  vake  vhere  r  is 
less  than  0.7^  in  particular,  there  is  only  a  sm}J.  difference  among  all 
five  sets  of  results.  For  sinplicity  of  calculation  it  vill  therefore 
be  assumed  that  the  distribution  given  by  f  j.rst-order  blast-vave  theory 
C8UI  be  used  for  both  uniform  flow  and  source  flow.  The  advantage  here 
is  that  this  result,  Eq,.  (4.l6),  involves  a  simple  aneOytlcal  expression 
that  is  a  function  only  of  r  and  7.  The  second-order  results  of 
Meyer  would  have  improved  the  acc\iracy  somewhat.  They  were  not  used, 
however,  because  of  the  added  complexity  of  calculatioii,  expecially  the 
inclusion  of  the  function  i|f  that  caimot  be  expressed  rnalytically. 

(5/  We  now  consider  the  available  solutions  or  the  longitudinal  variation 
of  the  surface  pressure  Pg  .  Result^  from  the  first-  and  second-orier 

blast-wave  theory,  Eqs.  (4.15)  and  (4.14)  respectively,  are  compared  in 
Fig.  23  SLgalnst  the  more  accurate  numerical  results  of  Inouye.  The 
corresponding  correction  c\jrves  are  .  ven  in  Fig.  c  4.  The  ..  rrection 
formula  of  Hall,  Eq.  (4.17),  bas  been  used  with  the  first-order  blast- 
wave  curve,  even  though  in  Hall’s  aneilysis  he  uses  the  Newton-Busemann 
pressure  formula  rather  than  the  blast -wave  expression  of  Eq.  (4.15). 

This  is  consistent  with  Hall's  theory,  however,  jince  in  deriving  the 
result  of  Eq.  (4.17)  for  the  pressure  ratio  Hall  implies  the  adoption  of 
the  first-order  theory  by  using  the  parabolic  bow  shock  wave  from  this 
theory. 

The  method-of -characteristics  results  in  Figs.  25  and  24  have  been 

obtained  by  extrapolation  and  interpolation  of  Inouye 's  results.  The 

calculations  of  Inouye  (1966a)  were  carried  out  rally  to  x  of  5«  The 

curve  for  urlform  flow  for  points  beyond  is  therefore  extrapolated  with 

the  use  of  the  blast-wa/e  correlation  method  suggested  by  Van  Rise  (1961). 

In  this  extrapolation,  we  have  used  as  a  g>iide  the  method-of -characteristics 

calculations  of  Feldman  (i960)  for  air  with  real-gas  effects  (J^  =  I8) 

n 

and  carrif  1  to  x  =  50.  Also,  since  Inouye' s  results  sure  for 
r^/L^  *  0,  0.02,  0.04,  and  O.O6,  the  correction  curve,  which  is  for 
r^/L^  =  0.01  (or  =  0.8),  is  obtained  by  graphicsil  Interpolation  of 
Inouye 's  correction  curves  (his  Fig.  Ub)  up  to  x/r^  ®  10.  Beyond  this 


point  a  linear  extrapolation  has  been  adopted  to  continue  the  linear 

trend  exhibit  -d  by  the  original  curve. 

The  coniparison  of  results  in  Fig.  25  shows  that  for  the  range  of 

stations  involved  in  this  study  (x  <  15 first-order  blast-vave  theory 

gives  a  more  accurate  prediction  of  the  surface  pressure  than  the  second- 

order  theory  for  both  uniform  flow  and  source  flow.  In  addition,  the 

correction  formula  of  Ha.ll  is  considerably  more  accurate  than  Meyer's 

second-order  result  from  Sq..  (h.lh).  Strictly  speaking,  the  numerical 

results  of  Inouye  are  superior  to  both  Meyer  and  Hall  in  terms  of  accuracy. 

For  the  present  study,  however,  ve  are  interested  in  a  set  of  relativels' 

simple  aneLLytlcal  expressions  that  afford  more  flexibility  vith  regard  to 

the  source-flow  parameter,  the  Mach  number  ,  and  the  longitudinal 

n 

station  x.  These  conditions  ao'e  met  by  the  combined  results  of  Hall, 

Eq.  (4.17),  and  first-order  blast-wave  theory,  Eq.  (4.15),  and  these  will 
be  used  in  the  c\irrent  work,  even  though  some  accuracy  is  sacrificed  in 
the  process. 

(6)  With  regard  to  the  lateral  distance  r  in  source  flow,  we  correct 
it  to  uniform  flow  according  to  the  simple  relation 


u'  s'  s 


(4.19) 


where  the  ratio  of  the  distances  to  the  shock  wave,  is  given  by 

Eq.  (4.l8).  This  correction  appeeurs  to  be  reasonable,  since  the  normalized 

A  . 

radial  distance  r  =  r/R  correlates  quite  well  the  pressure  results  for 
uniform  flow  and  source  flow.  For  the  ratio  ne  have  used  Eq. 

(4.l8)  rather  than  Meyer's  result  from  second-order  blast-wave  theory,  which 
can  be  calculated  from  Eq.  (4.9)*  This  is  because  Eq.  (4.l6)  has  been 
shown  by  Hall  to  be  in  better  agsreement  with  experimental  data. 

This  conpletes  the  discussion  of  the  essentials  of  the  correction 
scheme  for  the  Pitot  pressure  in  the  Inviscid  wake.  Some  of  the  details 
of  the  method  will  be  considered  in  Section  5* 3,  following  a  generaLL 
description  of  the  data- rediict ion  procedure.  We  proceed  now  to  a  discussion 
of  the  corresponding  scheme  for  the  viscous  wake. 
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4.6  Souroe-Flow  Correction  for  Vircous  Wake 


The  general  character  of  the  viscous  vake  has  been  discussed  in 
Section  4.4.  It  vas  pointed  out  that  this  portion  of  the  wake  originates 
in  the  boundary  layer  over  the  hemispherical  nose.  This  boiandary  layer 
upon  separation  becomes  a  free  shesur  layer  that  subseq.uently  coalesces 
to  form  the  viscous  weke.  To  consider  t}"9  source-flow  ?fects  on  the 
development  of  the  viscoiis  wake  —  and  in  particular  on  the  Pitot  pressure  — 
it  is  necessary  first  to  study  the  effect  on  the  flow  in  the  boundary 
layer  and  the  free  shea:?  layer.  In  the  analysis  to  follow  it  will  be 
useful  to  divide  the  viscous-flow  region  into  two  parts;  (l)  the  boiandary 
lai'er  on  the  body  together  with  the  free  shear  layer,  and  (2)  the  viscous 
wake  downstream  of  the  neck.  This  division  is  desirable  because  the 
flow  differs  somewhat  in  these  two  subregions  and  different  approaches 
eu?e  used  to  study  the  source-flow  effects. 

Boxindaiy  Layer  and  Free  Shear  Layer;  We  first  consider  only  the  bc\mdary 
layer  over  the  hemispherical  nose  prior  to  separation.  At  the  end  of 
this  subsection,  it  will  be  shown  that  the  conclusion  reached  for  the 
boundary  layer  applies  equally  well  through  the  free  shear  layer. 

To  study  the  effect  on  the  boundary  layer,  we  use  the  source-flow 
results  of  Inouye  (1966c,)  for  the  inviscid  external  flow  as  an  input  to 
the  boundary-layer  solution  of  Cohen  and  Reshotko  (1956).  Parts  of 
both  these  solutions  that  are  pertinent  to  the  present  study  will  now 
be  described  in  brief. 

The  work  of  Cohen  and  Reshotko  gives  self -similar  solutions  for 

two-dimensional  boundary-layer  flows  with  arbitrary  pressure  gradient 

euid  Mach  number  but  with  the  Prandtl  number  Pr  restricted  to  1.0. 

The  study  also  has  as  a  parametric  variable  the  quantity  S  =  h  /h.  -  1, 

ex 

where  the  subscript  "w”  refers  to  conditions  at  the  surface  of  the 

model.  For  a  perfect  gas  S  is  equal  to  (T  /T.  -  l)  and  thus  is 

^  ex 

a  measxuTe  of  the  ratio  of  model  ten5>erature  to  free-stream  total  ten5>era- 
ture. 

The  inclusion  of  a  non -zero  pressure  gradient  means  that  the  solu¬ 
tion  if.  not  restricted  to  a  flat  plate  and  can  be  applied  to  the  case 


of  the  hemispherical  nose  londer  study  in  the  present  vork,  provided 
that  local  similarity  is  assumed.  Also,  the  two-dimensional  boundary- 
layer  transformations  used  by  Cohen  and  Reshojko  must  be  modified  to  account 
for  the  axial  symmetry  in  the  flow  over  the  hemisphere.  The  restriction 
on  the  Prandtl  number  is  not  a  serious  one  for  the  current  study  in  which 
only  comparative  information  is  being  soioght.  The  ratio  of  the  thick¬ 
nesses  of  the  velocity  and  thermal  boundary  layers  is  associated  with 
the  Prandtl  number  (see  Schlichting,  i960,  p.  30l).  It  is  reasonable  to 
assume  that  the  source  flow  ai'fects  these  layers  to  about  the  same  extent, 
regardless  of  the  ratio  of  the  thicknesses.  We  shall  therefore  use  the 
results  for  Pr  =  1.0  for  the  present  case  where  the  Prandtl  number  is 
approximately  0.7.  As  regards  the  parameter  S  ,  it  is  assinned  to  have 

W 

the  value  -1.0,  that  is,  the  body  is  assumed  to  be  "cold”.  This  is 

Justifiable  because  the  initial  teii5)erature  of  the  model  does  not  exceed 

one-fifth,  and  for  most  runs  is  closer  to  one-tenth,  the  value  of  the 

free-stream  total  temperature.  Sine-'  the  duration  of  the  run  is  “hoar'-, 

the  model  does  not  heat  up  significantly  in  the  course  of  the  run. 

The  primary  input  to  the  solution  of  Cohen  and  Reshotko  is  the 

static  pressure  in  the  inviscid  stream  external  to  the  boundary  layer, 

or  equivalently,  the  surface  pressure  pg  on  the  body  of  interest. 

‘^ex 

For  the  case  of  the  hemispherical  nose  under  study  here,  this  information 

can  be  ootalnsd  from  Inouye’s  surface-pressure  results,  which  eore  gJven 

in  Fig.  25  for  uniform  flow  =  O)  and  source  flow  =  0.02). 

According  to  Inouye,  although  in  the  subsonic -transonic  region  the 

difference  [(p^  )  -(po  )„V(Po  )„  is  of  the  same  order  of  magnitude 

ex  ex  ex 

as  the  source-flow  parameter  the  difference  in  the  supersonic 

region  soon  becomes  an  order  of  magnitude  larger  than  For 

exan^ple,  the  results  of  Fig.  25  show  that  at  9  =  Jt/2  (supersonic  region) 

the  pressinre  ratio  Po  /(p*.  )  ^or  uniform  flow  and  source  flo\-7  is 

‘^ex  ^2  " 

respectively  0.056  and  0.050,  giving  a  difference  of  about  1??6. 

With  thio  large  difference  in  Pg  ,  it  is  reasonabxe  to  expect  that  the 

ex 

development  of  the  boundary  layer  will  be  significantly  affected  by  the 
source  flow.  The  magnitude  of  the  effect  can  be  assessed  by  substituting 
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both  sets  of  surface -pres  sure  resvilts  into  the  solution  of  Cohen  and 
Reshotko  and  couputlng  the  boundary-layer  growth  in  the  two  cases. 

In  contrast  to  the  relatively  large  difference  in  surface  pressure, 
bhe  total  pressure  along  the  bodj-  streaiallne  is  the  saiae  for  uniform 
flow  and  source  flow.  This  results  from  the  ar- sumption  of  issntropic 
flow  in  the  invltcid  stream  in  both  cases.  In  either  case  the  total 
pressure  corresponds  to  the  suagnat Ion-point  pressure,  which  is  the  same 
in  uniform  flow  and  source  flow  since  the  fluid  has  traversed  the  normal 
portion  of  the  bow  shock  wave  and  conditions  ahead  of  it  sure  identical 
in  both  cases. 

On  the  bSLsis  of  thrsse  considerations  of  siurface  pressure  sind  total 
press\are  in  tne  Inviscid  stream,  we  sure  now  in  a  position  to  study  the 
effect  of  the  source  flow  on  the  boundsury  layer.  The  details  of  the 
calculations  and  the  vsuriouF  assumptions  are  given  in  Appendix  D.  Here 
we  skSLll  merely  outline  the  approach  and  summarize  the  resuluc: 

(1)  In  the  solution  of  Cohen  sind  Reshotko,  sui  importsint  psurametrlc 
•'airiable  is  the  preesure-gradient  psurameter  P,  which  depends  on  cotidi- 
tir's  in  the  ex.^ma]  flow.  In  the  present  case,  by  virtue  of  the 
isentropic  condition,  the  external  flow  is  corapletely  specified  by  the 
pressure  I’esults  of  Fig.  2$,  It  is  thus  possible  to  use  these  data  to 
ca3.c\jlate  the  values  of  P  at  9  ®  for  uniform  flow  and  source  flow. 
It  is  found  that  the  difference  (p  -p  )/p  is  only  9.0066  for  the 
case  with  r^/l^  =  0.02.  The  ’'ariation  ii:  p  is  thus  negligibly  small, 
particularly  when  compared  with  the  relatively  large  difference  of  0.17 
in  the  surface  pressure. 

(2)  Certain  ratios  of  flow  quantities,  when  expressed  as  functions  ux 
a  boundary-layer  «;■  a?JLarity  variable  tj,  are  relatively  insensitive  to 
changes  in  p.  Since  the  difference  in  P  between  uniform  flow  and 
source  flow  is  small,  one  can  assume  that  these  ratios  are  the  same  for 
the  two  flows.  Examples  of  ratios  that  satisfy  this  condition  are  the 
total-enthalpy  ratio  b>(./hj.  and  the  velocity  ratio  V/V^^. 

(5)  Using  the  identity  in  the  ratios  of  velocity  and  total  enthalpy, 
one  can  calculate  the  temperature  ratios  and  T^/T^  as 


H 


functions  of  ij.  For  the  case  vith  =  0.02,  the  fractional 

difference  (T  -T  )/T  at  <p  =  is  about  0.05  throughout  almost 
the  entire  boundary  layer. 

(4)  The  temperature  ratios  calculated  in  (3)  can  be  used  to  estimate 
the  difference  In  the  lateral  extent  of  the  boundary  layer  in  the  physical 
plane.  At  cp  =  n/2,  the  ratio  6  /b  of  the  boundary-layer  thickness 
turns  out  to  be  given  by  the  relation 


8^/8^  =  +  4.75  (r^/L^)  ,  for  V  =  Jt/2  . 


(4.20) 


The  above  results  have  been  calculated  for  the  boundaiy  layer  at 
cp  =  jt/2.  It  is  assumed  that  the  boundary  layer  has  not  separated  up 
to  this  point,  but  that  it  does  so  shortly  thereafter,  so  that  the  results 
describe  the  conditions  Just  before  separation.  The  process  of  fluid 
septu  >tlon  is  complicated,  and  no  attempt  vlU  be  made  to  anedyze  it 
here.  We  shall  sliiQ)ly  assimie  that  neither  the  location  of  separation  nor 
the  local  separation  mechanism  is  significantly  different  for  imlform 
flow  and  source  flow.  Whereas  the  profile  shapes  of  specific  flow 
quantities  before  and  after  separation  undoubtedly  differ  considerably 
from  each  other^  it  seems  reasonable  to  assinne  that  the  functional 
relationship  in  each  case  is  independent  of  the  slightly  different  local 
conditions  that  exist  in  \mlform  flow  and  source  flow.  This  functional 
relationship  is  therefore  assu’’  to  be  the  same  for  the  two  flows.  With 
this  assumption,  the  boundary-layer  thickness  8  in  uniform  flow  and 
source  flow  bears  the  seme  relationship  to  each  other  after  separation 
as  before.  Equation  (4.20)  then  provides  the  ratio  of  the  thickness  of 
the  free  shear  layer  Just  after  separation. 

The  application  of  Eq.  (4.20)  can  be  extended  to  the  entire  free 
shear  layer  through  the  following  argument.  The  fluid  in  the  free  shear 
layer  (from  separation  to  recompression  at  the  neck)  is  dominated  by  the 
adjoining  recirculation  zone,  in  which  the  static  pressure,  or  base 
pressure,  is  essentially  constant.  The  free  shear  lEiyer  thus  undergoes 
a  constant-pressure  mixing  process.  The  assumption  of  a  constant-pressure 
process  has  in  fact  been  used  in  base-flow  and  near-wake-flow  models 
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by  various  investigators,  such  as  Denison  and  Baum  (1965)  and  Reeves 
and  Lees  (196^).  To  the  extent  that  this  assumption  is  valid,  it  is 
reasonable  to  stipulate  ths.t  the  ratio  of  the  cross-sectional  area  of  the 
shear  layer  at  separation  to  its  area  at  the  neck  is  the  same  for  both 
uniform  flow  and  soiirce  flow.  We  base  this  assumption  (regarding  the 
ratio  of  areas)  on  the  fact  that  significant  area  changes  in  the  shear 
layer  occur  only  if  there  is  a  marked  variaticn  in  the  pressure  field, 
as  for  example  in  an  expanding  flow.  In  the  present  Instance  of  the 
shear  layer  along  the  recirculation  zone,  however,  the  pressure  is 
essentially  constant  in  eaci^  of  the  two  cases  of  uniform  flow  and  source 
flow.  The  cross-sectional  area  of  the  layer  can  thus  be  expected  to 
remain  basically  unchanged  from  the  separation  point  to  the  neck.  When 
mixing  la  included,  the  shear  layer  does  vindergo  an  area  Increase,  owing 
to  the  entrainment  of  fluid  from  the  outside.  It  seems  reasonable  to 
assvime,  however,  that  the  rate  of  area  Increase  is  the  same  for  both 
uniform  flow  and  source  flow.  This  then  leads  to  the  identity  in  area 
ratio  for  the  two  cases.  On  this  basis  Eq,.  (k.20)  is  applicable  through¬ 
out  the  entire  free  ahetr*  layer  as  well. 

The  boundary  layer  and  the  free  shear  layer  are  annular  in  shape 
at  separation.  As  the  shear  layer  flows  towards  the  axis,  its  thickness 
Increases,  until  eventually  it  coalesces  Just  ahead  of  the  neck.  This 
change  in  cross-sect l<>.ial  geometry  brings  about  a  modification  of  the 
ratio  of  the  lateral  dimension  of  the  viscous  wake  at  the  neck.  The 
ratio  of  the  areas  of  the  annular  boundary  layer  at  cp  =  jt/2  is 


(2r,+6  )6 
'  b  s'  s  ~ 


(4.21) 


for  (2r. +6  )  »  16  -6  ).  Since  the  ratio  of  the  areas  of  the  boundary 
'  b  u  s  u 

layer  is  assumed  to  apply  as  well  to  the  areas  of  the  coalesced  shear 
layer  at  the  neck,  then  we  have 


o 


1  + 


o _ o 

>u 

o 


('+.22) 


vhere  r_  =  r_  /d  denotes  the  radius  of  the  viscous  vake  at  the  neck, 
o 

Thus,  from  purely  geometric  considerations,  the  effect  of  the  soxirce 
liow  on  the  lateral  dimension  of  the  viscous  vake  is  approximately  half 
the  effect  on  the  boundary  layer. 

We  are  now  in  a  position  to  consider  the  viscous  vaJce  v  .■er  down- 
strejtm  of  the  neck. 


Viscous  Vako  Proper;  The  problem  of  deriving  a  source-flow  correction 
for  the  V..SCOUS  wake  will  be  considered  in  two  parts.  The  first  is  to 
determine  the  extent  to  which  the  source  flow  influences  wake  growth. 

Prom  this  consideration  we  devise  a  correction  for  the  lateral  dimension 
of  the  wake.  The  second  aspect  of  the  problem  is  to  obtain  a  local 
magnitude  correction  for  the  measured  Pitot  pressure. 

The  growth  of  the  viscous  wake  is  governed  both  by  expansion  due  to 
the  Imposed  inviscid-pressiire  envelope  and  by  laminar  or  turbulent 
diffusion  of  various  flow  q,\iantities  from  the  inner  to  the  outer  wake. 
Since  one  of  the  major  points  of  difference  between  the  two  portions  of 
the  wake  is  that  th*^  inner  wake  is  much  hotter,  an  Important  fluid 
property  to  be  considered  is  the  static  entheilpy.  Lees  and  Hromas  (I962) 
have  studied  the  ^.roblem  of  wake  growth  on  a  semi-empirical  basis,  taking 
into  account  both  expansion  and  turbulent  diffusion  of  static  en.halpy. 
They  assume  that  the  development  of  the  viscous  wake  can  be  divided 
approximately  into  three  stages  for  which  either  expansion  or  diffusion, 
or  both,  is  Important.  Expansion  without  diffusion  governs  the  initial 
stage  of  wake  growth  up  to  the  point  at  which  the  imposed  inviscid 
pressure  —  and  hence  the  pressure  all  the  way  through  the  viscous  wake  — 
is  four  times  the  free-stream  static  pressure,  p^^  ,  at  the  nose  of 
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the  model.  Downstream  from  this  point  to  where  the  pressure  has  decayed 

to  p^  ,  both  expansion  and  diffusion  are  inportant.  Beyond  this  region 
n 

diffusion  becomes  the  sole  mechanism.  {The  term  "diffusion'’  here  refers 

to  turbulent  diffusion  only,  since  the  laminar  counterpart  may  be  assimied 

to  be  insignificant  by  conparlson. ) 

The  farthest  station  encountered  in  the  present  experiments 

(x  ^  1^)  is  about  midway  through  the  second  region  considered  by  Lees 

auid  Kromas.  According  to  their  model,  therefore,  both  expansion  and 

diffusion  should  be  included  in  devising  a  source-flcw  correction  scheme. 

When  their  results  are  considered  in  detail,  however,  it  is  found  that 

within  the  present  region  of  Interest  (that  is,  up  to  x  *  15)  the 

contribution  to  wake  growth  from  diffusion  is  relatively  small.  By 

assuming  that  the  expansion  and  diffusion  mechanisms  can  be  treated  as 

vincoupled,  one  can  devise  a  correction  scheme  that  accounts  for  source- 

flow  effects  due  to  expansion  alone.  The  diffusion  contribution  by 

virtue  of  its  Insignificant  part  in  the  overall  pictvire  of  wake  growth, 

is  sissimied  to  be  the  same  in  both  tiniform  and  source  flows. 

We  now  consider  in  some  detail  parts  of  the  solution  of  Lees  and 

Hromeie,  who  have  worked  out  the  uniform-flow  problem  for  M.  s  22,  At 

■^■n 

this  tiach  number,  the  'rrake  presstire  is  4p^  at  x  »  cuad  approaches 

n 

Pj^  at  X  of  about  50.  The  turbulent-diffusion  calculations  were 
n 

carried  out  for  two  cases;  (l)  presst’re  held  constant  at  p^  down- 

n 

stream  of  x  =  7.5,  and  (2)  pressure  gradually  decreasing  from  4p^ 

n 

to  p^  ,  starting  from  the  same  axial  location.  The  results  of  their 

calculations  are  shown  in  Fig.  26  ,  in  which  r^  =  r^/d  is  the  normalized 

radius  of  the  viscous  wake,  the  subscript  "i”  denoting  the  value  at 

X  =  7.5,  which  corresponds  to  the  station  at  which  diffusion  is  first 

Included.  In  case  (l),  ?_/?_  remains  less  than  1.2  up  to  x  of 

r 

*  The  c\urves  in  this  figure  sure  deduced  from  Fig.  of  Lykoudis  (1965) 
rather  than  from  Lees  and  Hromas  because  the  scales  used  in  Lykoudis 's 
work  60%  easier  to  read. 
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about  25,  vhereas  the  growth  is  considerably  faster  over  the  same 
distance  when  expansion  is  taken  into  account  in  case  (2).  The  relatively 
small  wake  growth  of  case  (l)  may  be  interpreted  as  indirect  evidence 
that  diffusion  does  not  contribute  significantly  to  the  overall  growth 
upstream  of  x  =  25- 

In  the  application  of  the  foregoing  results  -  worked  out  for 

=22  —  to  the  present  case  of  =  17,  the  longitudinal  distances 
n  n 

at  which  the  pressure  is  4p-  and  p  are  now  different.  To  cedculate 

n  n 

the  pertinent  distances  in  the  present  case  we  can  use  the  pressTore 
formula  from  second-order  blast-wave  theo'  r,  Eq..  (14-.14),  this  being  pre¬ 
ferred  over  the  first-order  expression  of  Eq..  (4.15)  because  the  former 

is  more  acciarate  for  the  larger  distances  Involved  here.  For  =  17, 

n 

it  is  found  that  the  pressure  is  4p^  at  x  =  4. 5  while  p^  occurs 

n  2 

at  X  =  30.  The  corresponding  dlstEmces  for  M-  =  22  were  x  =  7*5 

■^n 

and  X  =  50.  Since  for  this  earlier  case  turbulent  diffusion  was 
relatively  unimportant  upstream  of  about  25  body  diameters,  we  may  assume 
that  in  the  present  Instance  it  is  likewise  uniniportant  for  x  <  I5. 

Hence  in  considering  the  data  from  the  current  experiments,  in  which 
X  varies  from  3  to  15,  we  shall  account  for  the  source -flow  effects 
associated  with  the  expansion  process  but  not  those  connected  with 
dlffiision. 

We  now  consider  the  expansion  process  that  accompanies  the  decay 

of  static  pressiare  along  the  wake.  At  any  given  longitudinal  station, 

the  static  pressiare  within  the  viscous  wake  is  assumed  to  be  uniform 

and  eq,ual  to  the  press\ire  of  the  inviscid  flow  Just  outside  the  wake. 

The  longitudinal  variation  of  this  pressure,  for  both  uniform  flow  and 

soTorce  flow,  is  a  known  q.uantity  either  from  the  numerical  solution 

of  Inouye  or  from  first-order  blast-wave  theory  (see  Section  4.5)*  In 

addition  the  ratio  of  the  areas  of  the  viscous  wake  at  the  neck, 

A  /a  ,  is  given  by  Eq,s.  (4.20)  and  (4.22)  as 
”s  u 


r(f(  )si2 

o 


(?f«  ) 

f  'u. 
o 


=  1  +  4.75 


(4.23) 
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The  problem  Is  now  to  derive  an  expression  for  the  area  ratio  A  /A^ 
for  any  station  farther  downstream.  This  will  then  provide  a  means  for 
correcting  the  lateral  position  within  the  wake  from  source  flow  to 
uniform  flow.  The  expression  for  this  area  ratio  should  be  in  terms 
of  the  static  pressure  and  flow  q,uantlties  deduclble  therefrom. 

To  derive  the  desired  expression,  consider  an  elemental  annular 
streamtube  within  the  viscous  wake.  If  £A  denotes  the  cross-sectioned 
eirea  of  this  streamtube  at  any  station  x  and  £A*  Is  the  fictitious 
sonic  area  associated  with  then  one  cem  write  the  area  relation^ 


In  particular,  for  and  6k* 

o  o 


Just  downstream  of  the  neck  we  have 


o 


(4.24b) 


These  equalities  core  point  relations  euid  do  not  require  the  assumption 
of  Isentroplc  flow.  In  fact  the  flow  within  the  viscous  wake  In  general 
Is  not  Isentroplc,  and  the  total  pressure  and  sonic  area  are  not  invariant 
along  any  given  streamtube.  We  shall  however  make  the  slnplifying 
assumption  that  their  variations  for  uniform  flow  and  source  flow  are 
Identical  insofeo:  eus  ratios  between  corresponding  longitudinal  stations 
are  concerned.  That  is,  we  assume  that  for  a  given  pair  of  longitudinal 
stations  x^  and  x^  within  the  viscous  wake  one  may  write 


In  this  equation  and  in  the  discussion  to  follow,  the  subscript  ”2" 
that  denotes  "downstream  of  the  bow  sihock  wave"  is  omitted  so  as  to 
render  the  subscript  sltviatlon  more  easily  manageable. 
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(4.25) 


(p* /Pt.ls  '  '• 


V 


Thus  local  differences  in  pressure  and  other  thermodynamic  quantities 
(between  \miform  flow  and  source  flow)  are  assianed  to  have  an  insignificant 
effect  on  the  ratios  of  total  pressure  and  sonic  area  of  corresponding 
streamtubes.  If  Eqs.  (4.24)  are  now  written  for  both  uniform  flow  and 
soxurce  flow  to  form  the  ratio  AA^/AA^  and  Eqs.  (4.25)  are  used  to 
simplify  the  resilLt,  one  obtains  the  relation 


u 


{h.26) 


where  the  notation  Q  is  used  to  denote  the  ratio  Q/Q^* 

On  the  right-iiand  side  of  this  equation,  the  ratio  AA^  /AA^  is 

8  u 

equal  to  /A^  ,  which  is  given  by  Eq.  (4.25).  Its  value  is  a 
°s  u 

fxinction  of  the  source-flow  parameter  but  not  of  the  lateral 

location  of  the  streamtube  within  the  viscous  wake.  The  static  pressure 

is  uniform  at  a  particular  longitudinal  station,  so  that  the  ratios  p^ 

and  p  are  also  independent  of  lateral  position.  The  compound  ratio 
s 

in  general  varies  with  the  lateral  coordinate.  We  shall  assume 
however  that  at  corresponding  longitudinal  stations,  the  transverse 
distribution  of  the  normalized  Mach  number  M/M^  is  identical  for  xiniform 
flow  and  source  flow,  that  is, 

(M/M^)g  =  (H^Me^u  • 

The  rationede  behind  this  assxui5)'^ion  will  be  discussed  in  Appendix  E. 
Writing  this  eqxxation  for  and  x,  ve  then  obtain 

u  s 
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can  be  calculated  from  the  inviscid-pressure 


The  compound  ratio  St  /H 

u 

results,  since  at  the  edge  of  the  viscous  vaJce  vhe  flow  is  isentropic. 
When  this  relation  is  used  in  Eq.  (4.26),  the  entire  right-hand  side 
becomes  independent  of  the  lateral  coordinate  of  the  particular  stream- 
tube  under  consideration.  The  elemented  arer\s  can  thus  be  integrated 
to  give 


(4.29) 


Since  the  integration  can  be  carried  out  from  the  wake  axis  to  einy 
desired  position  r  within  the  viscous  wake,  the  ratio  applies 

not  only  to  the  overedl  area  of  the  viscous  wake  but  also  to  the  area 
of  radius  r.  Thus  the  desired  source-flow  correction  formula  for  the 
lateral  position  r  can  be  obtained  from  Eq.  (4.29) 


We  turn  now  to  an  outline  of  the  procedure  for  correcting  the 
magnitude  of  the  measured  Pitot  pressure.  The  scheme  is  analogous  in 
general  to  the  one  \ised  in  Section  4.^  for  the  Invlscld  wake. 

(1)  The  measured  Pitot  pressure  corresponds  to  the  source-  flow  value, 

denoted  by  p^  .  In  perticuleu:,  its  value  at  the  edge  of  the  viscous 

6  _ 

wake,  (p'  )  “  p'  (r^  ),  is  known  from  the  measurements  provided  r 

e  ^8  s  _  8 

is  foimd.  The  determination  of  r  will  be  considered  in  Section  5*5. 

6 

8  _ 

(2)  If  the  flow  at  the  viscous-wake  edge  r^  or  r^  •  is  assumed  to 

u  ®s 

be  isentropic,  the  total  pressvir«  p.  is  equal  to  the  known  stagnation 

% 

pressiire  at  the  nose,  this  pressure  being  the  same  for  both  uniform  flow 
and  source  flow.  The  Justification  for  the  isentropic -flow  assumption 


•will  be  given  following  this  outline 


(5)  For  the  source -flow  case,  the  knowledge  of  (p'  and  (p.  ) 

OS  u  s 

e  e 

determines  uniquely  the  ratios  of  flow  quantities  at  r^  .  In  particular, 

's 

the  quantities  -I  and  p  are  known. 

G  s 

(4)  Since  p  within  the  viscous  wake  is  equal  to  p  ,  its  determi- 

®  ®s 

nation  in  step  (3),  together  with  the  measurement  of  p^  ,  leads  to  the 

specification  of  the  Mach  number  within  the  entire  viscous  wake  in  source 
flow.  Thus  it  is  possible  to  calc\alate  the  quantity  (M/M  )  through- 
out  the  wake. 

(5)  The  major  assumption  that  relates  the  source- flow  results  to  those 

for  uniform  flow  concerns  the  Maeh-number  ratio  l^^e*  ^ 

deri'vation  of  the  formula  for  correcting  the  latereG.  distance,  this 
ratio  is  assumed  to  be  the  same  for  both  uniform  flow  eind  source  flow. 

The  basis  for  this  assumption  is  considered  in  Appendix  E. 

(6)  To  calculate  M  ,  we  observe  from  step  (2)  that  p.  is  the  same 

0  ^ 

U  0 

for  uniform  flow  end  sovirce  flow.  Only  one  other  flow  quantity  is  needed 
to  specify  the  flow,  and  for  this  the  static  pressure  p^  =  p^  is 

chosen.  The  -value  of  p  is  obtained  from  p  by  meems  of  the 

e  e 

u  s 

correction  formula  of  Hall,  Eq.  (4.17)» 

(7)  The  flow  within  the  viscous  wake  in  a  hypothetical  uniform  stream 

is  now  specified  by  p^  and  to  the  extent  that  the  desired  -value-. 

of  p'  can  be  calculated  with  the  use  of  a  routine  normal-shock  relation, 
u 

In  the  foregoing  procedure,  the  isentrqpic-flcw  assumption  is 
invoked  at  the  edge  of  the  viscous  wake.  The  flow  here  is  assumed  to 
have  expanded  isentropically  from  the  stagnation  point.  Strictly  speaking, 
there  are  losses  at  the  separation  point  and  during  recorapression  at  the 
wake  neck.  Although  these  losses  may  not  be  negligible  in  each  of  the 
two  cases  of  imiform  flow  and  source  flow,  we  reason  that  their  effect 
on  the  overall  correction  scheme  can  be  taken  to  be  negligible.  This 
can  be  seen  from  the  following  consideration.  Owing  to  separation  and 
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recompresBlon  losses,  the  actual  total  pressure  at  the  edge  Is  less  than 

th^  stagnation  pressure  and  can  be  represented  by  k,(p.  )  .  Since  the 

Ml  t2  n 

losses  are  not  large  euid  are  con?)arable  for  uniform  flow  and  source  flow, 
the  factor  takes  on  a  value  close  .  1  for  both  cases.  Now,  the 
total  pressure  at  the  edge  is  used  in  steps  (3)  and  (6)  in  a  ’’parallel" 
manner,  that  is,  once  for  uniform  flow  and  once  for  source  flow  in 
identical  sets  of  formulas.  It  would  anpear,  therefore,  that  using 
(p^  in  place  of  k^(p^  and  thereby  neglecting  the  factor  k^ 

in  the  calculations  will  lead  to  only  a  small  error  in  the  overall  results. 

The  analysis  of  the  source-flow  correction  scheme  for  the  viscous 
wake  is  now  complete.  The  application  of  the  results  of  this  and  the 
preceding  section  to  the  experimental  data  will  be  considered  in  Section 
5.3.3,  which  also  presents  some  of  the  intermediate  steps. 
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5.  DATA  REDUCTION  AND  RESULTS 


5 • 1  Preliminary  Remarks 

The  theory  te’  ind  the  method  for  correcting  for  transducer- 
galvanometer  lag  has  been  treatei  in  Chapter  3.  The  corresponding  theory 
for  source-flow  effects  has  been  presented  in  Chapter  4.  Consideration 
will  now  be  given  to  the  problem  of  implementing  these  theories  and 
applying  them  to  the  experimental  data. 

Although  data  from  both  the  vertically  and  axially  traversing  probes 
will  be  considered,  the  correction  will  be  applied  only  to  the  results 
from  the  former.  The  reasons  for  omitting  the  corrections  in  th ,  case  of 
the  axially  traversing  probe  are  discussed  in  Section  5*2. 

We  shall  present  the  results  from  the  two  probes  separately  in 
Sections  5*2  and  5 *3*  In  each  section  the  salient  points  in  the  data- 
reductlon  procedure  are  discussed  in  some  detail.  Because  only  Pitot- 
pressure  data  were  obtained,  it  is  not  possible  to  deauce  additional 
information  regarding  the  wake  flow.  For  this  reason,  no  attempt  will 
be  made  to  discuss  the  experimental  resuJ.ts  in  this  context.  Some  con¬ 
sideration  will  be  given,  however,  to  the  magnitude  of  the  corrections  to 
which  the  experimental  data  were  subjected. 

5*2  Results  from  Axially  Traversing  Probe 
5.2.1  General  Discussion 

A  relatively  small  number  of  runs  were  made  with  the  axially  traversing 
probe.  The  results  from  these  runs  will  be  pi'esented  without  correction 
for  c isponse  lag.  As  to  source-flow  effects,  instead  of  applying  the 
detailed  correction  scheme  of  Chapter  4,  we  adopt  a  simpler  and  more 
intuitive  approach  to  account  for  the  effects  in  this  case. 

Response  corrections  are  considered  unnecessary  since  the  time  rates 
of  change  in  the  pressure  signal  are  substantially  smaller  than  those 
encountered  with  the  vertically  traversing  probe.  Two  factors  contribute 
to  this  effect:  (l)  the  spatial  gradients  in  the  wake  flow  are  less 
severe  in  the  longitudinal  direction  than  in  the  traverse  direction,  and 
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(2)  thrt  speed  of  35  fps  of  the  axially  traversing  prohe  is  small  cotirpared 
vith  the  speed  of  about  60  fps  for  the  vertically  traversing  probe,  so 
th;it  the  spatleil  gradients  do  not  appear  as  rapidly  changing  time  signals 
to  the  transducers.  For  this  reason,  it  vas  thought  that  only  minor 
signal  distortions  were  Involved.  It  was  therefore  considered  unnecessary 
to  correct  for  the  time  lag  due  to  a  slow  galvanometer. 

An  aspect  of  the  resjionse  proble-  that  cannot  be  Ignored,  however, 
is  the  anomalous  behavior  of  some  of  the  wafer  transducers  (see  Appendix 
a).  When  a  step  pressure  of  known  magnitude  is  applied  to  these  trans¬ 
ducers,  they  register  initially  anywhere  from  0.6  to  1.5  of  the  expected 
pressure  Increase.  These  measured  values,  based  on  static  calibration, 
do  not  reflect  the  expected  oscillatory  response  that  is  characteristic 
of  geilvanometers;  that  is.  If  the  galvanometer  overshoot  were  Included, 
the  value  of  1.3  would  be  even  higher.  Thus  the  st.isltlvlty  of  the 
galvai'.ometers  as  determined  by  static  calibration  is  not  applicable  when 
recording  time-'v'arylng  signals. 

In  the  case  of  the  axially  traversing  probe  this  behavior  of  th-* 
transducers  has  led  to  difficvilties  that  ui.d  not  arise  with  the  vertically 
traversing  probe.  In  the  latter  instance,  it  was  possible  to  ascertain 
experimentally  the  response  characteristics  of  the  particular  transducer 
involved,  priiarily  becaut  .•  this  transducer  responded  quite  accurately 
to  the  Imposed  pressiire  step.  The  value  of  k,  the  transducer  correction 
factor,  was  found  uo  be  1.04  (see  Appendix  A).  When  the  value  of  k 
departs  significantly  from  1.0,  however,  the  question  of  recovery  to  the 
asymptotic  value  becomes  a  difficult  problem  for  which  the  response 
study  does  not  provide  a  solution.  The  inconsistent  response  of  the 
transducers  explains  in  part  the  peculiar  shape  of  transverse  Pitot- 
preosure  profiles  that  were  previously  reported  (without  comment)  by 
Alllgood,  Kyser,  and  Tsao  (1963). 

In  view  of  the  serious  degree  of  uncertainty  associated  with  trans¬ 
ducer  performance,  we  consider  it  unwarreuited  to  attempt  a  detailed 
source-flow  correction  of  the  data  frcxn  the  axially  traversing  probe. 

For  the  same  reason  attempt  will  be  made  to  reduce  the  data  obtained 
froiT’  alx  the  Impact  tubes  in  the  rake.  Effort  will  be  concentrated 
.’uistead  on  deducing  the  variation  of  Pitot  pressiure  along  the  axis  of 
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the  vske,  this  heing  the  most  valuable  information  obtainable  from  the 
data  under  these  circumstances.  In  deducing  this  variation,  an  intuitive 
method  of  correction  wiiJ.  be  used  to  account  for  source-flow  effects. 

This  is  discussed  in  the  next  section. 

5.2.2  Data  Reduction  and  Results 

Tunnel  runs  with  the  axially  traversing  probe  were  made  (both  with 
and  without  the  sphere)  for  each  of  the  three  sets  of  operating  conditions 
listed  in  Section  2.7.  In  the  wake  runs,  the  sphere  was  positioned  at 
station  F,  and  the  useful  traverse  of  the  probe  was  from  x  -  4  to 
X  -  9*  Figure  27  shows  the  oscillograph  records  of  two  runs  made  at  the 
same  operating  conditions,  with  and  without  the  model.  The  traces 
corresponding  to  the  impact  tubes  located  at  y  =  ±  0.094  are  marked  on 
the  records.  (The  coordinate  y  denotes  the  normalired  vertical  distance 
with  respect  to  the  measured  tunnel  axis,  values  abor/e  the  axis  being 
positive. )  These  two  tubes  are  the  closest  to  the  idealized  wake  axis. 

In  actuality  the  sphere  falls  sligntly  during  the  run;  if  exact  position 
correspondence  is  desired,  appropriate  adjustment  must  be  made  to  account 
for  the  sphere’s  free  fall. 

Since  the  stagnation  pressvire  on  the  sphere  cannot  be  measured 
directly,  this  inlurmation  is  obtai  .ed  from  a  corresponding  calibration 
run  (see  Section  2.7  for  details).  Ordinarily  this  stagnation  pressure 
is  used  to  normalize  the  wake  Pitot  pressure  so  as  to  account  for  the 
gradual  decay  in  the  t\mnel  flow.  In  the  present  instance,  we  are  inter¬ 
ested  also  in  making  use  of  any  readily  available  means  of  accounting 
for  source -flow  effects  in  the  timnel,  other  than  making  direct  use  of 
the  detailed  correction  scheme  devised  in  Chapter  4.  How,  the  tunnel- 
eii5)ty  traces  of  Fig.  27(b)  give  the  longitudinal  \'arxation  of  the  Pitot 
pressure  in  the  free  stream  of  the  tunnel  source  flow.  An  intuitive 
approach  is  therefore  to  use  this  free-stream  variation  to  correct  for 
the  source -flow  effects.  Given  tlriat  this  is  the  only  information  avail¬ 
able,  we  can  show  that  both  tunnel  decay  and  source-flow  effects  can  best 
be  taken  into  accoimt  if  the  wake  Pitot  pressure  is  simply  divided  by 
the  corresponding  free-stream  Pitot  pressure  at  the  same  locati  :i.  This 
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relatively  stralghtfojrward  data-reduction  procedure  vas  followed  in 
obtaining  the  longitudinal  variation  of  the  Pitot  pressure  near  the  wake 
axis  (r  =  O.l).  The  results  of  three  runs,  correspo.'  '.ing  to  the  three 
basic  Reynolds  numbers  for  the  experiments,  are  presented  in  Fig.  28. 


5*3  Results  from  Vertically  Traversing  Probe 
^.3<1  Data  Reduction 

Figvire  29  shows  the  oscillograph  record  for  a  wake  run  made  with 
the  vertically  traversing  probe.  The  traces  that  are  labelled  correspond, 
from  top  to  bottom,  to  the  arc-chamber  pressure,  two  channels  of  wake 
Pftot  pressure  measured  with  a  stationary  probe  (see  below),  photocell 
output,  th';  wake  Pitot-pressure  profile  and  the  associated  mechanical 
noise  measured  with  the  ti'aversing-probe  transducers,  and  the  pulsed 
signal  that  gives  the  instant  at  which  the  sphere  is  released  into  the 
stream.  The  function  of  the  stationery  probe,  whose  position  in  the  wake 

is  shown  in  the  accompanying 

Tunnel  axis 

sketch,  is  to  monitor  the 
random  fluctuations  in  the 
tunnel  flow.  This  is  done 
by  using  the  Pitot  pressure 
from  the  stationary  probe 
in  conjunction  with  the 
stagnation  pressure  taken 
from  a  calibration  run 
(refer  to  Section  2.7  for 
the  method  of  obtaining  the 
latter  pressure).  Ey  assuming 
that  the  ratio  of  these  two 
pressures  is  essentially 
constant  during  the  run,  we 
can  incorporate  the  random 
fluctuations  tliat  are  present  in  the  woke  Pitot  press  re  into  the  stag¬ 
nation  pressure.  The  stagnation  pressure  thus  modified  is  then  used  to 
normalize  the  wake  Pitot-pressure  profile. 


2.25" 


Top  View 


r 


Tunnel  axis 
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The  photocell  output  gives  the  position  of  the  prohe  with  respect 
to  the  tunnel  during  probe  traverse.  Before  release  the  sphere  is 
positioned  at  the  center  of  the  timnel  (to  witv.in  about  0.05  inches); 
after  its  release  it  is  assumed  to  fall  from  gravity  effects.  Thus  with 
an  appropriate  correction  for  free  fall,  the  photocell  output  then  gives 
the  position  of  the  probe  with  respect  to  the  sphere,  if  for  the  moment 
we  assume  that  there  is  no  response  lag  in  the  system  associatea  with 
the  photocell.  To  proceed  from  this  to  a  calculation  of  the  lateral- 
distance  coordinate  for  the  Pitot-pressure  profil  the  oscillograph 
trace  must  first  be  corrected  for  transducer-galvanometer  response  lag. 
(This  correction  will  be  discussed  in  Section  5«5*2.) 

After  the  gravity  and  response  corrections  had  been  applied,  it  was 
foxmd  that  the  axis  of  the  wake,  taken  to  be  the  minimum  of  the  Pitot- 
pressure  profile,  is  at  a  calculated  position  of  y  of  approximately 
+0.2  inches.  This  discrepancy  occurred  consistently  throughout  the 
reduced  data.  A  probable  source  of  error  is  that  some  response  lag  does 
exist  in  the  photocell  galvanometer,  even  though  its  response  is  relatively 
fast  (flat  to  1,000  cps).  We  assume  that  it  is  permissible  at  this  point 
to  shift  the  entire  profile  unifoimily  so  that  the  minimum  corresponds  to 
y  =  0.  After  this  shift  the  lateral  dimension  of  the  wake  is  then 
finally  corrected  for  source-flow  effects  (see  Section  5»5»5)» 

5.3«2  Response  Correction 

In  Section  3.6  an  example  was  given  of  the  application  of  the 
correction  scheme  to  a  specially  constructed  input-output  pair  of  curves. 

In  that  instance,  we  were  able  to  recover  accurately  the  input  from  the 
output  through  a  senes  of  successive  approximations.  When  the  correction 
scheme  was  applied  to  the  measurea  Pitot -pres sure  profiles,  however,  some 
difficulty  was  encountered  at  first  in  deducing  the  correct  input  for  a 
given  recorded  output.  If  the  d  >.ta  that  have  been  subjected  to  mechanical- 
noise  correction  and  subsequent  normalization  by  the  stagnation  pressure 
are  used  directly,  the  input  ciurves  deduced  from  the  correction  scheme 
exhibit  large  oscillations  throughout  the  entire  curve. 

An  examination  of  the  output  data  shows  tliat  the  points  do  not  form 
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a  smooth  c\irve,  owing  to  the  fact  that  it  is  not  possible  to  eliminate 
the  mechanical  noise  completely  through  the  use  of  the  accelerometer  out¬ 
put.  A  factor  that  also  contributes  to  the  unevenness  is  that  the  random 
fluctviations  in  the  tunnel  flow  cannot  be  tedcen  out  congjietely  by  normal¬ 
ization.  Since  the  galvanometer  functions  basically  as  a  low-pass  filter, 
a  small  oscillation  in  the  output  is  interpreted  by  the  correction  scheme 
as  having  origiziated  from  a  much  larger  oscillation.  In  addition,  the 
input-output  relationship  is  not  unique;  that  is,  it  is  possible  for 
more  than  one  input  to  have  the  identical  output.  This  aspect  of  the 
system,  which  is  reflected  in  the  correction  scheme,  apparently  leads  to 
instability  problems  in  the  iteration  procedure  once  oscillation  sets  in. 

To  remedy  the  problem  of  exaggerated  noise  in  the  input,  the  output 
curve  is  first  smoothed  before  it  is  subjected  to  the  correction  scheme. 
The  cvtrve- smoothing  formulas  used  are  steuidard  least-squares  formulas 
(see  Hildebrand,  1956,  Chapter  especially  pp.  29I-296).  The  degree 
of  the  polynomial  and  the  number  of  points  employed  are  dependent  on  the 
shape  of  the  particular  portion  of  the  curve  under  consideration,  and  on 
how  much  smoothing  is  required.  The  curve-smoothing  routine  that  was 
eventually  incoiporated  in  the  data-reduction  program  consists  of  third- 
degree  formulas  that  apply  to  five,  seven,  or  nine  points.  Application 
of  the  formulas  is  repeated  \intil  the  deduced  input  becomes  free  of  un¬ 
wanted  oscillations.  Apart  from  the  foregoing  modification  00  the 
recorded  output,  however,  the  basic  correction  scheme  is  identical  to 
the  method  outlined  in  Sectlcxi  5.6. 

5.5.3  Source-Flow  Correction 

In  the  theory  of  source-flow  correction  presented  in  Chapter  !«■, 
perfect-^s  relations  were  used  throughout  fche  derivation.  Strictly 
speaking,  the  perfect-gas  assumption  is  at  times  incorrect,  since  the 
temperature  In  parts  of  the  wake  is  sufficiently  high  that  real-gas 
effects  must  be  taken  Into  consideration.  Although  the  temperature  dis¬ 
tribution  in  the  wake  is  not  given  a  priori,  it  is  nevertheless  known 
that  the  viscous  wake  is  notably  l">t,  particularly  in  the  initial  stages. 
The  viscous-wake  temperature  is  likely  to  be  of  the  order  of  the 
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free-stream  total  ten^jcrature,  ■which  ranges  approximately  from  1,250 

to  5,T50®K  in  the  experiments  reported  here  (see  Section  2.7)*  In  this 

range  of  temperatures  in  nitrogen,  whose  characteristic  temperature  for 

vibration  is  5,590®K,  the  vibrational  contributions  to  the  ■values  of  the 

specific  heats  is  considerable.  For  example,  at  a  temperature  of  2,850®K, 

c  is  0.89R.  Thus  the  possible  influence  of  caloric  imperfections 

vib 

in  the  gas  must  be  recognized. 

Since  the  correction  method  itself  is  approximate  in  nature,  no 
attempt  has  been  made  to  use  exact  calculations  for  flow  properties. 

Instead  the  alternative  approach  is  adopted  whereby  it  is  shown  that 
C'-iloric  imperfections  liave  negligible  effect  on  the  corrected  results. 

To  show  this  we  first  make  the  engineering  assumption  tl.at  perfect-gas 
formulas  may  be  used  locally  to  relate  flow  quantities.  The  value  of  y 
in  these  formulas  is  •taken  to  be  the  value  appropriate  to  the  local  flow. 
Typical  correction  calculations  are  then  made  for  selected  ■values  of  y 
within  a  sufficiently  inclusive  range.  The  limiting  •values  are  taken 
to  be  1.4  and  1.5,  corresponding  respectively  to  the  ■value  for  a  calorically 
and  thermally  perfect  gas  and  to  the  ■value  of  9/T  (approximately)  when 
the  vibrational  modes  are  fully  excited.  The  results  of  these  calcula¬ 
tions  differ  by  less  than  0.1^.  From  this  we  conclude  that  the  correction 
method  is  insensitive  to  caloric  imperfections,  and  that  the  ■value  of 
1.4  may  be  used  for  y  throughout  the  calculations  without  seriously 
affecting  tne  results. 

With  these  preliminary  considerations,  we  now  proceed  to  the  appli¬ 
cation  of  the  correction  scheme  to  the  results  from  a  typical  run.  We 
begin  by  using  the  experimental  pressure  da'ta  to  deduce  the  lateral 
■variation  of  static  pressure  thro\:ighout  the  wake  (at  the  given  longi¬ 
tudinal  station).  This  reduces  essentially  to  finding  the  pressiire 

(pp  )  at  the  edge  of  the  viscous  wake,  since  in  the  viscous  wake  th'* 
e  ® 

pressure  is  assumed  to  be  cor.  .teuit  and  equal  to  (pp  )  ,  and  the  pressure 

e 

in  the  inviscid  wake  is  given  through  the  ranio  Pp  (r  )/(pp  )  from 

^8  £.8 
8  e 

the  first-order  blast-wave  result  of  Eq.  (4.16). 

It  is  possible  to  deduce  (pp  )  =  Pp  (f  )  from  experimental 

w  S  £  0 

e  8 


can 


pressure  data  provided  the  location  of  the  viscous-vaJte  edge,  r  , 

s 

be  ascertained.  For  the  near  stations  (stations  A,  B,  C)  vhere  the 

viscous  vake  is  laminar,  r^  can  be  found  by  observing  the  change  of 

s 

shape  of  the  Pitot-pressure  profile  in  going  from  the  viscous  vake  to 

the  inviscid  vake.  The  pressure  in  the  laminar  viscous  vake  follows  a 

Gaussian  variation,  vhereas  in  the  inviscid  vake  the  profile  exhibits  a 

parabolic  shape.  The  edge  of  the  viscous  vake,  r  ,  is  taken  to  be 

“s 

the  point  vhere  the  pressure  ceases  to  follov  the  Gaussian  variation. 

•  * 

Thus  for  the  profile  shown  in  Fig.  30,  r^  is  taken  to  be  O.56.  For 

“s 

the  more  distant  stations  (stations  D  through  H),  the  profile  does  not 

exhibit  an  obvious  change  of  shape,  and  it  is  not  possible  to  locate  the 

edge  from  the  mere  examination  of  the  profile.  It  is  then  necessary  to 

assume  that  the  value  of  r  can  be  determined  from  a  linear  extra 

e 

*  s  * 

polation  of  a  plot  of  r^  versus  x  for  the  near  stations  (stations 

s 

A,  B,  C). 

With  Tg  known,  the  Pitot  pressure  (r^  )  is  then  obtained 
s  s 

from  the  met  svired  profile.  The  pressure  (p^  )  at  the  edge  can  now  be 

e  ^ 

deduced  from  p*  (r_  )  and  the  local  total  pressure  [p.  (r  )]  .  This 
.  °s  ®s  ^2  ®s  ® 

latter  quantity  is  equal  "io  the  stagnation  pressure  on  the  model,  by 

virtue  of  the  assvoned  Isentroplc  condition  along  the  streamline  emanating 

from  the  stagnation  point  (see  Section  k,6  for  details).  We  first 

calculate  the  Mach  number  (Mg  )g  from  the  normal-shock  relation  (with 

unnecess6n:y  subscripts  omitted) 


(7+l)t^' 

JL. 

7-1 

r 

_(7-1)I^+2^ 

27^^-(7-l) 

(5.1) 


This  method  of  determination  has  been  used  successfully  by  McCarthy 

(1962). 

This  is  dor.e  by  first  rearranging  Eq.  (5»l)  to  give  an  equation  in 
and  then  solving  for  the  root  by  any  standard  method,  such  as 
Newton's  method,  that  is  applicable  to  algebraic  equations. 
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The  static  pressure  (p„  )  is  then  obtained  from  the  Rayleigh  Pitot 

b  s 

Q 

formula 


K 

’(7+l)^^' 

_2L 

7-1 

7+1 

Pg  “ 

2 

_27m|-(7-1). 

(5.2) 


With  reference  to  the  step-by-step  outline  of  the  correction  schemes 
set  forth  in  Sections  4.5  and  4.6,  we  now  proceed  to  some  of  the  addi¬ 
tional  details. 

Invlscld  Wake;  In  deference  to  the  iimixied  lateral  range  of  validity  of 
the  assun^jtion  concerning  the  identity  in  the  total-pressure  distribution 
for  xiniform  flow  and  source  flow,  the  correction  scheme  is  applied  to 
data  up  to  r  =  0.8  only.  (Since  the  procedural  steps  to  follow  do  not 
necessarily  correspond  to  the  ones  in  Section  4.5,  a  different  numbering 
system  is  used  here. ) 

(a)  The  pressure  ratio  Pg(r)/p2  is  the  same  for  both  flows  and  is 

given  by  the  first-order  blast-wave  result  of  Eq,.  (4.l6).  From  this 

equation  the  static -pressure  distribution  in  source  flow,  Pg  (r),  is 

s 

readily  obtained. 

(b)  The  transverse  distributions  of  static  pressure  and  of  the  measured 
Pitot  pressure  are  combined  to  give  the  total-pressure  distribution 

(p.  )  in  source  flow.  In  these  calculations  we  make  use  of  Eqs.  (5.l) 
and  (5.2),  in  reverse  order,  again  using  the  Mach  number  as  an  inter¬ 
mediate  variable. 

(c)  For  the  portion  of  the  Invlscid  wake  up  to  r  =  O.B,  the  transverse 
distribution  of  total  pressure  is  the  same  for  uniform  flow  and  source 
flow.  Thus  the  distribution  in  imiform  flow  (p.  )  is  known  from  step 
(b). 

(d)  The  static -pres sure  distribution  in  uniform  flow  Pg  is  given  by 

u 

Eq.  (4.l6)  provided  (pg  )  is  known.  The  latter  quantity  is  calculated 

e 
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from  (p„  )  with  the  use  of  Hall's  correction  formula  of  Eq.  (4.17)* 
t-  s 
e 

(e)  The  required  Pitot-pressure  distribution  in  uniform  flovj,  p'  (r), 

u 

is  then  deduced  from  [p.  (r)]  from  step  (c)  and  p„  (r)  from 

2  u 

step  (d),  again  using  the  normal-shock  relations  with  the  Mach  number 
as  an  intermediate  flow  parameter. 

(f)  The  correction  for  the  experimentally  deduced  lateral  position  in 

the  inviscid  wake  is  simply  a  scaling  down  of  the  coordinate  r  with 

s 

the  ratio  R,  /R  ,  that  is,  r  =  (R  /R  )r  .  For  the  ratio  R  /R  we 
Vs'  'u'Vs  s  u's 

use  HaLl's  source-flow  solution  in  Eq.  (4.l8)  for  the  shock-v^ve  shape. 

To  ascertain  which  portion  of  the  Inviscid  profile  is  to  be  sub¬ 
jected  to  the  foregoing  correction  (that  is,  which  portion  satisfies  the 
condition  r  <  0.8),  we  need  the  value  of  either  R  or  R  .  For 
cases  where  the  bow  shock  wave  Intersects  the  wall  boundary  layer  down¬ 
stream  of  the  probing  station,  it  is  possible  to  determine  R  from  the 
Pitot -pres  siu:-e  trace  on  the  oscillograph  record.  In  this  situation  the 
transducer,  in  crossing  the  shock  wave,  experiences  a  marked  pressure 
Jimp  which  is  clearly  reflected  in  the  galvanometer  output.  For  stations 
where  the  intersection  occurs  ahead  of  the  station,  the  interaction  between 
shock  wave  and  boundary  layer  renders  it  inposslble  to  determine  R 

s 

experimentally.  We  then  use  the  following  empirical  formula  for  R^ 
from  Van  Rise  (I96C); 

R^d  0.98{x)°*^  .  (5.5) 

Viscous  Wake; 

(a)  The  static  pressure  within  the  viscous  wake  at  a  given  longitudinal 

station  is  uniform  and  equal  to  Pg  at  that  station.  It  has  been  shown 

e 

above  how  the  value  of  (p_  )  is  deduced  from  experimental  pressure 

e  ® 

data.  The  corresponding  value  of  (p^  for  uniform  flow  is  determined 
as  in  step  (d)  for  the  inviscid  wake. 

(b)  The  distribution  p_  (r„)  =  (p„  )„  is  used  with  the  measured 

£  S  c  s 

s  e 
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Pitot-pressure  distribution  p*  (r  )  to  calculate  the  variation  of  Mach 

°s 

number,  ),  according  to  Eq,.  (%l). 

s 

(c)  In  the  viscous  vake,  the  normalized  Mach-number  distribution 

Mg(r)/M2  is  assumed  to  be  the  same  for  lu.iform  flow  and  source  flow, 
e 

Since  the  distribution  for  .  oi  'ce  flow  is  known  from  step  (b)  and  since 
the  value  of  (M^  )  has  already  been  calculated  in  step  (e)  for  the 

P 

inviscid  wake,  the  distribution  (r^)  is  fully  determined. 

u 

(d)  The  Mach  number  (r  )  is  then  used  with  the  static  pressure 

u 

^2  ^^u^  ~  ^^2  ^u  step  (a)  to  calculate  the  required  Pitot-pressure 

distribution  p’  (r  )  in  uniform  flow, 
u 

(e)  The  lateral  position  coordinate  r  in  the  viscous  wake  is  corrected 

s 

according  to  Eq.  (4.30)  to  give  r^  in  uniform  flow. 

5.5.4  Results 

For  a  comparison  of  the  data  before  and  after  the  successive 
corrections  for  response  lag  and  source-flow  effects,  a  set  of  results 
from  a  typical  run  is  given  in  Fig.  51.  Shown  in  this  plot  are  (l)  the 
uncorrected  data  with  appropriate  normalization,  (2)  the  results  corrected 
for  response  lag,  and  (5)  the  final  results  after  the  profile  of  (2)  has 
been  properly  centered  and  then  subjected  to  source-flow  corrections. 

The  difference  between  adjacent  curves  is  considerable,  especially  between 
curves  (l)  and  (2),  whei’e  the  shift  in  time  (or  equivalently  the  lateral 
coordinate)  leads  to  large  magnitude  corrections  at  a  given  position 
within  the  profile.  For  the  particular  run  illustrated  in  Fig.  51,  the 
source-flow  correction  amounts  to  about  which  is  about  ten  times  the 
magnitude  of  the  source-flow  parameter  ~  O.OO'jQ.  For  certain  of 

the  more  distant  stations,  the  source-flow  correction  is  as  much  as  50^. 

Representative  data  from  selected  runs,  after  appropriate  correction, 
are  presented  in  Figs.  52  through  55*  The  dashed  curves  are  the  results 
corrected  for  response  lag,  and  includes  a  subsequent  centering  of  the 
profile.  The  application  of  the  source-flow  correction  to  the  dashed 
curves  then  leads  to  the  final  results  given  by  the  solid  curves.  The 
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results  of  Figs.  52,  55,  and  54  correspond  to  runs  made  during  the  final 
phase  of  the  testing  program,  in  which  only  the  near  stations  (stations 
A,  B;  C,  and  d)  were  probed.  The  results  of  Fig.  55  pertain  to  an 
earlier  series  of  runs  that  covered  the  entire  range  of  stations  con¬ 
sidered  in  Section  2.7,  the  farthest  position  being  at  x/d  =  15.6 
(station  H). 

The  later  run3  were  confined  to  the  lear  stations  for  two  reasons. 
First,  because  the  Pitot  pressure  decreases  with  distance  from  the  model, 
there  is  a  corresponding  increase  in  the  relative  level  of  mechanical 
noise  that  is  apparent  to  the  galvanometer.  The  recorded  date  for  the 
distant  stations  must  therefore  be  subjected  to  increased  curve  smoothing, 
and  in  the  process  some  accuracy  is  unavoidably  sacrificed.  The  second 
and  more  important  reason  for  favoring  the  near  stations  is  that  earlier 
experience  showed  that  the  inviscid  stream  in  the  tunnel  is  too  small  to 
permit  satisfactory  testing  at  the  more  distant  stations.  When  the  sphere 
is  placed  far  upstream  into  the  nozzle,  the  Interaction  between  the  bow 
shock  wave  and  the  tunnel  boxmdary  layer  appears  to  have  a  sizable  effect 
on  the  wake  at  the  probing  station.  The  profi?.e  in  this  instance  exhibits 
a  trough-like  shape  that  is  suggestive  of  turbulent  wakes  (refer  to  Fig. 

55  for  profiles  of  the  distant  stations).  Data  from  ^ther  sources,  such 
as  Vas,  Murman,  and  Bogdonoff  (1965)  and  McCarthy  (1962),  show  that  the 
wake  is  not  predominantly  turbulent  at  these  stations,  however.  We 
surmise  that  the  tunnel  boundary  layer  has  "infected"  a  considerable 
portion  of  the  inviscid  wake,  if  not  part  of  the  viscous  wake  as  well. 

For  this  reason  some  uncertainty  exists  in  the  data  for  stations  E  and 
beyond. 

With  regard  to  the  source-flow  correction,  we  have  used  the  scher.es 

derived  in  Sections  4.5  and  4.6  for  the  distant  stations  as  well,  even 

though  these  schemes  are  basically  for  the  situation  in  which  turbulent 

diffusion  is  uninportant.  This  has  been  done  for  want  of  a  better  method. 

* 

The  results  presented  are  grouped  according  to  Reynolds  number. 


The  Reynolds  niimber  is  based  on  the  sphere  diameter  and  on  calculated 
free-stream  conditions  at  the  nose  of  the  model  at  20  milliseconds 
after  the  start  of  the  run. 
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The  three  groups  of  runs  in  Fig.  52,  55,  and  54  correspond  to  the  three 
basic  Reynolds  numbers  at  vhich  the  tests  vere  performed  ( see  Section 
2.7).  The  nominal  test-section  Reynolds  number  for  the  runs  of  Fig.  55 
is  approximately  20,000  and  is  essentially  the  same  as  the  lov- Reynolds - 
number  series  of  Fig.  52.  Nevertheless,  tvo  separate  plots  have  been 
used  because  the  runs  in  Fig.  55  (runs  459  through  47l)  vere  madt  prior 
to  certain  modifications  to  the  tunnel  arc-chamber  and  collector  assembly. 
The  changes  in  operating  conditions  and  flow-decay  characteristics  that 
accompanied  the  modifications  necessitated  recalibration  of  the  tuniiel 
flow.  It  was  thought  desirable,  therefore,  to  consider  runs  459  through 
471  as  a  separate  series.  Since  the  flow  in  the  tunnel  after  modification 
decayed  at  a  slower  rate,  it  seems  reasonable  to  expect  that  the  data 
from  the  post-modification  runs  are  of  better  quality. 

Although  the  usable  traverse  of  the  probe  is  centered  roughly  about 
the  wake  axis  (approxiiiately  four  inches  to  either  side),  we  have  pre¬ 
sented  only  the  data  from  the  upper  half  of  the  profiles.  If  the  whole 
profile  is  considered,  some  asymmetry  is  noticeable  between  the  upper 
and  lower  halves,  even  after  appropriate  normalization  and  response 
correction.  We  presume  that  the  asymmetry  is  not  due  to  system  response 
lag,  since  the  correction  scheme  has  been  shown  in  the  exan5)le  of  Section 
5.6  to  be  capable  of  removing  the  associated  shape  distortion.  The  most 
plausible  explanation  for  the  asicnmetry  is  that  the  effects  of  flow  decay 
are  not  con5)letely  accounted  for  by  the  process  of  nojTnalization.  Given 
this  circumstance,  it  seems  reasonable  to  eiiq)hasize  one-half  of  the  pro¬ 
file.  The  upper  half  was  chosen  because  the  information  on  this  side  is 
more  complete.  This  choice  also  minimizes  interference  of  the  probed 
region  from  the  rest  of  the  traversing  probe,  should  such  interference 
exist. 

5.4  Concluding  Remarks 

In  Sections  5*2.2  and  5«5.4  we  have  presented  the  Pitot-pressure 
results  obtained  with  the  two  traversing  probes.  No  extensive  discussion 
was  given  to  interpret  the  results  because  of  lack  of  experimental  know- 
leuge  of  other  quantities  such  as  velocity  and  total  tenperature.  Since 


79 


verall  flow  picture  in  the  wake  behind  a  sphere  is  relatively  com¬ 
plex  is  not  possible  to  make  stn5)lifylng  assumptions  and  deduce  the 
variation  of  other  flow  q,uantitiea  from  the  meastired  Pitot  pressure  alone. 
A  case  in  point  concerns  the  assungition  f often  made  in  hypersonic  inviscid 
flc'.')  that  the  velocity  is  approximately  constant  and  equal  to  the  fr>^e- 
stream  val\ie.  This  assungjtion  is  invalid  in  the  viscous  wake,  which 
unfort 'onately  is  also  the  main  region  of  interest  in  wake  studies.  In 
the  viscous  wake,  it  is  thus  not  possible  to  infer  density  variations 
from  Pitot-pressure  results. 

A  considerable  portion  of  the  research  effort  in  this  study  was 
directed  toward  developing  the  rapidly  traversing  probes  and  toward 
perfecting  the  technique  of  testing  behind  a  support-free  sphere.  From 
the  standpoint  of  mechanical  operation,  these  techniques  have  been  shown 
to  be  feasible.  The  success  of  the  overall  scheme  of  mapping  Pitot 
pressure  is  licaipered  somewhat,  however,  by  inadequate  instrumentation 
response  and  by  source-flow  effects  in  the  tunnel  stream.  Both  these 
factors  can  lead  to  substantial  errors  in  the  recorded  data.  It  is 
possible,  however,  to  correct  for  these  effects  with  the  use  of  approxi¬ 
mate  schemes  devised  herein.  When  these  corrections  are  taken  into 
account,  the  '.;<perimei'tHl  systems  oind  techniques  -:volved  in  the  study 
appear  to  me.  .’t  consideration  as  a  means  of  'tudying  the  near  wake  of  a 
cold,  support-free  sphere. 
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APPENDIX  A 


EXPERIMENTAL  DETERMINATION  OF  PRESSURE-mSTRUMENTATION 
TRANSFER  FUNCTION 

Ip  Section  3.3,  a  discussion  vas  gi"en  of  the  transfer  function 
G(s)  of  the  transducer-galvanometer  system.  The  method  cf  deducing  the 
result  of  Eq.  (3.3)  is  now  described.  The  experimental  procedure  is 
represented  schematically  in  Fig.  36.  A  pressure  step  of  suitable  size, 
generated  by  a  sinple  shock  tube,  is  impc  ed  on  the  pressure  transducer. 
The  response  of  the  transducer  as  seen  at  the  amplifier  output  i^  recorded 
on  an  oscilloscope,  whereas  the  galvanometer  output  is  recorded  on  the 
oscillograph  that  is  a  normal  part  of  the  pressure-recording  system.  It 
is  thus  possible  to  consider  the  response  of  the  overall  system  in  two 
stages,  a  transducer-anplifier  combination  and  the  galvanometer.  Separate 
transfer  functions  are  fouind  for  these  two  sub-systems.  The  overall 
transfer  function  is  then  the  product  of  the  two  component  functions. 

The  shock-tube  arrangement  for  generating  the  pressure-step  input 
to  the  transducer  is  shown  in  Fig.  57*  The  transducer  is  mounted  onto 
the  end  plate,  the  inlet  tube  of  the  transducer  being  flush  with  the 
inside  of  the  plate.  As  many  as  three  transducers  of  the  wafer  type  can 
be  accommodated  in  the  same  test.  When  the  shock  wave  is  reflected  from 
the  end  plate,  at  time  t  =  t^,  say,  the  transducer  experiences  a 
pressure  junp  of  magnitude  p^  -  p^^.  The  pressure  level  then  remains  at 
p^  until  the  disturbances  that  result  from  the  interaction  of  the  con¬ 
tact  surface  with  the  reflected  shock  wave  reach  the  transducer,  say  at 
t  =  tg.  The  time  interval  of  interest  is  thus  t^  -  t^  at  the  input, 
since  it  is  during  this  interval  that  the  pressure  differential  imposed 
on  the  tran'^ducer  conforms  to  a  step  input. 

It  is  possible  to  relate  P^/P^  example, 

Liepmann  and  Roshko,  IS  p.  388).  For  the  present  experiments,  two 
sets  of  conditions  were  used;  these  are  as  follows: 
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P4/P1 

P5/P1 

M 

sw 

9-3.7 

6.60 

1.581 

24.30 

12.50 

1.890 

(The  test  gas  in  there  experiments  vas  air;  the  initial  teii^jeratures 

and  Tj^  were  assun'jd  to  be  equal.  )  A  pressure  Jurq)  p^  -  p^  of  any 
desired  size  could  then  be  obtained  by  choosing  appropriate  values  of 
p^  '  '.d  pj^.  Two  sets  of  conditions  are  useful  because  for  certain 
values  of  p^  -  p^,  the  pressure  level  of  p^  or  Pj^  calculated  from 
one  set  may  be  difficiilt  to  control  accui^tely,  in  which  case  the  other 
set  of  conditions  is  usually  found  to  give  more  sat^sfactory  results. 

The  weeiker  shock  wave  does^  however,  yield  a  longer  interval  -  t^ 
(approximately  4  milliseconds  conqjared  to  about  2.5  milliseconds  for  the 
case  of  the  stronger  shock  wave)  and  is  preferred  where  a  choice  exists. 

Typical  response  curves  are  shown  in  Fig.  58.  The  response  curve 
of  Fig.  38(a)  pertains  to  the  first  stage,  the  transducer-auiplif ier 
combination,  while  the  curve  of  Fig.  58(b)  relates  to  the  entire  system 
including  the  galvanometer.  In  the  oscilloscope  record  of  Fig.  58(a), 
the  upper  trace  is  the  signal  from  a  pressure  transducer  located  about 
one  and  a  half  feet  upstream  of  the  end  plate.  This  signal  serves  to 
trigger  the  oscilloscope  before  the  arrival  of  the  shock  wave  at  the 
transducer  under  test.  It  is  thus  possible  to  record  the  response  at 
time  t  =  t^. 


The  response  curve  at  the  amplifier  output  may  be  idealized,  as  in 
Fig.  39(a),  to  a  step  of  height  k  and  with  a  frontal  plane  that  is 
inclined  rather  than  vertical.  The  quantity  6  is  the  time  taken  for 
the  output  to  reach  the  value  k  and  in  the  present  instance  is  taken 
to  be  the  rise  time  of  the  transducer-amplifier  stage.  The  output  level 
k  is  a  normalized  quantity  and  its  value  is  given  by  the  magnitude  of 
the  output,  based  on  static  calibration,  normeilized  with  respect  to  the 


The  pressure  gage  and  manometers  used  to  measure  p^^  amd  pj^  were 
such  that  reading  errors  (percent)  were  more  serious  for  certain 
values  of  the  pressure  them  for  others. 
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input  pressure  p^  -  p^.  (The  latter  quantity  is  thus  shown  as  1  in 
Fig.  39(a)-)  In  the  sense  that  k  relates  the  dynamic  behavior  of  the 
system  to  its  static  calibration  factor,  it  is  appropriate  to  define  k 
as  the  "dynamic  calibration  factor".  With  the  foregoing  definitions  of 
k  and  6,  it  can  be  seen  that  if  the  system  is  linear,  the  shape  of 
the  normalized  response  curve  will  be  indtpendeut  of  the  magnitude  of 
thi.  input  pressure  p,.  -  p^. 

The  overall  system  was  first  checked  for  linearity  by  ingjosing 
pressure  steps  that  ranged  from  0.1  to  1.5  psi  differential.  For  this 
range  of  pressure  differentials,  which  corresponded  to  the  levels  en¬ 
countered  in  the  wake  tests,  there  was  no  detectable  dependence  of  the 
normalized  response  on  the  magnitude  of  the  input  pressure.  According 
to  the  response  studies  of  Smotherman  and  Maddox  (1963),  the  rise  time 
&  does  depend  on  the  magnitude  of  the  input  pressure,  which  implies  that 
the  transducer  is  essentially  nonlinear.  However,  the  variation  in  6 
is  significant  only  when  relatively  low  pressures  are  involved,  as  encoun¬ 
tered  for  example  in  measuring  static  pressure.  For  the  present  appli¬ 
cation  to  Pitot  pressures,  the  transducer  was  assumed  to  have  linear 
response  and  a  correspondingly  constant  rise  time.  Aside  from  the 
possible  small  discrepancies  in  the  input  independence  of  6,  the 
linearity  of  the  overall  pressure  system  was  confirmed  by  experiment. 

The  pressure  step  was  then  held  fixed  at  1.0  psi  f c  convenience  and  a 
sequence  of  ten  tests  were  performed. 

An  unforeseen  result  of  these  tests  was  the  discovery  that  initially 
the  pressure  level  (p^-p^^)  measured  by  the  transducer,  based  on  the 
static  calibration  factor,  did  not  correspond  to  the  theoretically 
expected  value  of  the  pressure  jumped  generated  by  the  shock-wave 
reflection.  That  is,  the  dynamic  calibration  factor  k  wes  found  to  be 
different  from  unity.  After  a  check  of  the  relations  uceu  for 
and  of  the  associated  calculations,  the  possibility  of 

erroneous  computation  was  ruled  out.  Also,  several  transducers  used  in 
the  same  test  consistently  measured  different  values  for  (P^'Pi)^^  with 
typicil  values  of  k  ranging  from  0.6  to  13.  This  suggested  that  the 
discrepancy  was  in  the  pressure-recording  system.  Subsequent  checks  for 
possible  anqjlifier  or  galvanometer  malfunction  showed  that  none  of  the 


system  cotiQ>onents  external  to  the  transducer  were  at  fault.  Furthermore, 
it  was  established  from  the  oscillograph  record  that  pressure  traces 
tended  to  settle  to  the  correct  values  eventua." ly,  indicating  that  the 
discrepancy  was  associated  with  the  dynamic  behavior  of  t;  ^  prr  ire 
transducer.  This  anomalous  response  characteristic  is  attributed  to 
faulty  fabrication  of  the  transducers,  most  probably  in  the  construction 
and  positioning  of  the  sensing  coils  and  the  spacing  of  the  pole  pieces. 
The  reason  for  this  hypothesis  is  that  there  appears  to  be  a  definite 
correlation  between  the  performance  of  a  parti -alar  transducer  and  how 
well  it  balances  in  the  an^ilifier  circuit. 

Because  of  the  unpredictable  nature  of  the  response,  the  transducers 
must  be  calibrated  dynamically  on  an  individual  basis  if  they  are  to  be 
used  for  measurement.  For  the  transducer  actually  used  in  the  vertically 
traversing  probe,  the  dynamic  calibration  factor  k  was  found  to  be 
1.04.  An  average  value  for  the  rise  time  B  was  found  to  be  0.211  milli¬ 
seconds. 

The  transfer  function  for  the  transducer-arnplifier  stage  can  now  be 
calculated  directly  by  taking  the  Laplace  transform  of  the  idealized 
response  curve  of  Fig.  39(a)‘  The  evaluation  of  the  Laplace  transform 
is  sin5)lified  if  the  response  diagram  is  taken  to  be  the  difference 
f '  -  f",  where  f  and  f"  are  shown  in  Fig.  59(b)  to  be  two  parallel 
lines  of  slope  k/6  whose  t-intercepts  are  0  and  6  respectively. 

The  Laplace  "transform  of  f '  and  f"  are  standard  expressions.  The 
transfer  function  of  the  transducer-amplifier  stage,  Gj^(s),  is  thus 
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_  (k/6)(l/s^)-(k/8)(e'^^/s^)  _  k  ,  -sBx 
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The  transfer  function  for  the  galvanometer  stage  was  taken  to  be 
the  standard  function  for  a  second-order  system,  which  is  clriaracteristic 
of  galvanometers  in  general.  The  response  curve  of  Fig.  38(b)  provides 
the  necessary  check.  It  is  possible  to  deduce  from  this  curve  the 


natural  undamped  frequency  o)^  and  the  damping  factor  ^  f«.r  the 
galvanometer  in  question.  These  values  correspond  closely  to  the  manu¬ 
facturer's  specifications.  The  response  curve  of  7ig.  36(b)  represents 
the  overall  system  response,  which  includes  the  transducer-amplifier 
stage.  Since  the  rise  time  6  is  small  compared  with  the  overall  system 
rise  time,  the  matching  with  a  second-order  system  would  be  valid.  The 
standard  second-order  transfer  function  is  given  by 
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The  overall  system  transfer  function  is  thus 


G(s)  =  Gj^Cs)  GgCs) 
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APPE2JDIX  B 


THE  Z-TRANSFORM  OF  THE  TRANSFER  FUNCTION 

The  system  transfer  fiinction  for  continuous -data  signals  has  been 
found  in  Appendix  A  as 
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In  Section  3.4,  a  fictitious  data-sampling  process  was  introduced  into 
the  transducer-galvanometer  system  to  describe  tl*e  data-reduction  procedure 
of  reading  discrete  data  points  off  the  oscillograph  record.  In  the 
mathematical  model  for  the  resulting  sampled-data  system  (discussed  in 
Section  3.5)>  a  fictitious  data  hold  was  inserted  to  interpolate  the 

data  between  adjacent  ?Fjnpling  instants  t  =  0,  T,  2T,  5T, .  For 

the  transducer-gal-vanometer  system  with  the  transfer  function  of  Eq. 

(B.l),  a  first-order  data  hold  Q(z,s)  of  Eq.  (3.10)  was  found  to  give 
the  best  results.  The  transfer  fxmctlon  Including  the  data  hold  is  then 


(5(s,z)  =  G(s)  Q(z,s) 

(B.2) 

six/) 
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From  the  theory  of  sampled-data  systems,  the  z-transform  of  S(z, s)  is 
given  by  (see  Ragazzlnl  and  Franklin,  1958^  P*  57) 


2(z)  =  ^  Residue 

poles  of 

G(z,p) 
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where  p  is  a  dummy  variable  used  in  place  of  the  Laplace  variable  s. 
The  function  ?r(z,p)  of  Eq.  (B.2)  has  a  third-order  pole  at 
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p  =  0,  and  simple  poles  at  p  =  -o.,  -a,  where  a  is  the  complex 
conjugate  of  a,  and  a  =  ^  +  i  Vl-C^)^  since  ^  <  !•  Equation 

(B.3)  mai’’  now  he  written 
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?t(z)  -  Z  .  fles  z(l-z‘^f 


p=0, -a, -a 
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p=o,^a,-a  I  jPi  p+a)( p^)(  l-e^^z"^] 


the  last  step  being  permissible  since  the  factor  that  is  not  considered 
in  the  evaluation  of  the  residues  does  not  contain  the  variable  p 
explicitly. 

Residue  at  p  =  0;  Define  ^Cp)  as 
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Then  the  residue  at  the  third -order  pole  at  p  =  0  is 


Res  (O)  -  i  ^-^1 
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where  c^  =  -2^6  -6/2  , 

c^  =  (2^-»T)6  +  6^/2  , 


and  6  =  cu  6,  T  =  T 
n  '  n 


R-sldues  at — £_  f  -  -Q«‘  It  is  convenient  to  evaluate  the  residues 

at  p  =  -a  and  p  =  -a  simultaneously,  because  of  the  con,1ugate  nature 
of  the  roots  a  and  a.  We  then  have 

Res(-a)  +  Res(-a)  =  _ _ ^  i  ^  1  -  . 

lp=-a  7Fp^<r)(i-eP^z-^)!p=.a 
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vhere  dj_  =:  2e"^  cos  XT  , 

^  , 

fQ  =  -Sin  30  -  e^^  sin  (X8-"')  , 

f^  s  -e  sin  (XT-50)  +  sin  (XT+X5-30)  , 

and  X  =  Vi-5^'  ^ 

0  =  tan  (X/5),  from  the  notation  a  =  co  e^^  . 

n 

Combining  Eqs.  (B.4),  (B.6),  and  (b.?).  ve  obtain 

(5(z)  >  ,_k  °-i^  °o  •*• 

XTB  1  +  dj^z  ^  +  dgZ~^ 


(B.7) 


(B.8) 


vhere  m^^^  =  Xc^  +  f^  , 

“O  “  ■*■  -  2f3  +  f^  , 

"i  =  ^0^2  +  ^‘^I'^i  +  ^0  ■  ^ 

“2  =  ^^1^2  ■*•  ^1  • 

The  system  parameters  vere  found  experimentally  to  have  the  values 
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k  =  1.04  , 

6  =  0.211  milliseconds  , 

i  =  0.6266  , 

=  1900  radians/sec. 

The  saii5)ling  period  T  vas  taken  to  be  O.25  milliseconds,  even  though 
during  data  reduction,  data  vere  read  off  the  oscillograph  record  every 
one-half  millisecond.  A  value  of  0.25  milliseconds  for  T  vas  found  to 
give  more  accurate  resvilts  coirpared  to  the  case  vith  T  =  O.5  milliseconds 
because  vith  the  sampling  rate  doubled,  the  first-order  or  linear  data 
hold  reconstructs  the  sanpled  data  into  a  curve  that  resembles  more 
closely  the  original  continuous  data  curve.  The  information  at  the 
sanpling  instants  betveen  data  readouts  vas  obtained  vith  the  use  of 
interpolation  formvilas  (see  discussion  at  the  end  of  Section  5.5)» 

With  the  values  of  the  system  parameters  thus  specified,  it  is 
possible  to  evaluate  the  constants  that  appear  in  Eq.  (B.8).  An  electronic 
computer  vas  used  to  perform  the  arithmetic  involved  in  transforming  the 
quotient  into  a  series  in  z"^: 

S(z)  =  r_^z  +  Tq  +  +  ^2^”^  ■*■  *"  *  (B*9) 

The  values  of  a  representative  number  of  r^'s  are  listed  in  the  table 
at  the  end  of  this  appendix.  It  can  be  seen  from  Eq.  (5.1J)  that  the 
n 

sirni  Y  represents  the  response  of  the  system  at  t  =  nT  to  a 

i^l  ^ 

unit-step  input.  Thus  the  alternation  in  sign  of  successive  sequences 
of  8  or  9  r^'s  indicates  that  ojcillations  (of  approximate  period  16t, 
say)  are  present  in  the  calculated  response  output.  The  general  trend 
of  decreasing  coefficients  beyond  r^  inplies  that  the  oscillation 
eventually  damps  out.  These  tvo  observations  are  in  keeping  vith  the 
characteristic  behavior  of  an  underdamped  system. 

The  fact  that  the  coefficient  r  ^  is  nonzero  means  that  the  out¬ 
put  has  a  deflection  at  t  =  -T,  vhlch  is  physically  inadmissible.  This 
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behavior  is  due  to  the  inclusion  of  the  first-order  data  hold  that 
effects  a  linear  interpolation  between  the  zero  input  at  t  =  -T  and 
the  unit-step  input  at  t  =  0.  For  a  similar  reason,  r^  is  not  zero  as 
it  should  be  in  principle.  In  the  application  of  the  coefficients  r^'s 
in  the  correction  scheme  (see  Appendix  c),  r  ^  was  taken  to  be  zero, 
whereas  the  nonzero  value  of  r^  was  retained.  This  approach  was  con¬ 
sidered  reasonable,  since  the  output  ceLLculated  in  this  manner,  for  a 
uiiit-step  input,  agreed  well  with  the  experimentally  recorded  output 
discussed  in  Appendix  A. 
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Coefficierts  in  Series  Expansion  of  ^(z) 


i 

^i 

i 

^i 

-1 

-6. 4792689, -03 

25 

1.4417189,-04 

0 

2.3199572,-02 

26 

-5.8938240,-06 

1 

9.9651717,-02 

27 

-8.7658587,-05 

2 

2.0891119,-01 

28 

-1.1811693,-04 

2.3429669,-01 

29 

-1.1520153,-04 

4 

2.0919678,-01 

30 

-9.4369901,-05 

5 

1.6044702,-01 

31 

-6.7135251,-05 

6 

1.0679144,-01 

32 

-4.0914622,-05 

7 

5.9384394,-02 

33 

-1.9628658,-05 

8 

2.3333841,-02 

34 

-4.6157902,-06 

o 

-4.3871823,-04 

35 

4.4327616,-06 

10 

-1.3474022,-02 

36 

8.6825478,-06 

11 

-1.8413414,-02 

37 

9.5770539,-06 

12 

-1.8064358,-02 

38 

8.4721437,-06 

13 

-1. 4857^^7,  -02 

39 

6.4491078,-06 

14 

-1.0609776,-02 

40 

4.2573953,-06 

15 

-6.4970768,-03 

Ji-l 

2.3384101,-06 

16 

-3.1450905,-03 

42 

8. 9003754,  .-07 

IT 

-7.7192826,-04 

43 

-5.7138223,-08 

18 

6.6547992,-04 

44 

-5.6988905,-07 

19 

1.3470357,-03 

45 

-7.5752768,-07 

20 

1.4980708,-03 

46 

-7.3458316,-07 

21 

1.3313658,-03 

47 

-5.9934902,-07 

22 

1.0172732,-03 

-4.2476382,-07 

23 

6.7432266,-04 

49 

-2.5761347,-07 

24 

3.7268871,-04 

50 

-1.2245572,-07 

APPENDIX  C 


DERIVATION  OF  RESPONSE  CORRECTION  SCHEME 


In  Section  3.5#  ve  arrived  at  expressions  for  che  z-transform 
Pou^(^)  "the  pressure  output  and  G(z)  for  the  system  treuisfer 

function.  These  expressions  etre  given  in  Eqs.  (5.7)  and  (3.12),,  respec¬ 
tively.  The  problem  in  Section  3.6  Is  to  deduce  the  input  z-transform 

/  \  * 

P^(z;  from  the  relation 

'■cut'’")  ' 

It  is  convenient  in  what  follo!*s  to  use  the  sin^pler  notation 


X(a)  -  Pj^(2)  , 


Equation  (C.l)  then  becomes 


with 


Y(z)  =  G(z)  •  X(z)  , 

-1  -2 

X(z)  =  Xq  +  x^z  +  x^z  -!•  ... 

•"  *2 

Y(z)  =  Xq  ^1^'^ 

~/  \  -1  -2 

G(  z )  =  r_^z  +  r^  +  r^z  +  r^z  +  • • • 


(C.2) 


(C.5) 


The  coefficients  x  ,  y  are  the  values  of  th’  input  eind  output, 
respectively,  at  t  =  ml  (see  Fig.  i*0).  The  term  r_^  in  ^z)  is 
a  spurious  effect  of  the  inclusion  of  the  first-order  data  hold,  since 
r  ^  is  an  order  of  meignitude  smaller  than  r^  (see  values  listed 


This  relation  is  identical  to  Eq^. 
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Appendix  B)^  the  deletion  of  the  z  term  in  ?i(z)  does  not  cause  any 
serious  error  but  does,  on  the  othei  hand,  simplify  the  presentation  of 
the  ensuing  analysis.  The  coefficient  r^  is  retained,  however,  because 
it  does  not  cause  any  added  difficulty  in  the  algebra. 

Substituting  Eqs.  (C.3)  into  Eq.  (C.2)  and  equating  the  coefficients 


of  z  ,  ve  obtain 


7  =■  /  r.x  .  . 


(C,4) 


Let  x',  x"  be  respectively  the  i^^  and  (i+l)^^  .pproximation  of  x. 
Then  from  Eq.  (C.4)  we  have 


y*  =  j  r.x'  .  , 

^  jfe  J"--' 


V**  s  \  V** 


(C.5) 


where  y'  and  y  are  the  outputs  corresponding  to  the  x'  and  x" 
sequences  respectively.  With  the  notation  e  =  x”  -  x',  Eqs.  (C.5) 
can  be  rewritter  as 


y"  -  y'  "  T'  r.e  .  . 

^m  ^m  J  m-J 


(C.6) 


Now  from  Eq.  (C.4)  and  Fig.  40,  it  can  be  seen  that  y  =  5^  r. 
is  the  response  of  the  system  at  t  =  nff  to  a  unit-step  input 
at  t  =  0.  For  a  system  with  a  dynamic  calibration  factor  k  =  1.04, 


we  have 


lim  ^  r.  =  k  =  1.C4  . 

m  -» 00  ^ 


(C.7) 


Consider  the  values  of  r  through  r.,^  from  the  table  of  Appendix  B: 
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0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


0.0232 

0.0997 

0.2089 

0.2545 

0.2092 

0.1604 

0.1068 

0.0594 

0.0255 

-  0.0004 

-  0.0X55 


From  this  ve  find  that  the  sum 


6 

5 


rj  has  a  value  of  1.0425.  Thus 


Tq  through 


r^  roughly  speaking  accoxjnt  for  the  correct  response  that 

corresponds  to  the  limiting  vaJue  given  hy  Eq..  (C.7).  The  remaining 

rj*s  represent  the  oscillations  that  occur  in  the  output  of  an  vuider- 

danqped  system.  If  the  errors  €  tiiroughout  Eq,.  (C.6)  are  assumed  to 

A 

be  of  the  same  order  of  magnitude,  the  net  contribution  of  ,  r.e  , 
vill  then  be  approximately  zero,  and  from  Eq.  (C.6)  we  can  write 


v”  -  v’  = 


/  r.e  , 
3%  •1”-^ 


(C.8) 


The  iiiq)ortant  assunq)tion  is  now  made  that  the  errors  e  ,  for 

m- 

J  =  C  through  6,  form  a  linear  sequence  ( see  sketch  on  next  page); 
that  is,  they  are  expressible  in  the  form 


c 


m-J 


-  ^0  -  (3)^ 


(C.9) 
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’m-6  ^111-5  ^111-4  ^ni-3  in-2  n-l  m 

^6  "■5  ""3  ""2  ^1  "‘o 


It  is  then  possible  to  use  the  idea  of  a  weighted  average  to  determine 
which  of  the  errors  most  accurately  satisfies  the  relation 


e  /  r.  =  r.e  .  . 

^  J  jro  J 


(C.IO) 


Letting  and  solving  for  v  from  Eqs.  (C.9)  and  (C.^C), 

ve  obtain 


6 

Z 


^  "■j 
j=o 


=  5.46 


(c.n) 


Thus  if  V  is  approximated  by  5#  becomes  Eq..  (C.8) 

may  be  approximated  as 


^m  ^m  ^m-5  J 


Since  e  =  x"  -  x',  this  equation  can  be  rearreuiged  to  give 


(C.12) 
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(C.13) 


x"  =  X 


•  + 

n 


C3 "  yA+3 


J=o 


vhere  (iii-3)  has  been  replaced  by  n.  On  the  right-hand  side  of  Eq. 
(C.l3)f  y”  is  the  output  corresponding  to  the  improved  input  x".  It 
is  then  reasonable  to  assume  that  if  the  recorded  output  y  is  used  in 
place  of  y”,  the  quantity  x"  calculated  from  Eq.  (C.l3)  is  closer 
to  the  actual  input  x  and  therefore  represents  an  improvement  on  x ' . 
We  can  thus  vrite  the  final  result  as 


(C.14) 


In  deriving  this  equation,  the  error  sequence  used  goes  back  to  (m-6), 
or  (n-3);  the  derivation  is  thus  applicable  only  for  n  >  3. 

For  n  -  2f  the  e  corresponding  to  r^  is  zero  (refer  to  sketch). 
The  summation  in  Eq.  (C.IO)  is  carried  to  J  =  3  only,  resulting  in 
V  =  3.06  ^  3.  Hence  the  only  modification  required  in  Eq.  (C.l4)  is  in 
the  summation  of  r^'s  as  follows; 


3 

^  r 


(C.I3) 


In  a  similar  manner,  expressions  can  be  derived  for  n  =  1  and  n  =  0 
as  follows; 
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+ 


(C.16) 


The  initial  guess  for  the  input  is  taken  to  be  the  recorded  pressure 
output  shifted  backvjard  in  time  by  1  millisecond.  The  correction  scheme 
vorks  well  with  this  guessed  input  except  where  jump  discontinuities 
occur.  Consider  for  example  the  input  shape  of  Fig.  15  and  the  corres¬ 
ponding  output.  In  the  neighborhood  of  the  point  t  =  5  milliseconds, 
the  shifted  output  differs  significantly  from  the  input,  and  the  cor¬ 
rection  scheme  may  break  down  here.  To  remedy  the  situation,  the  portion 
of  the  shifted  output  that  causes  the  difficulty  is  deleted  and  replaced 
by  a  linearly  extrapolated  segment,  shown  by  the  dashed  line  in  Fig.  15* 

A  junp  discontinuity  is  provided  for,  the  point  at  which  it  is  placed 
in  the  guessed  output  being  detemiined  fairly  accurately  from  where  the 
sharp  rise  occurs  in  the  output.  The  guessed  input  for  the  segment  to 
the  left  of  t  =  5  milliseconds  may  be  any  straight  line  with  approxi¬ 
mately  the  same  slope  as  that  of  the  recorded  output  in  that  interval. 

The  indirect  correction  scheme  derived  here  is  not  unique.  It  is 
conceivable  that  two  slightly  different  inputs  will  give  the  same  output. 
For  instance,  the  same  signal  superiii5>osed  on  two  carrier  voltages  of 
different,  but  high,  frequencies  will  probably  produce  such  an  effect. 

It  is  then  not  possible  to  differentiate  between  the  two  cases  when  the 
inputs  are  to  be  calculated  from  the  identical  output.  A  case  in  point 
is  when  the  recorded  output  contains  small -amplitude  oscillations  that 
are  superimposed  on  the  main  Pitot-pressure  curv^e.  These  are  in  actuality 
due  to  random  variations  in  the  tunnel  flow  or  to  ramnents  of  the  mechan¬ 
ical  noise  that  are  not  completely  removed  by  the  monitoring  accelerometer 
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output.  When  the  recorded  output  Is  subjected  to  the  correction  scheme 
however^  the  oscillations  are  taken  to  be  the  result  of  a  filtering 
process  of  large -an^jlitude  oscillations  in  the  original  pressure  trace. 
Successive  iiqjrovement  on  the  guessed  input  then  introduces  an  unstable 
iterative  procedvire.  This  difficulty  was  in  fact  encountered  in  the 
data-reduction  procedure.  To  eliminate  the  spurious  effect^  a  certain 
amount  of  curve  smoothing  of  the  recorded  output  was  found  to  be  necessary 
before  it  could  be  used  as  the  initial  guessed  input.  This  problem  is 
discussed  in  some  detail  in  Section  5*3*2.  Provided  that  this  Instability 
condition  is  taken  into  consideration^  the  indirect  correction  scheme 
derived  herein  is  Judged  to  be  adequate  for  the  present  application. 
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APPENDIX  D 


SOURCE-FLOW  EFFECTS  ON  THE  BOUNDARY  LAYER 


The  boundary-layer  transformations  for  two-dimensional  flow  used  by 
Cohen  and  Reshotko  (1956)  are  based  on  Stewartson's  transformation  and 
can  be  written 


> 


V  = 


dn 


where  s  and  n  are  rondimensionalized  distances  tangential  and  normal 
to  the  body  surface;  &,  p,  and  p  are  respectively  the  speed  of  sound, 
pressure,  and  density  of  the  fluid,  and  is  a  constant  of  propor¬ 
tionality  in  the  linear  /iscosity  law  (n/ia.  )  =  X  (T/T.  ).  The  sub- 

ex  ex 

script  "t"  refers  to  total  conditions,  while  "ex"  pertains  to  the 
local  inviscid  flow  outside  the  boundary  layer  (external  flow).  The 
value  of  X^  is  determined  by  the  condition  that  the  viscosity  at  the 
model  temperature  given  by  the  linear  law  should  match  the  value  from 
Sutherland's  law.  Thus  X^  assuii.es  the  valu'' 


X  = 

V 


where  T  is  the  Sutherland  constant,  which  for  nitrogen  is  102. 

SU  7  0  1 

(see  Chapman  and  Cowling,  1952). 

For  axially  symmetrical  flow,  as  for  exan^ile  over  the  hemispherical 
nose  under  study,  the  transformation  must  be  modified  as  follows  (see 
for  exan^ile  Moore,  .1964), 
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I  = 


} 


(D.la) 


^  = 


(D.lb) 


vhere  r  is  the  nondlmenslonallzed  distance  measured  from  the  axis  of 
symmetry. 

The  Flmllarlty  variable  used  in  the  boundary-layp.r  solutions  is 
defined  by 


1 


(D.2) 


vhere  v.  is  the  kinematic  viscosity  at  T.  .  and  U  is  the 

u  6X 

ex  ex 

transformed  velocity  in  the  external  invJscid  stream,  related  to  the 

physical  velocity  by: 


U 


ex 


a 


ex 


ex 


(D.3) 


The  quantity  m  is  the  exponent  in  the  power-lav  description  of  the 
inviscld  flow 


,  (D.4) 

Where  C  is  a  constant.  In  the  present  application,  the  inviscld  flcp-j 
cannot  ue  described  by  this  equation  with  a  single  value  of  m  that  is 
applicable  over  the  entire  hemispherical  nose.  It  is  assumed  neverthe¬ 
less  that  the  power-lav  description  applies  locally,  so  that  at  a  parti¬ 
cular  point  of  Interest,  the  appropriate  value  of  m  to  take  is  the  one 
for  which  Eq.  (D. 4)  describes  the  flow  in  the  immediate  neighborhood  of 
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r 


1 


I 


r 


the  point. 

An  ingjortant  parameter  in  the  solutions  is  the  pres sure -gradient 
parameter  p,  vhich  is  defined  in  terms  of  m  as 


P  = 


2m 

m+1 


(D.5) 


A  second  parameter  is  the  enthaipy  function  at  the  model  surlace  given 

by  -  !•  For  the  present  experiments,  is  assumed  to  be 

^ex 

negligible  in  con^jarison  vJth  h.  ,  so  that  S  takes  on  the  value 

u  V 

ex 

-1.0. 

It  is  now  possible  to  show  that  the  values  of  P  for  uniform  flow 
and  source  flow  do  not  differ  significantly  from  each  other.  To  calculate 
P,  or  equi-valently  m,  we  see  from  Eq.  (D.4)  that  m  is  given  by  the 
expression  m  =  ^  which  can  be  evaluated  graphically 

from,  a  plot  of  versus  |.  This  plot  in  turn  can  be  deduced  from 

the  surface -pressure  distributions  in  Fig.  25  by  application  of  the 
transformations  (D.la)  and  (D.3),  since  by  virtue  of  the  isentrppic  con¬ 
dition  in  the  inviscid  stream  the  pressure  ratio  con5)letely  specifies 
the  flow.  The  values  of  P  thus  calculated  for  uniform  flow  =  O) 

and  source  flow  =  0.02)  at  <p  =  n/2  are  1.6t4  and  I.685 

respectively.  The  difference  (P  -P  )/p  is  O.OO66,  which  is  sub- 
stantially  smaller  than  the  source-flow  parameter 

We  now  examine  the  solutions  of  Cohen  and  Reshotko  for  the  total- 
enthalpy  ratio  h. /h.  and  the  velocity  ratio  V/V  particularly 

V*  U  ©X 

ox 

with  regard  to  their  rariation  with  the  pres sure -gradient  parameter  p. 
These  solutions  are  given  in  Fig.  4l  for  non-negative  values  of  p,  that 
is,  for  cases  where  the  pressure  gradient  is  favorable.  For  values  of 
P  approaching  the  limiting  value  of  2,  both  ratios  are  insensitive  to 
changes  in  p.  This  is  particularly  true  of  the  total-enthalpy  ratio, 
and  Lees  (1956)  for  example  invoked  the  associated  invariance  of  total- 
enthalpy  gradient  at  the  wall  to  justify  the  omission  of  the  latter  in 
con5)uting  the  laminar  heat  transfer  on  cold  blunt-nosed  bodies.  In  the 
present  instance,  since  no  actual  solution  of  the  boundary-layer  problem 
is  being  atten^ted,  it  is  not  necessary  to  neglect  the  pressure  gradient 
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m 


altogether.  V/e  can  assume,  however,  that  the  distribution  of  the  total- 
enthalpy  ratio  as  a  function  of  the  similarity  variable  4  is  the  same 
for  uniform  flow  and  source  •^low.  In  view  of  the  small  difference  of 
0.0066  between  P,,  and  p  ,  particularly  when  these  are  close  to  the 
limiting  value  of  2,  we  shall  extend  the  assumption  to  include  the 
invariance  of  the  velocity  ratio  from  uniform  flov;  to  source  flow. 

With  the  assumption  of  the  identity  of  the  ratios  of  total  enthalpy 

and  velocity,  it  is  now  possible  to  compare  the  uniform-flow  and  source- 

flow  distributions  of  the  static-tenperature  ratio  t/T  .  Starting  with 

\/ht  using  the  energy  equation  h.  =  ^  V  +  7Rr/(/-l),  one  obtains 

^ex  ^ 

for  T/T^^  the  expression 


fp  T  ^  ex^  h 


ex 


ex  'V 


(D.6) 


ex 


ex 


On  the  right-hand  side  of  this  equation  the  ratios  h. /h.  and  V/V 

V  Xt  GX 

ex 

are  taken  to  be  the  same  for  uniform  flow  and  source  flow,  and  are 
ootalnable  from  the  solution  of  Cohen  and  Reshotko.  The  values  for 
M  and  M  at  cp  =  n/2,  deduced  from  the  pressure  results  of  Fig. 

“X  €X 

u  s 

25,  are  respectively  2.8l  and  2.93.  The  temperature  distributions 

T, /T„„  and  T  /t^  thus  calculated  from  Eq.  (D.6)  are  shown  in  Fig. 
u  s 

42.  ThJ  curves  are  very  close  to  each  other  and  are  identical  for  values 
of  Tj  larger  than  2.5*  For  smaller  values  of  t|,  the  difference  is  of 
the  same  order  as  the  source-flow  parameter  =  0.02).  Within  the 

accuracy  of  the  present  study,  this  difference  is  negligible  and  the 
temperature-ratio  distribution  can  be  taken  to  be  identical  in  the  trans¬ 
formed  Tj-varlable;  that  is. 


for  any  given  t| 


(D.7) 


We  are  now  in  a  position  to  consider  how  the  boundary  layer  in  the 
physical  plane  differs  from  uniform  flow  to  source  flow.  The  transformation 
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from  the  similarity  variable  t)  back  to  the  physical  coordinate  n  is 
accomplished  with  the  use  of  Eqs.  (D.lb)  and  (D.2),  with  the  result 


n 


^  Pt 

ex  ex 

a  p 
ex  ex 


ex 


m+1  U 


ex 


1/2 


/  -  T - 

^0  r  t 

ex 


(D.8) 


On  the  right-hand  side  of  this  equation,  all  quantities  except  the 
integral  have  been  used  before  in  the  determination  of  P  and  are  there¬ 
fore  known.  By  virtue  of  the  relation  (D.?)#  the  ratio  of  the  integral 
for  uniform  flow  and  source  flow  can  be  written 


for  any  given 


(D.9) 


Here  we  have  used  the  condition  that  at  cp  =  :t/2 


r 

s 


+  n 

s 


r,  +  n 
b  u 


r 

u 


since  r,  is  much  larger  than  the  difference  n  -  n  . 
b  “  s  u 

With  these  considerations,  one  can  now  calculate  the  ratio  n  'n 

from  Eq.  (D.8).  Because  the  ratio  of  the  temperature  integrals  in  Eq. 

(D.9)  is  independent  of  i],  the  ratio  n  /n  is  the  same  for  the  entire 

s  u 

boundary  layer.  Hence  it  is  necessary  to  consider  only  the  ratio  of  the 
"boundary-layer  thickness"  6,  with  the  result 


5 

^  =  1.095  ,  for  r^/L^  =0.02  . 

u 
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For  application  to  cared  where  r^/l.^  is  smaller'  than  0.02,  a  linear 
interpolation  ia  assumed,  leading  to  the  final  result 

S  j- 

-g-  =  1  +  1^.75  ^  at  cp  =  jt/2  . 

u  o 
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APPENDIX  E 

COMPAKTSON  OF  NORMALIZED  MACH-NUMBER  PROFILE  FOR  UNIFORM 
FLOW  AND  SOURCE  FLOW 

In  the  development  of  the  source-flow  correction  scheme  for  the 

VISCOUS  wake  (Section  4.6),  it  was  assuned  that  the  distribution  of  the 

Mach-number  ratio  M/M^  was  the  same  for  both  uniform  flow  and  source 

flow.  This  condition,  Eq.  (4.27),  was  used  in  two  instances:  (l)  in 

reducing  Eq.  (4.26)  to  the  final  result  of  Eq.  (4.29)  for  wake-growth 

correction,  and  (2)  in  relating  the  uniform  flow  to  the  source  flow  in 

* 

tne  correction  for  the  mgnitude  the  Pitot  pressure.  Although  we 
cannot  verify  directly  the  accuracy  of  the  assumed  conaition  (M/M  )  = 

(m/M^)^,  we  can  nevertheless  demonstrate  its  plausibility  the  basis 
of  the  results  of  some  calculations  pertaining  to  a  related  exan^ile. 

The  related  proolem  may  be  stated  as  follows:  We  assume  that  at 
some  initial  longitudinal  position,  close  to  the  wake  neck,  the  distri¬ 
bution  of  M/Mg  is  the  same  for  uniform  flo\J  and  soi’rce  flow,  that  is. 


f  M(r2]  _ 

I  \  ls,A  ■  \  Ju,A  ^ 


for  0  <  r  <  1  , 


(E.l) 


where  f  in  this  appendix  denotes  and  the  subscri^it  "A"  refers 

to  the  initial  station.  This  Mach-number  distribution  can  be  deduced 
?rom  information  that  is  considered  given  in  the  context  of  the  present 
study.  The  problem  then  consists  in  showing  that  the  similarity  in  the 
normalized  Mach-number  profile,  Eq.  (E.l),  is  preserved  In  the  flow 
downstream  of  the  initial  station. 

Some  of  the  details  of  the  foregoing  related  problem  are  now  c..  .- 
sidered.  For  the  initial  station,  we  use  station  A  wiuh  an  x/d  of 
3.8,  this  being  the  nearest  station  studied  experimentt-lly.  The  Mach- 
number  distribution  M  .(r)  is  deduced  from  the  measured  Pitot -pres sure 

S  fPi 

Specifically,  this  refers  to  step  (5)  on  p.  65. 
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distribution  (r)  and  the  static  pressure  p(r)=p(r)=p  , 

the  last  quantity  being  taken  from  the  method-of -characteristics  results 
* 

of  Fig.  23.  The  use  of  experimental  data  for  the  Pitot-pressure  profile 
ensures  that  the  calculated  Mach-number  file  has  the  proper  distri¬ 
bution,  even  though  the  magnitude  may  not  necessarily  be  correct.  Since 

M  is  by  definition  M  .(l),  the  lefthand  side  of  Eq.  (E.l)  is 

s,A 

now  known.  The  Mach-number  distribution  for  uniform  flow  then 


follows  from  Eq.  (E.l)  with  the  determination  of  M 


This  quantity 


is  deduced  from  the  static  and  total  pressures  at  the  point  r  =  1.  The 
distributions  M  .(r)  and  M  .(f)  thus  calculated  are  listed  in 
col^umns  (2)  and  (3),  respectively,  of  the  table  at  the  end  of  this 
appendix. 

We  now  proceed  downstream  from  the  initial  station  and  shcnv  that 

the  similarity  in  the  distributions  of  the  Mach-number  ratio  is  presear^ed 

for  the  range  of  longitudinal  stations  of  interest.  To  establish  this, 

it  is  sufficient  merely  to  show  that  the  similarity  holds  at  the  most 

distant  station  tested,  that  is,  at  station  H  with  an  x/d  of 

To  this  end,  we  deduce  independently  both  the  uniform-flow  and  source- 

flow  distribution  of  the  Mach-number  ratio  at  station  H.  The  dis'iribution 

[M(f)/M  ]  t,  for  the  source-flow  case  can  be  cai  collated  in  a  manner 
6  n 

sinLilar  to  the  corresponding  distribution  at  station  A,  again  using  the 
roeasvired  Pitot -pressvire  profile  to  ensure  that  the  Mach-number  distribution 
adopted  is  reeilistic.  As  for  deducing  the  uniform- flow  distribution,  we 
note  that  at  station  H,  the  only  quantity  known  thus  far  is  the  static 
pressure  from  the  method-of-ch£iracteristics  results  of  Fig.  25.  To 
obtain  the  Mach-number  information,  a  second  pressure  quantity  is  needed, 
and  for  this  we  choose  the  toteil  pressure.  The  distribution  p.  (f ) 

The  fact  that  the  weike  static  pressure  (p^  in  figure)  is  given 

e 

as  the  ratio  Pg  /p^  in  this  figure  does  not  introduce  any 
“^e  n 

difficulties.  Since  only  pressure  ratios  are  involved  in  the  calcu¬ 
lations,  the  quantity  p^  cancels  out  in  all  cases. 


V 


can  be  calculated  as  follows.  For  each  of  the  ises  of  uniform  flew 

and  source  flow  at  station  A  and  of  soxirce  flow  at  station  H,  the  Mach- 

number  distribution  and  static  pressure  are  known.  It  is  thus  possible 

to  calculate  the  respective  total -pressure  distributions,  that  is,  tie 

distributions  p  Cr),  p  (r),  and  p  (r).  With  the  use  of 
u,j  s,A  s,H 

Eq.  (4.25)  that  relates  these  three  quantities  and  p  (r),  this 

u,H 

latter  distribution  can  then  be  calculated.  The  required  Mach-number 

distribution  [M(r)/M  ]  „  follows  immediately  from  p.  (r)  and 

®  u,H 

^^e  * 

u,H 

The  distributions  M  „(r)  and  M  „(r)  calculated  in  the  fore- 

n  n 

going  are  listed  in  coliimns  (4)  and  (5)  of  the  table,  while  the  ratios 
are  given  in  columns  (6)  and  (?)•  We  see  that  whereas  the  Mach  r umber 
Itself  differs  considerably  for  uniform  flow  and  source  flow,  the  "atios 
differ  only  slightly  for  the  two  cases.  From  this  the  conclusion  can 
be  drawn  that  the  similarity  of  the  distributions  of  Ifech-number  ratio 
is  indeed  preserved  in  the  flow  downstream  of  the  initial  station. 
Strictly  speaking,  this  conclusion  applies  ily  to  the  related  problem 
that  uses  Eq.  (E.l)  as  the  starting  point.  The  result  does  not  provide 
direct  proof  for  the  assvraption  throughout  the  entire 

viscous  wake.  We  note,  however,  that  the  basic  assumption  of  the  related 
problem,  Eq.  (E.l),  is  itself  in  line  with  the  results  of  calculations 
pertaining  to  the  development  of  the  boundary  layer.  It  was  showi  in 
Appendix  D,  Fig.  4l,  that  the  temperature  ratio  T/T  did  not  diffei 
significantly  for  uniform  flow  and  source  flow.  This  iTr5)lies  that  the 
difference  in  the  Mach-number  ratio  M/M^  ^  is  even  smaller,  thereby 
justifying  the  assumption  of  Eq.  (E.l).  On  this  basis,  it  appears 
reasonable  to  use  the  result  of  the  related  problem  to  substemtiate  the 
assunptiC'.i  (M/Mg)^  =  (m/M^)^  of  Section  4.6. 
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Test  stations. 


(a)  Sphere  suspension. 


(b)  Side  view  of  system  for  accelerating  mass  (loaded  position). 
Figure  Sphere  suspension  and  release  mechanisms. 
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Figure  7.  Trigger  assembiy. 
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Figure  II.  Axially  traversing  probe  mounted  In  tunnel. 


(a)  System  for  continuous  signals. 


(b)  Idsol  input  somplsr  and  rsprsssntotion  by  unit- impulse  train. 


(c)  Sampied-doto  system. 

Figure  13.  The  transducer-galvanometer  system  as  a  sampled-data 
control  system. 
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Sampling  period  T  =  0.25  milliseconds 


Figure  15.  Comparison  of  reconstructed  input  with 
given  (constructed)  input. 
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'Virtual  source 


Flow  conditions  assumed  uniform 
across  plane  surfocc 


\ 


(b)  One-dimensional  approximation  of  flow  from  virtual  source. 

Rgure  17.  Virtual -source  representation  of  contouring  effect  of 

boundary  layer. 


Rgure  18.  Location  of  virtual  source. 
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Figure  19.  Radial  distribution  of  pressure  difference 
due  to  source- flow  ef^ts. 
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Figure  21.  Radial  distribution  of  total -pressure  ratio. 
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Figure  22.  Radial  distribution  of  static*  pressure  ratio. 
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Figure  23.  Surface  pressure  on  hemisphere- cylinder  model. 


FIgurt  29.  Osclllogroph  record  for  run  with  vorticoily 
trovtrving  probo. 


14 1 


r'  -»i 


2 


Run 

Sto 

x/d 

r|,/LQ 

RS|  d 

M|„ 

499 

A 

3.8 

0.0078 

2.I8X  10^ 

16.0 

498 

B 

5.1 

0.0080 

2.27  X  10^ 

15.9 

500 

C 

6.6 

0.0083 

2.19x10^ 

15.6 

Response  correction,  centered 


Figure  32.  Corrected  Rtot-pressure  profllet  for 
low  Reynolds  numbers. 
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5.1 

0.0080 

7.96  X  10^ 

17.6 

495 

C 

6.6 

0.0003 

8.02  X  10^ 

17.4 

Response  correction,  centered 


Figure  33.  Corrected  Pitot- pressure  profiles  for 
intermediate  Reynolds  numbers. 


145 


Run 

Sta 

x/d 

K 

517 
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0.0078 
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I.GI  X  10^ 

16.1 

514 
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6.6 

0.0083 

1.80x10^ 

15.9 
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0 

8.1 

0.0086 

1.84  X  10® 

15.8 

(Or-  - Response  correction,  centered 

I  -  Response  and  source-flow  corrections 


.8 


Figure  34.  Corrected  Pitot-pressure  profiles  for 
high  Reynolds  numbers. 
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0.0086 

2  .18  X  10^ 
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00089 

2.15x10^ 

16.8 
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II. 1 

0.0092 

2.19x10^ 

16.6 
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0.0096 

2.24  X  10^ 

16.4 

459 

H 

15.6 

0.0102 

2.04x10* 

15.6 

•  • 

RMponse 

correction,  cantered 

RespooM  and  sourc6*flow  corroctlons 


Figure  35.  Corrected  Pitot- pressure  profiies  for 
iow  Reynolds  numbers. 


Figure  36.  Schematic  diagram  of  experimental  procedure 
for  determining  transfer  function. 
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Figure  37  Shock-tube  arrangement  for  response  studies. 


(o^  Otcilloscopa  trace  of  amplifier  output. 


(b)  Oedllogroph  recorc  of  golvonometer  output 
(for  tm^  transducers). 


I  Figure  36.  Response  curves. 
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Rgure  42.  Comparison  of  tsmpsroture  profiles  in  uniform 
flow  and  source  flow. 
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